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ABSTRACT 

On ths Establishment and 
Evolution of Orbit-Orbit Resonances 

by 

Charles Pinnty Vodsr 

In ths solar systttD, thsrs exist several examples of 
gravitational resonance between two or more satellites or planets 
in which a special angle variable is observed to librate. Coldrelch 
has suggested that in the case of planetary satellites a tidally 
Induced torque acting on the satellites may have played an 
essential role in the establishment and subsequent evolution of the 
observed resonances. This proposal is thoroughly investigated as 
it applies to the three resonances among pairs of satellites of 
Saturn and is shown to be a plausible mechanic for their establish- 
ment but is less successful, in the Tital-Hyparlon case, in provid- 
ing a reasonable time scale for the damping of the amplitude of 
libretlon. 

The solution of the problem is reached In three etages. first, 
a theoretical description of transition is developed for a simple 
time dependent pendulum plus constant applied torque. The 
evolution of the system through the various phases Ci.«. pcsitiva 
rotation, negative rotation and libration) ia described in terns of 
ths motion of the extremes or "roots" of the momentum variable in 
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the complex plane. A transition phase is defined and equations of 
motion of these xoots are derived frost which a lowest order estimate 
of the probability for transition fro# a rotation into libratlon is 
obtained. 

Second, the bodv gravitational interaction is expanded and 
reduced to a one dimensional tins independent Hamiltonian which 
accurately describes the motior. of the resonance variable in the 
absence of tides - if the satellites' inclinations and eccentricities 
are relatively small and if the perturbations in tha semimajor 
axes during each phase of its evolution are also small. The effect 
of the tides is then introduced by redefining the orbital elements 
in such a way as to rscover the Hamiltonian formulation, the 
important difference being that it is now time dependent. 

The theoretical approach outlined for simple pendulum systems 
is ther. applied to eccentric ity dependent resonances. The dependence 
of the probability for transition into libratlon is obtained as a 
function of the mean eccentricity and the mechanism governing 
transition in various limits is discussed. The damping of the 
amplitude of libratlon as a function of the tidal charge in the 
orbital parameters (principally semimajor axis) is found via tha 
action integral. 

Finally, the theoretical model developed is then applied to 
the Saturn resonances and found to agree with the recent work of 
Allan, Greenberg and Sinclair. In addition, a proposal by R. <5. 
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Hipkin that the moon may have been trapped in an orblt-orbit 
resonance with another planet in the past is examined and found to 
be untenable. 
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1.1 INTRODUCTION 

In the solar system there are several instances in which the 
ratio of the mean motions of a pair of satsllites is very nearly a 
simple fraction. This kind of relationship is called a cosmien sua- 
bility. Some examples include the three satellite-satellite conr«n- 
surabillties of Saturn (Mimas-Tethys, 2ili Enceiadas-Dione, 2i *; 
Titan-Hyper ion, 4*3) and the 3t2 conmenvurabil ity of Neptune-Pluto. 
An equivalent statement of a commensurabll ity is the relation 
jX ♦ j # X» - 0, 

Where j and j * are integers and \ and X* are t' e respect We mean 
longitudes of the pair of satellites (or planets) . The obvious 
extsns Ion of the above relation to N bodies is 

N 

l j n X n - 0. Cl. l.D 

n*l 

Tht best known example of a cossnensur ability of three bodies 
involves the JI, JII, and Jill aatellitea of Jupiter throuoh the 
relation 


*JI " 3X JII + 2 *JIII 


Observed comaensurabil lties are not restricted to these orbit- 
orbit types. Another type involves the ratio of the orbital period 
to the rotational period of either the same or different bodies. 



«**• 
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Respective examples of this type are the spin-orbit interaction of 
Mercury (Goldreich and Peale, 1967) and the possible orbital 
coraaensur abilities of artificial satellites with the earth's slderlal 
day (Allan, 1967) . 

These special relationships would not be nearly as interesting 
if they did not have a physical basis for their existence. First of 
all, the coromensur ability relation is often not the physical 
var 'able which best describes the observations, examining the 
visual evidence more closely, we find that in many cases there exists 
a single resonance variable 4 which appears to librate about either 
mod'-) or mcd(?ir), and, for the two-body case, has the form 


linear function of perihelion (ffl) , 
and node (ft) of each body. 


The mechanism which mair tains this commensurability or resonanc# in 
every known instance involves a gravitational Interaction which is 
fairly well understood (Hagihara, 1972, pp328-52) . A study of the 
satellite-satellite interaction of the two-body resonance, after the 
two-tody gravitational potential has be^n expanded in terms of the 
orbital elements vi.lch describe the position of each body, reveals 
that 4 is that argument of a cosine function in the expansion which 
becaies very slowly varying for a nearly commensurate motion. 
Furthermore, this resonant term in the expansion often acts as a 
pendulum- like potential, being the dominant factor controlling the 
Ve r y long period motion of both satellites. This suggests that a 
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one-dimensional Hamiltonian might be derived as an approximation of 
this very long period behavior, and is the subject of chapter two. 

If the expansion (and Hamiltonian!) is valid for a range >t 4 which 
includes both the rotational and librational phases of the resonance 
variable, then, of course, 4 may execute either rotations or li- 
bra t ions, depending on the parameters of the system. Although the 
libratlon of the resonance variable can be explained in terms of the 
mutual gravi tational interactions between thr partners, it aeons 
unlikely that a state of libratlon could have existed sines the 
earliest stages of formation of the solar system. In other *rards, 
there should be some mechanism or mechanisms by which the partners 
evolved into their presently observed state. One interesting fact 
that Roy and Ovenden (1954) have shown Is that the high frequency 
of ccnooenrurabllities in the solar system cannot be assigned to a 
chance initial arrangement. 

There appears to be two basic solutions to the two questions* 

1) why so many commensur abilities? and 2) why so many libratlng 
resonance variables? One possibility is that a resonant or 
coimnensurable configuration is inherently more stable than a 
slightly off-resonant configuration. As Illustration, consider the 
Trojan asteroids which move in approximately the same orbit as 
Jupiter, clustered in two groups, 60® ahead and behind Jupiter 
(Brown and Shook, 1964, pp. 250-88). Their gravi tational interaction 
With Jupiter tends to maintain the one to one rommensurabil ity. 
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Imagine what wild happen to an asteroid with a nearly circular 
orbit which i* slightly larger or smaller than Jupiter's. Within a 
short time span that asteroid would make a close approach to Jupiter. 
If close ehough, the gravitational force of Jupiter could dominate 
that of the sun. and radically change the orbit of this asteroid, 
perhaps even removing it from the solar syston by changing its orbit 
frcm elliptic to hyperbolic, nr removing it through a collision. In 
this example, the lifetime of any asteroid not in resonance with 
Jupiter but having nearly the same orbit or a crossing orbit, would 
tend to be very short in comparison with the age of the solar system. 
Thus time, by Jupiter, endows the Trojans with a divine relationship.' 
A recent proposal by Ovonden (1972) , based on this idea of maximum 
stability, is that the high frequency of cammena'-jrabilities is a 
reflection of the evolution of the solar system towards a "Least 
Interaction Action" configuration, driven to its present state by 
purely conservative gravitational forces . 

The second possibility is that dissipative effects, which give 
rise to secular torques on the affected bodies, drive than towards 
a comcensurability with one or more other bodies, and that something, 
either in the nature of the f’seipetive mechanism or in its inter- 
action with the gravitations force, leads to transition into a 
libration of a particular resonance variable. Goldreich (1965) has 
suggested that the dissipative mechanism operating in satellite 
systems of the planets is the inelastic tidal respcr.se of Lae planet 
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to the pull of each satellite. Already a similar theory, when 
applied to the spin-orbit resonate of Mercury, has led to a 
satisfactory explanation of evolution and capture into libration 
(Goldreich and PeaLe, 1967). In this instance, capture is apparently 
caused by an asymmetry in the tidal torque acting on the spin of 
Mercury, as the velocity of the resonance variable 4 vanishes and 
then changes sign. An important difference in the orbit-orbit type 
of resonance is that the capture mechanism does not appear to depend 
in any important way on the details of the tidal interaction itself, 
as it does with the spin-orbit case. Recent numerical studies of 
Greenberg (1972) and Sinclair (1972) Indicate that the tidal 
mechanism does satisfactorily explain the Saturn resonances and that 
capture into libration is caused by the tidal torque acting through 
the gravitational interaction. The existence of a se/- 'ar torque of 
any significance acting on the satellites of the major planets has 
not been documented with corresponding visual evidence of a secular 
change in their orbital periods. TUs effect is apparently too 
small to be measurable at present. Still, an estimate of its 
magnitude has been inferred from the present sice of orbit of the 
innermost satellite (Goldreich. 1965). The best evidence for tidal 
friction comes from observations involving the period of our moon 
(Munk and Macdonald, i960, p. 190). In fact, the present rate of 
increase leads to something of a paradox in the age of the earth- 
moon system compared to the age of the earth, assuming a constant 
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dissipative mechanism* A novel proposal of R* Hipkin (in press) is 
that the moon may have been trapped in a cosnensurabillty with Venue 
at a fixed radius for a long enough period of time to resolve the 
time scale paradox. Unfortunate ly, the proposal does not appear to 
be feasible because of various factors discussed in section 4.2. 
Perhaps the solution is that. In the past, the dissipation function 
has been variable, as indicated by some paleontological evidence 
(Panne 11a . KacClintock and Thompson, 1968). 

At present, Coldrelch's tidal evolution hypothesis roust find its 
support through indirect evidence involving its consistency in 
explaining present observations and past history. Applied to each 
satellite-satellite resonance, it should lead to the results 1) that 
capture into the presently observed resonance is a reasonably 
probable event, and 2) that the time that this event took place was 
within the age of the solar system, given a reasonable estimate of 
the tidal torque. Allan (1969) has already shown for the Mimas- 
Tethys corasensurability, given tldally Induced torques acting on 
each resonance partner, that the evolution of t e orbital elements, 
including the amplitude of libration, could be followed backward in 

time to determine the initial values of the elements at the time of 

0 

capture, estimated to be about 2 * 10 years ago. Unfortunately, the 
approximations that Allan made for this case cannot be used for the 
remaining two cases. A major problem before us is to carry out a 
similar analysis of the other satellite-satellite resonances and 
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follow their evolution back in time. 

Attempting to understand how these three satellite-satellite 
resonances of Saturn evolved, in the context of Coldrelch's 
hypothesis, was certainly one of the major g^als of this thesis, 
although the first problem attempted was to determine the feasi- 
bility of Hipkin* s lunar resonance hypothesis. Lots of time passed 
before it was realized that each was governed by a one-dimensional 
Hamiltonian, although it turns out to be a much poorer approximation 
in Hipkin's lunar case. More time was spent determining how to best 
introduce the tidal torque into the Hamiltonian so as to preserve its 
canonical character. The resulting Hamiltonian is, of course, an 
explicit function of the time, and it is this explicit dependence 
which allows the system to evolve. Even more time elapsed before it 
was realized that no adequate analytical theory existed with which 
transition for even the simplest pendulum system was thoroughly 
explained. In the process, the scope of the thesis has broadened 
considerably and made it difficult to find some point to eno the 
affair and bring it to some conclusion. The exposition of this 
paper breaks down into three exerciscsi 1) development of a one- 
dimoncional Hamiltonian from the eatelllte-satellite interaction, 

2) development of transition theory for this pendulum-1 Ike 
Hamiltonian, and 3) applications. 

The development of an approximata description of the motion due 
to the resonance variable is a complex exercise using, for the most 
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part, well-known techniques of celestial mechanics and variational 
theory. The first step is to reduce the interaction to a one- 
dimensic .ltonian of the form 

H{x,*,t> - 1/2 lx + c(t)) 2 + blx,t)cos*. (1.1.3) 

Chapter tw> takes the spec it i«. rvaraple of a satellite- satellite 
gravitational interaction and outliner a procud .*re for expanding the 
interaction in terms of the orbital element*. The variational 
equations of motion of a canonical set of elements are also derived, 
A method for the elimination of the "short-period" terms in the 
interaction is sketched, along with a qualitative discussion of the 
approximations involved in reducing the system to one degree of 
freedom. The tidal interaction is further discussed and a method is 
proposed for introducing the tidal interaction into the tide-free 
Hamiltonian. In addition, a discussion of the other physical 
situations for which the above Hamiltonian is a good approximation 
of the motion is given. In chapter three, ^he analytic behavior of 
the Hamiltonian H(x,4,t) is discusser' in detail. First, the 
similarities to and differences from a simple pendulum are discussed 
for the tide-free case, and the possible notions of the system for 
different functional forms of b(x) are found using elementary 
analytical principles. Then the motion of the time-dependent system 
is obtained through an investigation of the motion of the "roots" of 
a polynomial in x. In the time- independent system, a pair of these 


9 

roots is exactly *he extremes of the motion of x, while in the time- 
dependent case, they at least bound the motion of x. This method of 
attack and its associated "picture" have rather wide application and 
can reduce many difficult probl«ms involving transitions between 
distinct phases to a tractable form. Capture criteria are 
specifically developed for two kinds of eccentricity -dependent 
resonance variables. 

The results are then applied in chapter four to the three 
satellite- satellite resonances of Saturn and to Hipkin's lunar- 
planetary resonance hypothesis. In the first exercise we find that 
the existence of a tidal ly Induced toroue does successfully explain 
the capture process. In fact, for two of the three examples 
discussed, the resonance variable automatically evolves into libra- 
tion. But the hypothesis is less successful in resolving the 
evolutionary time scale. The negative results of the second 
exercise uid have a sobering influence on those might over- 
estimate the importance of -e so nance phenomena . Fi: illy, the 
three-satellite coemensurability of Jupiter *n briefly discussed, 
and a probable history of its evolution is given. 

The material in chapter three concerning the • r criteria 
is different enough from other approaches to the problt require 

a lengthy introduction of its own, mainl' 'o understand nature 

of the approximations which will be used. First, t* ion of a 

simple pendulum with a time-dependent restoring t< i*. 
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investigated . Attention la focused on the "roots" . We discover that 
they are four in number and are points in the complex plane. One 
pair of these roots bounds the motion of the momentum variable, and 
each root ~an be uniquely labeled. Whether the angle variable $ 
executes positive rotations , 1 ibrations, or negative rotations is 
qualitatively determined by the relative position of the roots, along 
with the specification of these roots which bound the motion. This 
leads eventually to the precise definition of the trarsition pt*ase 
in terms of tne mutual motion of the roots a id of the momentum 
variable. To complement this picture, the flrat-order equations of 
motion of each root are derived, from which t* e analytical properties 
are deduced and a transition integral** is defined. Next, a constant 
torque term is added to this simple pendulum Harailtonirn and 
transformed to a new form very like (1.1.3) , except that the 
coefficient b is independent of x. Again the equations of motion 
of each root are derived, their soMon discussed, and a transition 
integral defined. These simpler system need to *,e clearly under- 
stood before we apply similar methods to the more complex c? 
certain amount of repetition is involved, but it is necessar' 
understand the scope of the theoretical approach taken. 
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Hie sign convention for the equations of lotion ±s the normal con- 
vention adopted in celestial mfw.anics a i Js useJ throughout this 
paper. Ordinarily, changing » • • sijn convention should result in 
replacing the ordinary Hrrilt -iiian by its negative counterpart to 

preserve the equations of moti- .. Thin should mean that th-.- kinetic 
2 

energy term (1/2 p ) should enter in equation (1.2.1a) with a minus 
sign. The above Hamiltonian is our subject of stJdy simply because 
of its similarity to the Hamiltonian develop d later from the orbit- 
orbit interaction. In that interaction, the equivalent kinetic 
energy term is negative-definite. 

Initially, say, the pendulum is executing positive rotations. 

If the coefficient b(t) slowly increases in magnitude with time, then 
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the velocity of the pendulum as It passes over the top will slowly 
decrease (Best, 1968) . Eventually, the pendulum will not have 
enough energy to pasB over the top, and thereafter will librate. An 
examination of the solution for the momentum variable p will suggest 
an eouivalent picture of the transition. If b't) were constant, then 
ve could fini a solution tor p in terms of the time by using the 
Hamiltonian to eliminat > 4. .he ~esulting integral solution ist 

r d P - * - v d-2-2) 

F ° /r (r )' 

where 

R(p> - b 2 - (h - 1/2 p 2 ) 2 . 

itie function R (p) is a quartic polynomial whose four roots are iven 

by 

Roots - i r 2iTb), (1 ,2.3) 

The motion of p is bounded by a pair of these roots, with p oscillat- 
ing back and forth between them with increasing time. Inspection of 
the equttion of motion for p reveals that these turning points of p 
(maxima anJ minima) occur when 

♦ - 2nw, or 4 » (2n + 1)*, 

n being an integer. In subsequent discussion, Trod(ff) and mod(2ir) 
shall d, signate 4 equal to 2n» and (2n + 1) v, rasper lively , If the 
pendulum is executing positive rotations (($) > 0 ) , then the negative 
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set of roots corresponds to the value of 4> at the top and bottom of 
its swine/. 

Therefore, the four roots can be completely labeled t y determin- 
ing the value of 4 (either aod(v) or ’2w) ) and the sign of (* 1 
at that root. Inspection of the equations of - tion reveals ttu-t 
the first set corresponds to 4 - mod (it) and the second set to 
♦ • rood (2n) , and chat they are completely specified by the set of 
labels p . p , p, and p_ . Factoring R (p) in terms of these 

Tfr IT* .jTT ▼ 

roots we have 

R(p) - 1/ 4 <P 2w+ - (P** - P^ (P*. * P>tP - P 2r J • (1.2.4a) 

and 

p >4 m /2(H + b) , p 2if± “ ±^2 (H - bj . b) 

Physically p is always real, forcing R(p) to be * 0. If all 
the rooti ire real, then the motion of p is bound between either 
p , p , or p . p_ . and corresponds to positive or negative 
rotation, respectively. The motion is bounded between the two r- 
roots or two 2ir-roots, or, o librates only if the opposite pair is 
complex. We shall adopt the convention that in the rotation pha«;e 
(all roots real) the ti-roots lie in‘erior to the 2w-roots, or b < 0. 
Diagrammetically, we can represent the three distinct states of the 
pendulum — positive rotation, negative rotation, and libration — by 
a graph of the relative locations of the roots ir, the complex plane 


I 
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FIGURE 1.2.1 FENDULUH 9TATE9 
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(•H figure* 1.2.1a, to, o) i Mao included \n than diagram* ara 
graph* of the function b(x) vtraua # In polar coordinate*. Th# 
pendulum llbrata* about 4 ■ nod(2*> If th* M < |b( and roUtea if 
M * ( b ( • Allowing b to too a function of tin* dnaa not qualitatively 
chary* hh# integral aolution given by (1.2.2) ucapt that both H and 
b ara tine-dependant, and tharafor* tha root* ara tlna-dapandant. 
Inataad of being fixad, tha roota now nova in tha c<mpu/ plana. 
Intuitively, wa aaa that paasaqa fren rotation to 1 Duration nuat 
lnvoiva tha notion of tha w-ronta toward tha origin and than out along 
tha imaginary axia. Thia pietnra of tha notion auggoata that w* 
look for aquation* of notion of tha roota thnaaalvaa, Since tha 
roota ara only function* of b(t) and n(x,f#t), tha notion of a 
giv-n root p r muat *ati**y 


U .2.5) 


Tha partial derivati-a* of p with raapact to b and N ara 
obtained fren tha roota thamaalvaa (l,2.4b) t vhila tha tin* 
derivative of H ia obtained fren ita partial tin* dorlvativa 

(1.2.id) i 


R-JMS-.- 


(1 «2,d) 







(1 


coe*) . 
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b) 


qualitative motion of each can ba easily determined, given tbs 

signs of and of the root. If |b(t)| increases with tine, and 

since we have adopted the convention that b(t) be negative, then the 

s-roota me vo toward the origin and the 2w-roota away from it. Thus 

the fluctuation in 4 caused by the pendulum force grows as 

transition to llbration is approached. One Interesting observation 

dp 

ia that the "-roots are stationary (-£—■ • 0) when 4 * med(w) or 
P - P^, But p itself is a minimum when it is at the root p^ and is 
stationary. If the root were not stationary whan p * p^, then it ie 
a simp. 1 * exercise to show that p would suffer an infinite acceler* 
at ion at this point. The analogous situation holds when p • p^. 

We could use (1.2.7a, b) to determine the flrat order (in 

flt 

secular behavior of the roots by Integrating them over one 
revolution and approximating the motion of x and 4 by replacing b(t) 
wich its mean value over the revolution. In this instance, the 
action integral represents a simpler method to obtain the secular 
motions. 

The above equations are uniquely useful in the transition phase. 
This phase will be defined by the conditions that 1) it starts at 
the Instant the "-roots coincide and become imaginary for some 
Initial values 2) it continues as the motion of p is 
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carried to the opposite 2*-root, and 3) it ends when p returns to 
the origin (the real part of tha p ff root) . This motion corresponds 
roughly to tha revolution in which 4 goee to zero, reverses sign as 
the pendulum moves backwards through the t tom and ana in ocas to 
aero near the top. xf Initially the pendulum executes positive 
rout ions and > o, then figure 1.2.2 represents the motions in 

tha complax p plans. 

The change In the w-roota during this transition phase can be 
obUlned by integrating the equation of motion. The result Is 

l/2(p 2 (f) - p 2 U>) - J** dt^| (1 ♦ cos*). (1.2.9} 

Notice that p 2 (lj * 0 from th# definition of transition phase. The 

integrand Is negative definite < 0), implying that pj(f) 1* 

negative or p^ff) Imaginary, Tha range of the Initial value 4^ la 

♦w * 4^ « 3s. The angle 4 then increases, passing through the 

bo t ton position of ths pendulum and finally raaehlng a maximum at 

4 ■ 3r. - 4* , where 44 1* a small positive angle of 

c c c 

0( W -V2 -M . 


4 then reverses sign, and the angle 4 decreases until it reaches a 
minimum at 4 f - * ♦ 64^(44^ > 0 and 0(|b|’ 3 ^ 2 

To simplify our problsm, lot's choose ~ to be constant, and 
demand that it be comparatively small. Also change the integration 
variable from t to 4* The result is 





FIGURE 1.2.2 TRANSITION PHASE 

T!ii» phase will be defined by C3#*Utions that 
1> it starts at the instant the * -root a coincide and became 
imaginary for some initial values 

Here we choose the minimum possible value of p^ to equal *t 

2) it continues a3 the motion of p is first towards the p 2lJ -root, 
where $ - 2s ard M2*-) > 0, and then back towards the point 
where ♦ vanishes at a maximum angle + c ) 

3) it continues as P reverses sign and che motion of p is carried 
to the opposite 2*-root where p again equals 2ff but M2w+) < 0; 

4) it ends when p returns to the origin (the real part of the 
p^-roo~) where p again vanishes at some angle near p « w. 





It should be emphasised that the above is a path-like integral in 
which the variation of p in going from p ^ to p f is determined by the 
transition phase diagram (fig. 1.2.2). During transition K(t) 
changes from one to a value slightly above it. Since at 

P • mod(w), k ( t) is very nearly a constant, no matter how long a 
time the pendulum spends near the top during transition. The 
contribution to the integral is small for p « * since the integrand 
vanishes. Accept for near the top, the motion of P is fast corpared 
to any change in H(t) or K(t), and the Integral in well-bounded. 

There is a stationary solution for p • 0, P • mod(r) or mod (2w) , 
and there exists the singular possibility of a sticking motion in 
which p slowly approaches the top and "sticks" there. Motions very 
near this singular event will have very long transition times, 
implying H(t) could change appreciably as p moves between the 
complex w-roots. The question is, how near? 

It turns out that thw set of motions which have a long 
transition time are restricted to an exponentially small set 
of 0(|b|” 3//2 of initial values of p i » To dmngtr&f this 

assertion, lot's determine the condition for which the change in 

H(t) in the time interval t * t * t. is of 0(H(t )). In addition, 

a d a 
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this particular calculation will be restricted to the event where p 
mves between the complex v-roota while * - 0+ and reverses siqn. 
The transition tire At - can be calculated wit easily by 

finding the analogous solution for + in terms of t. We find 


(1 ♦ K(t) cost) 


*c <1 ♦ Mt)cos*) 1/2 


- / tb <2H(t)) 1/2 d". 


( 1 . 2 . 10 ) 


The integral over + begins at an angle for which i is positive, 

continues to the angle ^ where ♦ vanishes and reverses sign, and 

ends as « moves back to the angle The integrand on the left 

hand side of (1.2.10) Is large only for angles very near * ■ mod(ir), 

and its value will tend to be Independent of the limits * and as 

a b 

long as they are not nearly equal to * . Therefore, the 

c 

calculation can be simplified by choosing ^ expending cos* 

about * « 3* - H and changing the limits of integration as 
follows^ 


♦ c - - 3w - 


The small differences, 6* c and 5*^, ere positive and 
sine- ^ Near the top of the pendulva h is very nearly equal 

to b(tj. This means that the function K(t) very nearly equals the 
valud one to 0(b 1 during transition, implying that the change 
in t(t) as the pendulum moves over tne top is of 0 < | b* 2, ~| 2 ) , A 
first order estimate can be obtained by choosing K(t) to be equal 


<■« - .i -r 


- ^WW****-' ■' 



the pendulum i J 
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th* s.*ll iibretion limit. Solvit? for i 


«4 * /2A4 - e“®» a r- 

c c 


1 3/2 _ 


-r — ii Ir 

3 T i 


a. 2.i5) 


Th« value of M4®) it related to p 2 tt) by (1. 2.4,6, 1.2.9, 

1.2.11) 

P*<4 > P * U e ) 

M V * - 2hT 5T a 2bTTT (1 * 2 - 16 > 

c c 

P*<4 r >. which depends on the Initial value 4^, can be 
approximated by 

1/2 P*<* > « || f 3w ' 4 *<= dt i c o »<V2»l l, +3 „ # 

* c dt n /!m«T 1 


<1.2. 17) 


Observe that (1.2.16) depends on the value of bU fl ) f not blt ft ) 

or b (t ) . The relationship between t and t can be found from an 
P c a 

inspection of (1.2.10) , in the light of the approximations so far 

invoked. Me see that each of the terms on the left hand side of 

1.2.10 are approximately equal. This implies that the integrals 

(obtaired from the right hand side) evaluated between t and t and 

c a 

and t also eaual. Explicitly 

/ tc (2H(t) l/2 dt i / t b(2M(t) 1/2 dt . 

*a *c 


The above integrals are approximated by the right hand side of 
a. 2.) 3). Given t^t^r * the important results arei 


<t e - t a > - (I 2 - '' ~ i-) 2/3 - 1)T m 1/2 t 
b(t c ) a 3/2 btt # ) . 

The next step is to obtain an approximate result for the 

integral in (1.2. 4. 7). Clearly, the only values of 4^ which 

correspond to a long transition time are near the value (3*i *■ 64®) . 

Therefore coM/2 can be expanded about 4 - 3ir. Prom (1,2.13), the 

transition time is rougly proportional to (ln| 6$ (t)+/2A + C 2 f(ti\ 2 ^ 2 , 

so that the integrand does tend to vanish if b(t) increases 

indefinitely. Since t change much more slowly than 4 for 4 not too 

near In - 6#®, we can approximate b(t) by its initial value bit^). 

The approximate solution for p 2 <4 ) is; 

u c 

P*<*c> 4 - ■ ‘^ 2 *"^ H3« - » 1 > 2 '«* C > 2 J- a. 2. 18) 


By contraction, 3w - 4, > 54 * implying that p. (4 ) is imaginary. 

1 C « c 

$4. can be eliminated using (1.2.15). Solving for (3* - ♦.), 
c * 


(3. - ♦j) “ ♦ l) l/ V a ( 


2/2* 


a. 2 . i9) 


Por largs values of a, the set of values of 4^ which lead to a long 
transition time is exponentially small compared to the full range 
of 4^. The above result agrees qualitatively with Best (1968) 
although he appears- to calculate a quite different parameter not 
nearly so well related to the initial conditions, fnyway, the 
important result is that the transition integral is well defined 


?r.'C*k£r- 4* ' 



FIGURE 1.2.3a 


Diagram of i 2 versus * for equation 1.2.20 where b ie constant 
(1) A value of r such that * ie, at some tine t*, a rotating 
variable. (2) A value of H for which * libratea. 
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except for initial values very close to the sticking notion. But 
the functional approximation used for the Integrand is only good to 

0( | b |l/2 ijbL,. 

Therefore, we can effectively ignore those ^ which 
have a long transition time, if all that is desired is a solution to 
the integral accurate to first order. Cf course, ir. any real 
physical system governed by (1.2.1), there are arbitrary 
fluctuations which would effectively eliminate the possibility of a 
sticking motion and inhibit transitions which take an exceptionally 
long time. 

A system more nearly related to the problem at >'*'nd is that jf 

dc 

a ^endulum subject to a constant applied torque, • The 
Hamiltonian in this case Is 

H(p,*,t) • 1/2 p 2 ♦ b(fcos* ♦ (1.2.20) 

We shall choose to be positive such that if the pendulum 
initially exacutes positive rotations it will be Blowed down by the 
torque, and i will eventually raverse sign. For the special case 
where b(t) • const., K(p,*,t) is a constant of the motion. Figure 
1.2.3a is a graph of i 2 versus $ for 1) a value of H such that ♦ is, 
at some time, a rotating variable, and 2) a value of H for which 
libratea. tn the first case, the graph reveals that the path of 
motion of i into zero is the same path it follows away from zero. 
Transition from rotation into libra tion cannot occur, except for 
the singular event of e sticking motion. Obviously, if capture is 
to occur, a non- time- symmetric term must be included in th|^ 




FIGURE I. 2. 3b 

* 2 

Diagram of ♦ versus 0 for a pendulum- Ilka system, subject to 

• 2 

a torque, which is a symmetric in time. ^ i» the kinetic energy Of 
a pendulum as it goes over the top for the last time, while is 
the maximum possible value. 6$ 2 is the kinetic energy after it has 
reversed direction and again approached the top. Capture occur - 

it «i 2 > i 2 . 
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Hamiltonian. Figure 1.2.3b shows how such a term breaks the time 

dc 

sysmetry. From the previous example we can deduce that if — is 

ct 

constant, |b(t) | must be an increasing function of the time for 
capture to occur (a 1 so see Sinclair, 1972), Also, we except that 
the criteria for capture wi-tl depend on the torque, the function b(t), 
its derivative, and on the initial conditions. Incident!/, in spin- 
orbit coupling, the time symmetry of (1.2.20) is broken by a 
velocity dependent torque (see Goldreich and Peale, 1966, and 
3.1.10-17). 

The above Hamiltonian lacks the simplicity necessary to express 
* as a function of p. Fortunately it can be transformed to a new 
Hamiltonian H(x,4,t) which has the requisite simplicity, defined ty 
(cf. 2.9.12) ( 

H(x,A,t) • 1/2 (x ♦ c(t)) 2 + bvt)co*$, (1.2.21a) 

where 

X *► c(t) • P, c(t) - ~ dt, 

*o 

and 


If - If - ’WW' 

b) 

II " *H“ " tx + c,t 

c) 

It • !if • H cos * + ft u + c> • 

d> 
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The vain* of c(t) can ba chosen such that * - *c in the absence of 
the pendulum force. The variable x then represents a fluctuation 'S&rv 
4 caused by the pendulum force. 

For c(t> and b(t) constant, the turning points in the motion of 
x occur at * - mod(w), mod(2x), and we can solve for x as a function 
of t f am was done foi the simple pendulum* 



to* 


U. 2. 22) 


R(X) - b 2 -(H - 1/2 (x ♦ C)V. 

Again, the quartic polynomial R(x) can be factored in tewus of its 
four roots, and the four roots uniquely labeled by the values of $ 
and sign(-$) for x equal to that root, Allowing c and b to be time 
dependent does not change this situation In the rotation phase, 
since the maximum and minimum of x still occur at modi*) or mod (2*), 
These roots are 


-c(t) i/2 (H + bit} ) 


(1,2,23a) 


x 2vt m * c(t> ±/2<R “ bCt)). b) 

The equations of motion for each of the roots are obtained in 
a manner analogous to that used earlier for the simple pendulum, 
except that the roots are dependent on three variables* c (t) , b(t). 


and The equations are 


dx 

dt 


ir± • 


dc, fJJSti 
dtSt^ + cCt) 
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(1,2,24a) 


dx . . dc, X * *2ifi 

dt 2wt - rrm 


The denominators are .qual to the value of (-}) evaluated at x - 
root. The equations differ from thos derived for the simple pendulum 
in the first term. Like those equations, th» motion of the roots 
is stationary whenever x » x^ + or x • *2r.t and ♦ * <> r ♦ ■ 2", 

respectively. As before, the fl-roots lie interior to the 2*-roots. 
Unlike those equations, however, the roots are not sjemoetric about 
the origin but about <-c (t) ) . Also, the motion of each pair of « 
and 2n roots is equal in magnitude and opposite in direction about 
this moving poi.t (-c(t)) (see 1,2.23), 

Let's choose x to move between and x ^ (positive 

rotation) If we Ignore the second cent in each equation, then 

moves toward the right and towards the left, implying that the 

fluctuation, Sx (Defi £x - x - x . ), qrows as the system 

max run 

approaches transition. The other pair of roots, besides separating, 
has a secular motion towards the left of O(^) (see figure 1.2.4). 

The transition phase begins, as before, when the two "-roots 
coincide at time and thereafter become complex. The equations of 
motion for the w-roots could be separated into their real and 
imaginary parts, but this procedure can be circumvented here by 
obeerving that x ff+ + c goes to zrro as H -> -b (t ) and then becomes 




FIGURE 1.2.4 POSITIVE ROTATION PHASE 
T*« arrows Indicate the relative velocity of each root. 



31 

imaginary (see 1.2.23a), Thua the real and imaginary parts obey the 
equations i 

x v Re x ♦ i I» x, (1,2.25a) 

Re x ff± » -c(t), b) 

1/2 fj- I*\ t - ^bt - «• ♦ <=*•*>• <=> 

2 

We should note that if x is complex, then to x is real and Bn x is 
positive definite. 

The related integral which determines the value of to x^ at 
time when x makes the second coincidence with Re x is 

1/2 Im 2 x (f ) - dt^-lx - Re x ) - 

n . ar n 

*i 

J tf dt|£<l ♦ cos*). U.2.26) 

ti dt 

From the previous example we expect that if to 2 * n (f) Is positive 
definite, then the n-rooto are still imaginary at time t f and * 
reverses sign, implying transition into llbration has occurred. But 
2 

if to x^(f) i* negative, it implies that the *-roots returned to the 
real axis before x reached Re x frnm the right, and the pendulum 
executes negative rotations. Therefore to x(f) - 0 corresponds to 
the sticking motion U •> 0-) which separates the transition into 
the libration phase from the transition into the negative rotation 
phase. Note that this occurs after the firet sign reversal of $■ in 
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which i -> 0+. This condition (that So x(f) vanish) is not completely 
accurate as shall presently be demonstrated. The important point to 
make here is that the above relation can still be used to find to 
lowest order the critical initial angle 4 ic which leads to this 
sticking motion. 

The description of the transition phase is more complex than 
that defined earlier (figure 1.2.2), T*. * important question to 

resolve is the relative motion of x with respect to Re for the 
period of time that the w -roots are complex. We should keep in mind 
that the major goal is to define the appropriate "transition 
integral" **iich can be approximated to first order in the small 
parameters. There a*-e two small parameters in this system of 
pendulum plus constant torque i 1) the first parameter is the ratio of 
the constant torque to the maximum pendulum torque and is small if 

!*»’ 1 !fl « tl . 2 .27a) 


2) the second parameter is the ratio of the relative change in b of 
0(b“ 1/2 to the initial value Mt^, and is small if 


b) 


2 

Since the equation of motion of Xm is already first order in 
these small parameters, we expect that the motion of x and 4 can be 
replaced by their zero -order motion in calculating the transition 
integral. Equivalently, the transition phase can be replaced by its 
lowest order approximation in defining the appropriate integral. 



The meaning of this statement shall became clearer as we proceed. 

The relative motion of x during transition can be discovered 
from the simpler case where b(t> is constant. First of all, the 
explicit time dependence of H(x,£,t) can be derived and is (from 
1 .2.21c, d ) t 

H(x,*,t) - const. - || 4 . U. 2. 28) 

The constant in this equation can be chose such that the argument 
of the radical in (1.2.23a) for the w-roots vanishes when 4 * 

The resulting equation for the ti - roots is* 

x vt - -•* /+2gf><* t ’ ♦>> a. 2. 29) 

and the tr-roota become complex for 4 > The minimum initial 

angle is mod(ir) and shall be chosen for this discussion to equal ». 

The initial angular velocity ^ must be 4 0| otherwise 4 would have 

previously vanished and reversed sign. Also, if $ vanishes at the 

moment the w-roots coincide, Lhen ^ must be equal to w. Unless 

- tr, the angle $ must increase for t > until 4 reaches a 

maximum 4 , at which time 4 vanishes, and 4 thereafter decreases, 
c 

Eventually $ returns to the value at a later time t^ and the 

ir-roots are thereafter real. From (fig. 1.2.3a), we find that when 

* 2 • 2 
4 returns to the value $ also returns to its initial value 

In the complex x-plane, 4 vanishes when x - Re x fl , and reverses 

sign as x moves to the right of Re x. After x moves to the right 



34 


of Re * decreases, tfe see that once * returns to the value 
x is still to the right of re since $ < 0 # and the ‘notion of sc is 
trapped between *+ and 2v+ and has entered the negative rotation 
phase. This sequence neglects the possibility of a sticking motion. 
Exactly how this notion of x would appear in the complex plane 
depends on th» relative magnitude and direction of the pendulum 


torque at ♦ » ^ as compared to the constant applied torque. The 
relative notion of x with respect to Re x^ for the time interval 
£ t * t^, when x^ is complex, is found from the equation of 
notion of e rh variable, x and Re x ff i 


dx 

dt 


-b sin*, 


d Re *» 
dt 


dc 

dt 


- (1.2,30) 


The two velocities (and the tw torques!) are equal for angles 

♦ (-s) and * (-2*) given by the relation 
o o 

»ln4 0 - b* 1 || . (1.2. 31) 

Claarly lb* 1 j|-| s 1 for th.se angles to exist. Turthenwre. the 
two stationary* solutions (* - constant) of the system correspond to 
these two angles. Figure 1,2.5a is a "physical” picture of these 
two stationary solutions, while figures 1.2.5b,c are their equivalent 
’■'presentations in the complex x-plane. Incidently , the distance 
of each of the roots from re X^ in figures 1.2.5b,c is obtained 
fran (1.2.23), and 


H - bcos* . 
o o 



b) equivalent picture of c) Equivalent picture of stable 

unstable stationary stationary solution. * o (2*) 

solution in complex x- is -2s such that cos* Q £ ♦!. 

plane, * (ir) is near s 
o 

such that cos* * -1. 

o 

FIGURE 1.2.5 


Location of stationary solutions of simple pendulum plus tormie. 
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Furthermore, fran U.2,30-311, t be relative notion of x is towards 
*2- lf 

dx(ti d Re x 

- . . . < 22 

dt dt # 

and towards Re otherwise. 

Careful consideration of these facts reveals that there are two 

qualitatively distinct "transition phases" involving the motion of 

the T-roots and the variable x. Figure 1.2.6a shows the relative 

notion of x during transition for the case where along 

with the equivalent picture for the real pendulum. This diagram Is 

very similar to the transition phase of a simple time-dependent 

pendulum without an applied torque (figure 1.2.2). On the other 

hand, if ♦. lies in the range » s ♦. s f , then figure 1.2.6b is a 
A l O 

picture of the motion during transition. The next question is, 
which of these diagrams is important? 

From (1.2.29), lm 2 x (♦) la most positive when ♦ reaches its 

M 

maximum value at $ - ♦ and t vanishes. Since H * bcos$ when $ 
c c 

2 

vanishes, X» x Tr ( ^ c J is al ®° 9 iven by 
lra 2 x ($ ) - 2|b| (cost + 1). 

H C C 

Thus ^ and ♦ are related by (1.2.29 ) t 

- V • |b|(co»* c ♦ 1). (1.2.32) 

The sticking motion i *> 0+ must correspond to a motion described 




a) Here the initial angle ^ lies in the range 
^ i 3ir. The diagram on the left is a 
description of transition in the corplex x- 
plane, while the diagram on the right is the 
equiv/ lent description of thw motion of a 


physical pendulw. 



b) In these diagrams, ^ lies in the range 

* * *i s * 0 (*> • 

FIGURE 1.2.6 

TRANSITION PHASE FOR PENDULUM PLUS TORQUE 


WHERE b - CONSTANT 
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by figure 1* 2.6b, in which the initial angle (sticking) lies in 

the range » s Otherwise, the pendulum has not just 

passed over the top prior to the coincidence of the r-roots. Recall 

that * Q M is the unstable equilibrium position. Since the right 

hand side of (1.2.32) is a maximum for 4 • 4 (w), it follows that 

c o 

the maximum amount of the w-roots can move off the real axle for 

tha set of transitions defined by figure 1 2.6b is given by the 

sticking motion. Since 4 Q * b~* Zo 2 x(st.) is of 0(lb _1 

and is effectively secona order in the small parameter associated 
dc 

*^ lth dt ‘ But th * ma * ijaus possible value of - is 2w, which 
corresponds to ^ motion given by figure 1.2.6a. The maximum of 
2m 2 x is therefore of 0(^> . 

dt 

This suggests that the following >roximations be eag>loyed to 
find the first order motion lin b~ X ||) of the s-rootst 1) neglect 
transitions imolving figure 1.2.6b for which ^ is in range 

* s s since this set of notions are of second order j 2) for 

_ db 

th * C “* e ST * °' the sticking motion where * -> 0- (and 

* ■ R « \ by the condition that Ira x - 0 when x - Re 3) neglect 
any effect connected with exceptionally lo.»g transition times. With 
these approximations, the description of the transition phase for 
the system of pendulum plus torque reduces to that for the simple 
pendulum (see figure 1.2.2). 

The next step is to find the first order approximation to each 
integral in (1.2.26). For the irst integrand, x + a equals ( 4 ). 
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If TT m constant, then 
at 

-/ tf «4f<* - *• *,) ■ *§!<♦, - V • (1.2.33) 

ti 

4^ lies in the range t s ^ S 3n while equals 3n to lowest older. 
2 

Also, the contribution to Im x(f) from the fi^st integral to 

be negative definite. 

The lowest order approximation to the second integral is 
identical to hat fo*»nd for the sample pendulum (1.2.17) and is 

i .-1/2 db t jd$cos 4/2\ 

' dt [i4^ * 4 * 3"; 3n * $ * *) 

-2|b|* 1/2 ||{3 ♦ sln^/2). (1.2.34) 

Thus the critical initial angle which separates transition into 
libration from escape into negative rotation satisfies the relation 

||(3» - » lc ) - 2| ♦ «iK4 it /2) 0 (1.2.3.';) 

to first order. Capture into librati' ' occurs f- -r ^ in the range 

dc 

$ ic , since the first Integra) propor c«onal to ~ is 

smaller than the second integral proportional to jr for in this 

range. On the other hand, if 4^ lies in the range $ ic * 4^* 3ir, 

then the pendulum has esep'^ into J\e negative rotation phase. If 

, as determined v y (1.2.35), is greater than ?.* . the implication 
ic 

is that the pe%iulum will inevitably enter the libration phase, 
independent of the initial conditions. 
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Often It happens that the a. ie > is unknown. In such an 
instance , a more valuable *-©o) would be a function which describes 
the probability that capture into libra “ion will occur, given 
probability densities for the initial argle of the system. The 
most physically reasonable assignment of pr'-b ,Ii.i y -s the following! 
If we Treasure the values of t an** ; tar from transition It ■* - «) , 
then for a fixed value y (-») the angle ♦ (-•) would be equally 
distributed ir r.e range ♦* * 4 (-•) s 4* ♦ 2s, where is 
arbitral/. Unfortunately, it is not clear how this stataaent 
t- ^slates in defining the probability associated with a given value 
of ^ at transition. For the special case b « constant, the 
translation is that is equally cistr ibuted in v nge 
' « * 3 ». 

Another parameter which ca. be assumed to be equally distributed 
in some closed range for this case is the value of i 2 as the 
pendulum moves over the top for the last time. Both Goldreich and 
Peele (1966) and Sinclair (1972) adopt this definition of 
probability in th - respective studies in which the time 'jyr - try 
is broken, in the former case by a velocity-dependent torque, and 
in the latter case by a time-dependent coefficient b(t). This 
second case exactly correspond- to th- example being discussed. 

Since it is always desirable to make contact with others 1 results, 

this definition of probability density will be adopted here. 

* 2 

Fro© figure 1.2.3b, :he value of + durinq the last passage of 



(x. ♦ cr. 
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At colt.cia*nca of tits x^ roots, J(r) v«ni shea and is thereafter 
imaginary. Therefore 


;2 


in I» ■ 
limit? 


• 3 


**" * -X* )** , (1.2.37J 

f i 

^♦d by letting 4 2 vanish. Thus, th* ma>lnun changt 

2 

and is obtain*! by evaluating In between th'i 
* 3« (where 4 *• Ov) i 


1/2 




U.2.3UI 


sous Is th-* xlnuin dscrtssi in Iis 2 *^ minus th* decrees* 


\n * 2 (, lc > " ‘/s - i/j u^lj" 


a.a.j») 


• 3T< 3 * * ♦ic' - W 1 ' 2 £<- 1 * 




Using (1.2. 3S) to slimlnsts th* explicit lependanc* on th* 

probability is 


Y 

c 


i -/a:! bf l jf)(|b |‘ i/3 ||>* 1 


(1.3.43) 


»ib 

? c is isro if ^ * o, while P c a pp lose he* its maximum value of v*.*ty 
whnn th* ratio of th* small parameters of ths system is of CU). 

Also, it should b* puir.t*d out that th* above formula does rot apply 
dc 

to th* s. tci/.I cats rrr ■ 0 sine* 0.2.35} is then invalid. Tn 
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section 3,1 v* shall i nd that th* estimate of apture probability 
for possUl* resonance associated with the Plmae-Tethye 
commeneurability agree with the numerical calculations of Sircleir. 

••fort concluding • hi ■ discussion, there ip one mrr* interesting 
feature wo shell investigate, related to transitlo \ fr>« th* 
positive rotation rnto the negstive rotation phase, Afsl it is 
described by the following! If we look a* the lysten of pendulum 
plus torque far from transition, at approximately equal us* inter- 
val* tafor* and after transition, ve observe that thnrs is « 
secular change In th* ***an value of x, o*. equivalently, in the 
iraan /slue of 4. That Is, if v* measure the mean value of | in the 
positive rotation phase far from transition and find that it is 
equal to, say, at time (where •> -•) , then the result 

of performing ■ ImUar measurement in the negetiv* rotation phaee 
at time t^ (where •> ♦••) shall b«i 

*•*" ' C, V' 


where Ax ie this secular change* Th* value of Am can be determined 

uetng 1) th* action Integral and 2) the condition (#<t ))„ - 

a pot, rot. 

-c(t # ), Th* time average of 4 at time is 

j<* * Odt * -<*> po ,. rot .- «<v- 




T< V 


<t .2.41) 


Applying the eacond condition, we find 
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rapidly lnoraaaaa for small cha nges in tha relative aaparation of 
tha v-roots. Therefore J can ba calculated at transition to deter* .ne 
tha aacular chans* in tha tims-dspsndsnt functions. At tiaa 

t i' H i * and 3 *■' 

V«.rot. ■ 0 ' /J* ** ’ 

« -4|b|t t )! X/a }\' loo. */2) - 2*0^) . (X.2.44) 

<* 3 po».rot. * • 8 ' blt i ) ' 1/2 - S,c( V * 

After tha psndvluro has mad* tha transition into negative 
rotation at time J can again ba calculated. 

’ IT *♦ ■ /?<-*>« • »"“*>' ,l - 2 - 4W 

Unless t f - t ± is exceptionally large, c(t 4 > * c(t f ), and 
b(t^ a b(t f i * Tha result is 

•Wot. ’ -Hvl” * 2,e( V' u ' 2 * 46) 

and J is noneero. Par from transition, tha naan value of x 

nag .rot. 

tends to Ax in tha negative rotation phase, and tha value of 

j corresponds to 2»Ax. Thus 

nag, rot, 

a* »•»•«*> 

" l 

For a simple pendulum, <-Ax) corresponds to tha delay in tha 
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evolution of <$> by the pendulum potential Cue to the constant applied 
torque. In the orbit-orbit interaction, the variable x is related to 
fluctuations in the orbital elements a,e, and I (see section 2.1). 

Thi* means that there ie a secular change in the orbital elements 
associated with passage through resonance, not entir ly connected to 
the tidal interaction. 

The purpose of this investigation has been to develop a 
description of transition along with foe* analytical tools vbich 
shall prove useful when applied to the Hamiltonian governing the 
orbit-orbit interaction. Before we proceed to discuss transition 
for the orbit-orbit case, the Hamiltonian 'hich approximates this 
interaction will be derived. 
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2.’ NEWTONIAN THEORY FOR PLANETARY SYSTOJS 

The development of a one-dimensional Hamiltonian consumes more 
pages than anyone in his right mind would want to read. This 
exercise la, for thr most part, a rehashing of old material to make 
it suit our own purposes, and breaks down into four stages. First 
Is the formulation of the many-body planetary problem in terms of a 
disturbing function, acting on each planet, *Mch ie distinct from 
the more common potential function. In the planetary problem, the 
principal interaction of the planets ic wi n the sun, and is the 
major factor in determining their orbits. The planet-planet {or 
satellite-satellite) interaction can be developed as a perturbation 
on the two-body pi, » ist-sun or satellite-planet orbit. Instead of 
considering porturbations on the coordinates, it is more useful to 
find the perturbations of the two-body "constants" of the motion. 

The second stage is the expansion of the disturbing function 
in terms of the Kepleiian elements of tha two-body orbit. It is 
amazing that the resulting expansion is of any use since it is so 
complex, but useful approximations can be more readily invoked with 
the disturbing function in this form. The terms in the expansion are 
classified, and tha relative importance of each clun is discussed. 

We determine the restrictions which must be imposed on these 
classes such that a singls hypothetical resonance variable will 
dominate the long-term behavior of the system. Also, a procedure is 



49 

outlined for the analytic elimination of the short-period ternu order 
by order in power* of a "email" expansion parameter. This serves 
two purposes. Pirat, it explicitly shows that the effect of such 
terms on the resonance is of secoj*! order. We see that these terms 
have a minor influence on the two-body resonances discussed later in 
14.1); on the other hand, in the lunar resonance problem u^acussed 
in (4."') the second-order mixing of short-period terms is important, 
because the sun can be e substantial indirect participant in the 
moon-earth-planet resonance. Second, we see that the terms in the 
expansion which look like a sum of pendulum-like potentials, all 
with different angles, is not entirely an accident of the expansion 
procedure — a 1 iteration of one of those arglos is implicitly 
possible. 

Once satisfied that the original interaction can be reduced to 
one which involves only one angle variable, we then show that the 
system of equations can be reduced, in most cases, to a single 
canonical set fx,+). The method used is similar to the first step 
in Delauney's solution to the sun-moon interaction {Brown 1960, p. 
140). The Hamiltonian derived is a constant of the motion in the 
absence of any dissipative tidil interactions. 

Our lact act is to interject the tJ .ea irto the tide-free 
Hamiltonian just developed. Essentjally this is accomplished by 
subtree' ing out the secula effect on the canonical variables. 
Pinally, we summarise the results in (2. 1C) Also, the relation 
between cur Hamiltonian and the si - and-orde^ equation of motion 
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derived independently by Allan and Sinclair is discussed, to Jrst;ll 
confidence that no serious flaws exist between beginriny and end. 

The beginning, of course, is Newton's laws of gravitation. 

The gravitational force between two bodies depends on their 
mass, shape, and the distance between them. Most celestial bodies 
approach sphericity, allowing the gravitational force to be 
approximated by that between two point masses. Given a net of n 
interacting point masses, the forces acting on the 1 ' th body are 
additive and individually derivable from a potential! 



-iV 1 * 


( 2 . 1 . 1 ) 


where operates on the coordinates of m^ and the potential V* is 



( 2 . 1 . 2 ) 


Por two bod ( ♦*! the path of motion each describes is a conic 
sect '.on, either an elllpi'', a hyperbola or a parabola. If the system 
is bound, the shape and siza of the ellipse are specified by Its 
eccentricity n and aemimajor axis a. Its orientation -J> resract 
to a reference frame is giver* jy the Euler anjies ft, l, w (cf. t. 
2.1.1). These symbols and their d finition* art peculiar to 
aetrwr'^any. ft is the "longitude of thj ascending * w is t. J 

"argument of pericenter , " whereas I is the "inclination". Although 
not a pt ileal angle, another frequently used "broken angle" ia 
Q, defined by 




FIGURE 2.1.2. ELLIPTIC VARIABLES 
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In most analytical work, M Is the most useful anomaly, since it is a 
linear function of 'he time In the absence of perturbations. The 
element t is the time of perihelion passage and is the sixth constant 
which fully specifies the two-body system. Instead of t, another 
choice for the sixth constant is , the epoch, which is defined b} 

nx * c - w (2.1.* 


The constant n is the "mean motion** and is related to the semimajor 
axis by 


2 3 

U • G(ia + m) - n a , 
o o 


(2.1.7) 


Which is recognized as Kepler's Third Law. 

The general problem of three interacting bodies Is stll*. unsolved. 
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In a planetary system, one body Is much more massive than any of the 
others, predominantly determining the path of motion of all the other 
bodies* This fact suggests that a zero-order solution for each 
secondary mass would be an ellipse whose focus is at the center of 
mass of the two-bxly system (primary and secondary) , the mutual 
interactions between the secondaries being ignored, The effect of 
the secondary interactions can be developed as perturbations on the 
zero-order ellipse, in which the six constants just described 
(a,e,l,c ,n,w) becoow variables. The method is known as "the 
variation of arbitrary constants" (Brouwer and Clemence, 1961a, pp. 
273-307). The six differential equations of the elements are first 
order in fime, compared to the three equations of the coordinates, 
which are seccnd order. 

The first step towards a solution la to expand the potential 
function (2.1.1) in terms of the orbital elements of the i'th 
disturbed body and the other n-2 disturbing bodies. To accomplish 
this goal, it is convenient to choose the primary mass as the 
coordinate origin and subtract the motion of the primary (nM caused 
by the disturbing mass. The result of such an operation is the 
equation of motion of the relative position vector r^t 


dr 


dt 


| ♦ V v* - + Sr*, 

2 O 


(2.1. da) 


where it is understood that r. ■ r .. The two functions V* and R* 

l io o 


are: 


m 

1-9 ± 


i)G 


b) 




1_ 

13 



c) 


R* is known as the disturbing function of the i*th body due to the 
action of the j'th body, and has the opposite sign from that 
commonly assigned a potential function. Its parts are called the 
Direct and Indirect terms respectively. In the absence of ary 
disturbing body excspt the primary, R 1 * 0 and (2.1.0a) describes the 
motion of the i*th mass with respect to the center of mass of the 
two-body system. 
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2.2 ANALYTICAL DEVELOPMENT OF THE DISTURBING FUNCTION 

The disturbing function acting on an; given secondary is the 
sun* of the individual distuzbing functions due to other secondaries. 
These stay be other satellites, planets, or even the sun itself if the 
given secondary be a Voon" of a planet. Our primary concern is to 
understand the two -body satellite-satellite resonances of Saturn in 
which the important perturbations involve a single two-body inter- 
action between the partners of the resonance. The critical 
development involves the expansion of A" 1 - | r ^ - r 2 ( 1 . The 
inverse separation, A 1 in terns of r^ # r^, and , the angle between 
r^ and r 2 , is 

a* 1 - l i-fcV P,(co«0>, (2.2.1) 

t«0 r > r > 1 

(Jackson, 1962, p. 62), where r > is the greater and r < is the lesser 
of r l and r 2 . If the orbits are coplanar, © is the difference in 
true logitudes of and r 2 

3 * L x - L 2 (2.2.2a) 

there 

L = f*v + a*f + &. b) 

if the orbits are also circular, then © reduces to © • - * 2 , 

where A is called the "mean longitude" and is defined by 
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where the.* superscript indicates the complex conjugate operation. 

The spherical harmonics can be related to the Euler angles {fi,I,f + w) 
through an explicit expansion in which the trigonometric relations 
between {$,$} and {f5,I,f ♦ w) are utilized (Xaula, 1966, pp. 30-35). 
Another approach is to use the group properties of the spherical 
harmonics under rotations (Iszak, 1964). The results are: 
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'W 01 " ** - 2p f a>^t ^ ■ *» - 2p^ co# <P’> 

tt ♦ ») I b' 

2 l pi (t - p) l 

The f i-^al expansion relates the true anomaly f and the distance 
r to the mean anomaly M. From equations (2.2.1), (2.2.4a) , and 
(2.2.5a) we see that the function which must be expanded is 
K—)"exp i (tf ) ). Hansen's coefficient? are defined by the 
relation: 

«• 

(j~) 8 exp 1 Itf) - J X* .(e) exp i (q + t)M. (2.2.8) 

<J - — ' 

The expansion of Harden' a coefficients in a power series in e is 

feirly complicated and can be found in Plunmer <1960, p. 44). Table 

(4.2.1) gives a few typical values to lowest order in e, while more 

extensive tabulations are published by Cayley (1961). The basic 

properties of (e) are 
s, t 

x? ,<«> - x' q «» <2.2.9.) 

X q (•> “V OcJ^I " b) 

s, w 

X 8,t <0) " V C> 

Using these expansions (2.2.1, 2.2.4, 2.2.5a, 2.2.6) we find 
the direct part of the disturbing function R* (2.1.8) to be 




(2.2.10) 


^ ~ ??1 ) } " ( ' 5 " 2p 2 ) V J ;- * + ®<-'W 

+q l T V q 2 T ‘2' 

This expansion can be reduced to a cosine series from the symmetry 
relations (2.2.7, 2.2.9). 

The expansion of the indirect part of R ia, from (2.1.8), 


R 21 * * V 2~I CO * 0 
r 2 


(2.1.11a) 


■ W 2“I l a ♦ ”) T *lap tt l ,r tap ( V 

e »,p#q i 2 


1 ♦ qi - 2p x . 1 + q 2 ZJ <€4) 

*1,1 - 2Pl 1 1,X -2,£ - 2p 2 ' 

e i ^imp 1 P2qiq2 ' 

and the sums are restricted to the terms where 

1 - lr m - (-1,0, lb Pi '*2 “ 


l2 .2.11b) 
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The qualitative dependence of R* for anall e arid I can be deduced 
from (2.1.6) and (2.1.8). Formally, R* is 


R 2 * l !i C lJ l " ® ’ 2p l' i I* * ® * 2pji e UJe hjl 

2 4 a> CX 1 1 a 2 2 °1 1 e 2 


COS ^top 1 P 2 q 1 q 2 (2.2.12) 

2 2 *< 

wher e C is a series in 1 , e and o • — and is formally presented 
by: 


o - i 


*1 


0 V 


- 0 


«r i «2 V2 - {2 - 2 - 13 ' 


The restriction on the sums in (2.2.13) are that l - 2p, m and q are 

constant (i.e. $ 0 . is fixed). The leading term is of 

“*PlPlP2<ii<42 

order n lmin , where l is the smallest value of l consistent with 
min 

the argument of the cosine and the range of values for t and p. 

2 2 

Terms in C which contain e or I are at least of order e or T 
smaller than the leading terms which only contain factors of a. 

The important point, which will be demonstrated later (2.6.10,11), 

is that any variation in C (due to a variation of the orbital 
2 2 

parameters) is of 0(e or I ) smaller than the corresponding 
variation of the leading factors. This implies that tW« 
variation of C with respect to X can be negler . r w.''.:r,ly. 


the secular term of the disturbing function is a polynomial series 
2 2 

in powers of a, I and e , and this term has no leading factors of 
e or I. 

Obviously, the expansion of is a very complicated function 
of the orbital parameters, and would not be vn* f useful if it were 
not the case that the perturbations produced by the disturbing 
body are small, and, furthermore, that the major portion of the 
variation of the factor multiplying the cosine function is 
determined by the leading factors of e and I, as demonstrated in 
(2. 6. 1.2). 

The expansion as outlined does have one defect. For the 
satellite-satellite interaction (or planet-planet) the ratio (^-) 
is not much different from one. Instead of expanding L 1 using 
(2.2.1), we can expand it directly in a cojino series in ■3. The 
result of this operation is 


r A’ 1 - 1/2 l bj (a)cos(j © ), 
j - -« L/ 


(2.2.14a) 


and the Laplace coefficient b ^ 2 (a) A * obtained from the integral 


°l/2 


r cos(j©)de 

ta > ‘ * 3 ^ / P l rf- 2, Vco.» ' - 1 - " 2< j.o 


b) 

The drawback with this approach is that S’ is a complicated 
function of the Euler angles, and a is a function of the mean 
anomalies. Subsequent expansions which reduce the disturbing 





function* to * function of th* element* U,e,If X,©,0) comprise a 
vary tadiou* exercise, a* anyone familiar with the topic wall know* 
<ct, PI une?ar, I960, pp. 133-48}. The resulting expression is formally 
identical to tha expansion given by (2.2.12,13) in that the leading 
factors of a and I in any given term can be factored out and the 
remaining sum of term* can be lumped together as the coefficient C. 

The only difference would be that C is not a power series in a, 2^ 
and e 2 , but a series in which the sum over a is expressed as 
functions of Laplace coefficients. 

The first approach outlined is usually reserved for the 
expansion cf the nan-spherical geopotential acting on an artificial 
satellite (cf. Allan, 1967). The ratio - enters in the expansion 
where R is the earth's mean radius and a is the semiraajor axis of 
the satellite, but the problem of convergence in power* of ~ is 
avoided because of the functional dependence on the appropriate 
multipole moment ir each term of the disturbing function. Rapid 
convergence of tne expansion is obtained, because the multipole 
moments decrease rapidly with the order of the moment. 

A mixture of the two procedures can be applied to the 
expansion of the disturbing function of a planet acting on the noon, 
with a definite reduction in xhe amount of labor usually required 
(Brown, 1960, p 252). Higher order Laplace coefficient* appear In 
the expansion and are defined by 


The series expansion of bj*(o) , Which can be used to numerically 
evaluate these coefficients, is given by 


b) 


b i t8) ■ 2 ° j TT .inj Pt *’* + 3 * l » a2> ’ 

2 

where the function F(s,s + j, j ♦ lr a ) it the hyper geometric series 

2. r rts + n)T(i ♦ 1 ♦ n) 2 n 

Hb.s ♦ J. j ♦ 1, a ) - l r<» * i ♦ n) ° 

n ■ 0 J 


(cf. Plussner, 1960, p. 158). 

The important point to emphasize is that the formal expansions 
of the satellite-satellite Interaction, the lunar -pi a notary inter- 
action, and the interaction of an artificial satellite with a non- 
rpherlcal geopotential lead to similar terms and that the behavior of 
these systems near or in a "resonance" where one of these terms 
dominates the disturbing function is essentially the same. 

The next series of operations involves the introduction of a 
Hamiltonian whose variables are functions of the orbital elements 
rather than of the coordinates and momenta. Eventually this 
Hamiltonian is reduced to a single degree of freedom and a constant 
of the motion for the tide-free case. To aid this process, let’s 
discuas the various types of terras in R* and the different types of 
resonance variables which can occur. 
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2.3 CLASSIFICATION SCHOIES 

The art of classification is to impose an order on a hodgepodge 
of material, be it physical objects or ideas or terms of a disturb* 
ing function, using criteria which suggest a meaningful way of 
thinking about that material. The first scheme presented classifies 
these terms using the relative oeriod corresponding to the cosine 
argument as the principal criterion. Its purpose is to suggest 
which kinds of terms can be handled using standard perturbatii a 
technioues and which cannot. Those terms which cannot be removed 
may be ignored if they satisfy the basic criterion that their 
corresponding coefficient of the cosine is much smaller than the 
coefficient of the given "resonant" term. Otherwise the system 
ca .not be approximated by a Hamiltonian having one degree of freedom. 
Xncidently, a "resonance term" is a term in the expansion of R whose 
cosine argument is nearly constant due to special values of the 
orbital angular parameters. This term is distinct from "secular" 
terms whose cosine arguments are Identically zero. We shall define 
a "resonance variable" as the argument of such a resonant term. 

The Hamiltonian contains an infinite number of terms which may 
or may not affect the resonance. A useful classification in treat- 
ing the "non- re sonant*, terms is the following: 

a) Short Period Terms. These are ttrras which explicitly 
contain > and have periods on the order of tor less than) 
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the period of disturbed or disturbing bodies. 

b) Long Period Terms. Those teems with arguments which do not 
ccntain X, but do contain a or are in this category. The 
largest terms in this class have two powers of e or I in the 
coefficient of the cosine. 

c) Secular Terms. The secular terras are those for which $50. 
Only even powers of the e and I occur as factors. If the 
Hamiltonian is reduced to terms of this class, then the 
action variables (defined later) become constants and the 
angle variables are linear functions of the time. 

There are still term# whose cosine argument is not some multiple of 
the resonant angle and which do not belong to any of the classes -Just 
defined but have long periods. They are of two types: 

a) Those terms which "almost" satisfy the comnsensurability 

condition 0, but which have different j 2 

than the resonance variable. If the integer pair is much 
larger than the i es^nant pair, then its coefficient will be 
much smaller than the resonant coefficient, and for this 
reason is ignored. If this is not the case, these terms 
must be directly compared with the term (or terms) containing 
the resonance variable. 

b) If an angle differs from the resonant angle by a function of 
the slow angle variables G or ft, then the period of that 
angle will be of the order of the slow variables near 
resonance. They can be neglected if either the coefficient 
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of such terms is such sms liar than the resonant tern or if 
the period of the slow variable is relatively fast compared 
to the libra t ion period of the resonance variable. 

There are several types of resonance angles which shall 
discussed in detail. Below is a list of these types, alonq with a 
list of specific examples. 

a) Purely Synodic: l - 2p. « - raj q. . « 0; 

X,s 1,2 

4 - - * 2 ). '’’he Trojan asteroids librate about the 

Lagrargian triangular points of Jupiter with a resonance 
angle of the above type (Brown and Shook, 1B64, ch. 9). 

Unfortunately, t^e perturbation expansion of the type 
developed n (2.1) cannot be used for the Trojan resonance 
problem. 

b) Simple e Type: t - 2p x 2 - m; q 2 - 0; 

4 - m(> i - > 2 > ♦ Most known resonances fall into this 

class. The Enceiadas-Dione (A^_ - 2A . + , ) , The 

En Di Oi 

Titan- Hyper ion (41^ - 3A Ti - 0^) , and the Neptune-Pluto 
(3A - 2 A - ft ) resonances are well-known examples. The 

P K r 

leading terms in these t^oes of resonances contain a factor 
of Actually, the class of observed e-types is 

restricted to t* e |q| ■ 1 case. 

c) Simple I type: - Oi i - 2p 2 « m: 

4 • (1 - 2p^ ) (f^ ♦ Mj) + - mA. . Ther' are no known 

examples of this type. 


d) Mixed 1 type: q x - q 2 • 0: 2p 2 > * 

(1 - 2p 2 ) (f 2 + w 2 ) + m(ft x - n 2 )}. The only example is the 
Mimas-Tethys resonance (2A Ki - 4A^ e + + ^ Te ) * 

Other types are the mixed e-type and mixed type whose definitions 
are obvious. No naturally occurring examples belong to either tyoe. 


M ■ * ■ 
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2.4 BMULtXm'S EQUATIONS 

The Hamiltonian which describes the motion of a disturbed body 
can be written as 

jj + y |, 4 p 

H - n ♦ At H - -1/2 v-v ♦ -S— 1 - , (2.4.1) 

co r 2*i 

where H q Is the Hamiltonian of the unperturbed two-body system of 
primary and secondary. (Note again that is the negative of its 
usual counterpart in ordinary mechanics. Of course, the sign has no 
significance except at tradition.) The canonical elements of H ue 
the conjugate momenta and position coordinates. H is not a con- 
stant of the notion since the coordinates of the disturbing body are 
contained explicitly in R. Since H q supposedly dominates R, the 
motion of the disturbed body can be described in terms of the 
variation of an instantaneous ellipse which determines its position 
and velocity. The simplest method for deriving a canonical set of 
conjugate action and angle variaLles is a net hod involving the 
Hamilton- Jacobi equation (of. Appendix A'. 

Several sets of canonical variables have been <?' **ived 
(Hagihara, 1970, pp. 526-555). Th<* set chosen for this discussion 
is known as the modified Delauney variables fL,I\Z ’ The 

angle variables have already been defined. The action variables are 
related to the Keplerlan elements by* 



and > is a linear function of time. 

The disturbing function coitainB many terms other than hue 
resonant :erm. it is important that tneir effect on the resortance 
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*2 


and the meaning of the appxt Umatlon where they an neglected be 
understood . acne times such *n approximaticn for certain com# in R 
i* unjustified, snd a procsdur* is outlined which determines th«ir 
offset on the resonance. 

In Me lunar-plaietary resonance problem, the tun ii ait 
iijro'tan* par tic ij^n t through a coupling or mixing of shcrt period 
term* 4 n H. Thorrfore, a method In discussed which successively 

• Umvnat** those tans* order by order in terms of a wall expansion 
parameter. The expansion parameter (or parameters) in th# e/telllte* 

• re.Utq c«i« is th* mass of tbs disturbing body y* (realty th# 

rat 1 j H-1 . The mass ratio ranges from 10" 3 to 10~ 8 for satellites 
o 

cf tie major plamrte. In the lunar-sclar problem, the expropriate 
j a -a Tie ter is the ratio of the mean motion* of the moon to Cite earth, 

r } -i 

, which i* of order 10 , Thu* the rate uf convergence of a 

U 

conventional perturbation expansion in the lunar theory i* 

'.ompsrat* vely alow. 


2.5 KLIM J RATIO* W THE #HC*T pmXOD TWK* 

The expansion t.ha disturbing function for two satellite*, 
developed in section ?.?# corUina an infinity of term* which were 
classified according to various rtchsmws in section 2.). If “he 
system is near a cossssniturabJll' o, the usual assumption is that tie 
secular term and the given raeonmee tern are relative) y more 
important than any of the othsrr. which occur in n . When restricted 
to th'a set of terms, the KamiH:»n .an eri be -reduced to a form very 
like zhet of a pnndul m (2. /.•> . :».!• r«»u< ;ted (Serb ) tcnUn allow* 

the »ngla vetiable to llbate and the eystem to uxhibit the property 
of r * eona rc* . Oi.e ml |ht question whether this resonance i heftomrfon 
really saint*. r erhaps it is Just a property of (1) -.he particular 
axpansion which conveniently ex) raises * as the sus o 1 pendulum-: iks 
potential# and (2) the given as'usption that one of tnem rfcerlnat*’* 
the motion over a long time sc a- a. We shall try to allay such doubts 
Rise, th* expended disturbing f. net ion is non-llnsar in M»# varies 
orbital elements. This non-lmtaril » in f should lead tc a 
"coupling* 1 of ter not directl, nvolv- *» the resc-nanc ♦. This 

indirect effect may have the rta'ie frequency as the resort*, t term 
and thus change its effective potential. Therefore, the principal 
subject dis missed in thle section will be this coupling effect.. TO 
explicitly display the coupling, a procedure similar to trovr.'s 
method will b« outlined, by which the short period toms in P ran be 
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•3 iminated gran t re Hamiltonian order by order In tanas of a * small " 
mansion parameter I Brown and Shook, 1964, oh. 6) . The coupling 
will te shown to be smaller by a factor of an eoqpansion parameter, 
although the parameter may bo different from that which occurs 
directly in the two-body disturbing function. Next, the conditions 
which determine whether the procedure can be applied to the long 
psriod terms are enumera ted . At 'his stage the Hamiltonian will 
hive been reduced to the *. ecu Aar terms, terms which contain a 
multiple of the resonance angle, and other turns which cannot be 
el Irena ted by a perturbation expansion because their cosine 
arguments are very slow functions of time* Finally, an explielt 
c tculaticn must be made to see whether the coefficient of the 
resonance ten la much larger titan the remaining terms. Zf the 
litter terms* effect on the partrer* of the resonance is comparatively 
small, then they can be ignored. 

To reduce the procedure to J ts essentials* the disturbing 
futtetion ot the first body will be momentarily restricted to a 
single two-body potential of the form (free 2.2.12) 

R(J, o J \w’> t) - n'| A jL (J.J*)cos4 i , (2.3.1) 

end other effects such as those due to the presence of other 
satellite*, othfjr planets, the sun, and the primary planet’s 
cblateness will be ignored. The subscript notation is here replaced 
by printed notation, where unpriced and primed variables refer to the 
first and terond bodies respectively. k**,w) Is shorthand for any 
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pair of son jugate action and anfUe variables of thm first body. The 
above disturbing function can be formally stfmrated Into two parts* 

R - R g ♦ R c . (5 .5.2) 

R shall Include the short per ic die terms and R all the remaining 

■ c 

terms. The Hamiltcnian is therefore (2.4.1) i 

H ♦ H o <L) ♦ R g ♦ R c . (2.5*3) 

Formally, R^ can bo written as 

R. - VZ A g (J,J')coa6 B# (2.5.4a) 

where 

♦ • j* ♦ j'X* ♦ (function of nodes and 

lonqltsde of pericenter}. b) 

Ths mean longitude X is also equal to (2.1.6,2.2.31 r 

X • nt + t . (2,5.5) 

Zt follows that the mean longitude is an explicit function of its 
conjugate action variable L. Ha make this point now to avoid 
confusion later on- 

The mext step is to maks a Kami) ton-Jaccbi transformation on the 
old Hamiltonian H(J, w) . with the demand that the new KAmlxnian 
generated, K(J,w), does not contain the short period terms to first 
order. The transformation f ran ths old set of action-angl^ 
variables can be accomplished with a generating fur ct ion 





75 

S(J,Wf J» V) defined by the relatione i 

J - — (J,vp J*,wM> w - <J,w* J’y), (2,5.6) 

a* w 

which satisfy the Itamil ton-0 aoobi equation 

UJ.Vf J',w') - H(J,vj J'.v ) ♦ ^ (J,*| J’,w') . (2.5.7) 

A simultaneous transformation will be made on the action-angle 
varzaolea of the second body and the cor x e spending Hamiltonian H* , 
where R* , H* , H' and S' are dnfinod by relations similar to 
(2. 5.1-7) . .Altiouqh H(J,w; J'w*) in the above equation is a function 
of the old variables belonging to the primed satellite, eventually 
the right hand side will be expanded In terms of the new variables 
of both .jatellitet. 

if "he perturbations due to the short period terms inH. 
are snal., then tie old Hamiltonian can be expanded in the 
difference betweer the old and new variables. A new generating 
function, S, can he defined which differs from 8 by the identity 
transformation 

G m Jw ♦ S . (2.5.8) 

Then, by (2.5.4), 5 is directly related bo the differences dJ, 6v 
defined by 

■tt « j * j «. - 22 

&w 


(2.5.9a) 
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Since J,w nave no <ocpllclt time dependence, § can be replaced by S in 
the H-J equation. Again, the tame transformation can be peu: formed or. 
S 1 with similar results. 

Neither in II nrr In H* does the time variable literally occur, 
yet both Hamilto.ilims are time dependent because cf the ajptatance of 
the "external variables* belonging to their respt* tlvc partr.ers. 
Explicitly, we find that the equation cf motion of H 1st 

dH dJ «K A dw 3H A dJ‘ )H A dw» 3H 

— (J.w, J .*» > “ 3t aj' + 5t 3w + 3t“ Tr 7 * dt~ Iw 7 


dJ 1 3H dv^ 

* dt 3J 1 " dt W # 


and since j- equals we have 


3H dJ' 3H dv» BH 

at " dt aj' dt aw* 


(2.5.1C) 


The implication of this last equation (2.5.10) is that the partial 
time derivative occurring in each H-J nqustion acts both on the 
explicit time dependence which may occur in the respoctiv* 
generating functions and on the "external variable*" which occur 
therein. Of course, B and 8* can be chosen such that the tine 
variable does not occur explicitly In either generating function. 


Therefore 


as 

at 



dw' as 
3t - 57 


(2.5. it) 
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7J 


To help us de tormina the beet c holes for 8, let's replace 
dw 1 

and jT“ their equivalents, using their respective equations <»f 
motion. We find: 
dH * 

H " - 5IT ft (4S) ' Vh * r# (2.5.12) 


at 


(«S) - 


J*,w' 


( 3<R C 




as 

aj* 


3{r; ^ *;> 
aw 1 


as 

3J' 


•) 


Ms expect that tr.e generating function 8 will be of O(u) since the 
short period tor, os which will be eliminated by 8 are of 0{»). 
Therefore, the s-jcund sum of tema on the left hand side of (2.5.12) 

2 *»; as 

i* of 0,v ), while the function - jp* i» of 0(p) a*.d may contain 

sheet- period t*ru. Therefore we shall demand that this function not 

appear in H. To see how this can be accomplished, expond H + R in 

o c 

a power series of the differences 6J, iw, 6j* and 4w' . In order to 
metch the variable dependence of £(J,v> J\w') and R(.J,wj 
expand P^ in 6w. The result of this e. pension is* 

K ♦ R ♦ I • H (E) ♦ 6R + ft (5 f 5| 3,5* > 4 6R 4 

O C 9 O to OS 

(J,5, J\V> + 6R g , (2.5.13) 

wh ure 61' , 6 r end $R are 
o c s 

jT »V 

«H o i B^l) - - jjp ♦ !/>-/•< «fc> ♦ Ofa» ) f (2.5.14s) 

dli 

fin * R (J,V, J' f w') - R (5,5; 5,5) ■ 

c - c 


«%• I - l T3T 41 ♦ 4w * ?jT to' + r^r «**• ♦ 0(ji )} 

J ,W*J (V 


b) 


JR 


4* 


- R - H <J,w/J\vM - I {— 1 fiw 4 0(» 3 l}. c> 

B w 

Tk 

The term 5 — 'L is similar to that found in (2.5.12), end the 

Hi 

factors (- ) are equal to n and n' r respectively (2.4.5). 

80 Car the equation for H, in terms of the new variables, J 01 


H(J,wj J'w'J- H (L) + 4H ♦ R (55, 5*5') ♦ 4R ♦ R (J,w r',w') 
o c c 0 

4 6% ♦ n||{J,5iJ',w # ) 4 

♦!*<«>■ 

If th* first order short period terms in H are to b« eliniutted to 

2 

0<uj, the following equation must vanish at least to 0 (i- )i 

n|^ 4 n'lyr e p'J * b (J,J' ) cos« a <5,w) - 0(p 2 ). (2. 5.15) 

TIM angle + (5,v') is a function of t)M new variables belonging to 
the unprimed partner and of the old variables belonging to the 
primed partner. Explicitly, 

• a (5 / w) « Jn ♦ J*n* (function of ,Ii* > (2.5.16) 

If E j « noun such that 

A W # J*> 

f(j,5/j',w’) ■ - ,# ^TT7?rTrrr ( w » w ' > (2.5.17) 

v in,n ) s 


*.***«*> 




w» >»« that tha difference between r. and r in the generating 
fur -t ion is « second order effect. In any w/.aa, the differences 6j, 
Cw can be derived in an ertireiy consistent manner from the generat- 
ing function giver by (2. £.17}. 7he only apparent drawback 1 m that 
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6% may contain terms of 0(p 2 ) like those in But if vt expanded 

the equations for (6J,dw) in terra of the new variables, we would 

,2 

find, for example, that auch products as 6 j — | i n the expanded 

«fl 2 

equation for dw would also lead to a second or<*er contribution to 

the secular and long period terms contained in (ft + R ) . 

o c 

The new Hamiltonian H is, to Otp ), 

H » fi + R ♦ 6 2 H * Sn + 6k ♦ *(?£), i2. 5. 18a) 

O O O C 8 0t 


where the second order remainders £R q and 6k* are given by f2.5.l4b,c), 
and « 2 H o , d(||) are 


« 2 h o - 1/2 7 ^r<|r > 2 + o<p 3 >. 


b) 


3 (R* ♦ R* ) 

c s 




Also, the unbarred variables are replaced by their barred counter- 
parts ~<i each of the reminders to obtain an expression for H 
accurate to 0(w 3 ). 

Some elementary conclusions can be drawn about the qualitative 

2 

nature of these remainders. Each of these terms is of 0(y' or jj'u). 

Obviously, the equivalent second order perturbations of th< primed 

2 

satellite due to the unpriaed are also of 0(y or p'p). It these 
contributions are to be negligible, then both » and p* »uit be 
relatively small compared to • GM^. We shall find in stetion 2.7 
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that In the ©rbi t-orMt problem, t«ou in » of C (p 3/ * 2 ) will be 
neglected In order to reduce the t'enUltenian to one dimensional form. 
Thus, the second order coupling la a ccsrjaxatively minor effect ir 
this problem. 

Koxt, the remainder iir/olvea pcocucta of short period term* 

in which can be rfexpressed with a single cosine factor* Ita 

argument will involvo the sum and difference of the angles which 

occur in a given jroduct. If the two angles art identical, the 

above will produce a second order contribution to the secular part of 

H. if the *jjt or difference bet we on two short period terms in a 

multiple of the j * so nance argument, "hen the above produces a 

aecrrd order contribution to the rooonance term in H. By a similar 

line of reasonlrg, t R c contains only short period terms, sines fa, 

1*, *tc,, are all short periodic, while R contains onlj long period 

c 

and resonance term*. On the other hanc, the remainders 4R # and 

are much like ‘V although tl*ese renainfari* contain many 

d3 

mere terms. Closer inspection of t (^*) reveals th*t only the 
products 


3R* 3P* 

* AL* an A 2. ASL 

IZrtF and 3J« fa* 


contain secular and long period terms* Therefore, the second order 
contribution to the secular and resonant terms in H are contained in 

3 2 H ^ 3R 3R* 3R 1 

% .* o J8.2 r r a 3S s 3S s 38 * t ... 

1/2 fi; 2 Si* + w J* w’ 3W 3J *** 3j* ^T r 5w r) * (2*. ,19) 


Wo can rhow that the contribution of the bracketed set of terms to 
the secular and resonant parts of H is equal to the following, to 
0<» 5 ' 2 >. 


c. cR s as a *i as SR i as . 

> l V + TF w m 5F 1ST* 


3R 


3* 


3R 


(2.5.20) 


n-zr *" * IF + + 


3w* 


V2, 


The first set of term* on each side U.e, ^ and T^jfa) agree since 

3*^ 

the only step has been to replace by ite equivalent fa (2.5*9). 

The next two sets of tern on each side can be shown to agree to 
0(v 5/2 ) by first replacing by (- ?~> . Here it is understood 

that the partial with respect to t acts only on the X and X* 
dependences occurring in ths co*;ine nreunent of each term. Taking 
the second jet of term* on the left hand side, let's rewrite it as 
followsi 


3V dS 3 *V 3s 
aF 3F " aT'^Tt gjT 

,3S f 33 


.38* 3 2 6 3 ,3S* as * 

3F5t - 7t <5F IF* . 


12.5.21) 


Ths product i|— r -jjr) can bs expressed as a cosine series. Some of 
these terms will presumably contain cosine arguments identically 
equal to sero (i.e, secular) or equal to the resonance variable. 

The partial with respect to t of the secular part is identically 
zero. The effect of the partial on the resonance term is to 
multiply it by 'In + J'n*). Pear resonance this frequence is of 
Oiv 1 ^ 2 ) , while the product of the disturbing functions is cl 0(t 2 ). 
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5/2 

Therefo v the contribution of this turn to the resonance is of Olji ), 

The remaining step in demonstrating the equivalence is to replace 

as* 3 2 s 3p * 

rjjv ay and by (~ |j~>. A* " hi * point, we can see that 

the lame operations can be opplied to the third o*»t of terms on each 
side in (2.5.20) to show that they alto agree. This second 
formulation is used in section 4.2 to find the second order 
contribution caused by the indirect action of the sun. 

There are a couple of drawbacks tc the procedure outlined, 
first, if we wish to determine the orbital elsnenta accurate to 
seccnd order, then a second H-J transformation mist be made on H, 
transforming the variables from 5, w, etc., to 5. w to eliminate the 
second order short period terms in H. Second, the secular and resonant 
terms are not restricted to the Hamiltonian H, but also occur in the 
differences tJ and 6w. This means that in comparing the observed. 

long period behavior of L wit its theoretical notion calculated to 

3S 

second order, we must include the second order contribution in 
Explicitly: 

2 

.Long period. - ,long period r _g_3 
'motion of V ‘ tami in: T ,y , M3X 

+ T33T ^ 

where L is found from solving tor the notion of H» We should also 
point out thi great power of the technique* first, the equations of 
taotion of the transformed variables are nUll canonical! Becond, 


j^fe***'****^ *»* 
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eliminated, and if R fw only contains resonant terms, then the coupl* period term Is 10* 1 •> 10~ 2 times smaller than the resonance term. 

lag of the long per iod terms will produce either secular or other 

long period terms. Unless the resonance variable is also like the 

long period tsrms being eliminated (j * j ' - 0) , the second order 

coupling will only affect the secular part of the Hamiltonian, 

It must toe emphasized that if the above procedure is to work, 
then each v must be 1) nonzero, and 2) large enough such that the 
perturbation expansion converges. Instead of u, the effective 
eagmnsion parameter for the long period terms is and it is this 
parameter which must be small. The secular tents in the two-body 
interaction usually do not satisfy thi 3 criterion* For example, the 
motions of the planetary perihelion ami nodes (which in some cases 
cannot be defined I) have periods 10 5 larger than the orbital period 
(Prouwer and Clmnence, 1961b, p. 46). Still, the procedure can be 
applied to at le.it t two of the satellite-satellite resonances of 
c*iturn. The reason is that the combined perturb® tins caused by the 
pJanet's 'blatencss, the sun, and the largest satellite Titan, lead 

to a relatively large motion in the per ic enter and node of the inner 
2 

satellite with periods 10 larger than their orbit periods (Jeffreys, 

1953). It should be pointed out that the least satisfactory case is 
tlie Titan-Hyporion resonance. The proyrade motion of the per ic enter 
of Titan is only about 0.5Vy®*r, while ~ts orbital period is 
approximately 15 days. The per ic enter motion of Hyperion is actually 
retrograde and caused entirely by the impressed resonance. In this 
case, we must appeal to the fact that the coefficient of the long 
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2.6 REDUCTION CP THE flDE-ndi HAMILTONIAN 
TO A CONSTANT OP THE MOTION 

have reached the stage where each Hamiltonian that describes 
the motion of the appropriate partner of the resonance can be 
explicitly reduced to one degree of freed- w. The two tide-free 
Hamiltonians which govern the motion of the first and second partners 
take the forms 

K(J,cr;J*,v/’) = A, (J f J' ) cob(*<P) 

and 

!!• (TV*' ;J,w) = f A;(J , ,J)coo{»0), 
where the angle 4 is 

$ - *A ♦ j'A’ + k*S + k*# 1 + if) + i'n*. 

(Again, ptrisnd and unprimed variables refer to the first and second 
partners respectively.) Recall that the two-body interaction must 
satisfy tbe criteria that it can be expanded in powers of a or 
Laplace coefficients, and that the hypothetical resonance dominates 
the motion over a long time scale. This means that all short period 
terms j r* 0, j* + 0) ar4 long period terms (j - j* * 0) are 
assumed negligible and that any terms in the disturbing function 
which nay have very long periods must ha vs relatively small 
coefficients compared to those for the variable <fr. For example, 
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there may exist terms which nearly satisfy the cosaensurability 
relation (1.1) for a different set of integers than occurs in 

These very long period terms are neg* ected in our approximation. 
Implicit in the expansion criteria is the fact that the two orbits 
cannot intersect or make very close approaches. Then the magnitude 
of the perturbation acting or. either satellite can be comparable to 
that of the primary. 

We can se* that the ocev ence of the angle variables in each 
Hamiltonian is restricted to the combination of angles which 
comprises $. The motion of the action variables is determined by the 
partial derivatives (acting cn the appropriate Hamiltonian) wit . 
reipect to their corresponding conjugate angle variables. Therefore 
their equations of motion must be proportional. Explicitly, the 
equations of motion of the first partner are 

dL . IB „ *\ r s\r.(n) 

*- * -kf *A,air.(nfl 

-It d« tVo 

dZ , . -if * 

d t dO. |,o 

Clearly, the variations of the action variables are proportional to 
each other. Therefore, a new variable x can be defined which 
simultaneously satisfies the relations 




Hamilton's equations reduces to 




-^(x.y.x-.y) 




where ♦ - y ♦ y*. 

Still# neither H nor H' is a constant of the potion. The 

reason is tnat the above Hamiltonians still contain an explicit time 

dependence because of the appearance of, say, the second partner's 

variables in the first partner's Hamiltonian. The next step is to 

establish some connection between the two sets of variables and 

find a new Hamiltonian which Lb a function of a single set of 

conjugate variables {x,$>. To do this, v* must c~.isider the form of 

the coefficients A^ and and determine how they are related. 

The secular part in each Hamiltonian corresponds to the y * 0 

term and includes H . Oh* coefficients A and A* can be formally 
o o o 

separated into two parts 

A (X,X* ) - 8(X) ♦ V (X,X* ) 

O 

A'(x',x) - s'(x') + 99 * (x* ,x) • 
o 

The function s(x) contains H , the secular terms due to the 
o 

oblatenr.ss of the primary, and the secular terms due to a host of 
other interactions, but not the secular part due to its resonance 
partner. The same holds for s'tx'). The mixed functions v(c#c*) 
and (v'(x',x) are the secular parts derived from interaction with 
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their respective partners, and are proportional to first order. This 
follow? erora the fact that the direct parts of the disturbing 
functions are proportional, and the indirect parts, to first order, 
have no secular t*erm <2. 1.8c). By inspection of 12. 1.8c), we find 
that v (x,x 1 ) is related to v'lx' ,x) by 

v’(x*,x) • J-vlXjX*). (2.6.5) 

Comparing similar coefficients A^(x,x*> and A^(x*,x), we find that 
they will also be proportional to first order, provided that the 
resonance variable $ is not contained in the indirect part of the 
disturbing function. From (2.2.10), the angles contained in the 
indirect part are restricted to the following sett 

q,q* * l* 1 - 2 P + <l)A - U - 2p* - q* )A* 

- q«- q'3' + (m - 1 + 2p)il 
-(m - 1 + ?p* )ix} , 

where the integers m,p,p’,q,q' are restricted to the valuesi 
td - (-1,0,1)! p ■ { 0,1 } ! p* - (0,l)i q,q* * 

• The Enceladas-Dione oomensurabillty, 

^En-Oi * " ^Di + '^Di » 

falls into the above class of resonance variables, dependent on The 
indirect part. This requires a slightly different approximation to 
reduce the Hamiltonian to one dimensional form. Except for cases 
such as this, however, the coefficients A^tx,* 1 ) and A^ t lx f ,x) are 
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proportional to first order. That is, 

A’ ?t (x*,x) = jj^A^x.x'). .C.6.6 

Since the coefficients are proportional, the equations of motion for 
x and x 1 must also be proportional t 

dx' m dx 
fi t " a 7TT 

Zn addition, the integration constant can be chosen such that x and 
x* are proportional. Explicitly, 

x'^rX. (2.6.7) 

To avoid confusion in taking partial derivatives, replace >' by ~ x 
in both H and H' . Making use of the above relationships, we see 
that the two Hamiltonians take the forms 

m 

H(x,x) - S(x) «► v(x,k) ♦ l A lx# X) cos (y ♦ y*) 
y-1 Y 


H* (x,x) 



The equation of motion for ♦(- y + y # ) is 


*1 


dt 


3s (x) _ ^ 3s lx) 

3x p ai 


— >vlx,£) 

ax 


- 7 (1“ + 2_)AlX,X)COS*. 

y-l a; 

The sum of the partial derivatives with respect to x and x cw be 
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replaced by the partial with respect to * if the bar on "x" is 
3 

dropped. This means that — acts on both x and x. The equations of 
motion for 4 and x are now derivable from the Hamiltonian H(x,$) 
given by 

!Ux f $) = 51 Af(x)cos(*(J) , (2.6.8) 

*0 


where A^x) - A^<x), y ? 0, 

and Aq(x) - s<x) + (x) ♦ v(x)r y - 0. 

R(x>t) is a constant of the motion, which can be verified by taking 
the ordinary time derivative of H(x,+) and relating that to a sum of 
partials. Explicitly, 

si-* $ Hr ■£■<>• <=- 6 - 9 > 

A slightly different approximation shall be made to establish a 
proportionality relationship between A^(x,x*) and A^Cx* ,x) for the 
Enr-Di case. From (2.2.10) and (2.2.11) we see that each part of 
each disturbing function contains the same leading factors, depending 
on the eccentricities and Inclinations of the two bodies. We have 
already noted that for small eccentricities or inclinations the 
fractional fluctuation in L(x) is of 0(e 2 or I 2 ) smaller than in 
rtx) or 2U), respectively. This means that these same leading 
factors dominate the total fluctuation in A^lx). 

We shall explicitly demonstrate this assertion for the y * 1 


term of a simple e-type resonance (is ♦ » jX + j *% + kfl) . From 
(2,2.12), A^x) 1st 



The function C(x) is a polynomial series in a, I 2 and e 2 , and terms 

2 2 2 2 
which contain fact -s of e or I are at least of order e or I 

smaller than the leading terms which contain factors of a. After 

expanding C(x) and e(x) to first order in x, we have 

coo - c( 0 > ♦ tiff ♦ i{rr 1^7 * * fr>* ♦ <h* 2 > 
k 2 M M , , IM „ 

e K (*) * H-kx - T )(—)]“ . .( 0 ) - k|k|e(0) f- 

O l O 

(2.6.11) 


The partials with respect to L, L* and r acting on each term in C 
reduce those terms which depend on L, L' or r by a factor of L, L* 
or T, respectively. Since the largest terms which depend cn r 

multiply a polynomial in a of the same order as the leading term, we 

C (0) 

see that each of the partials of C of 0( f ■■) . But in the 
|k| o 

expansion of e{0), the coefficient of t!* i.‘near term in x is of 

j( 

0 ( 0 (0 j_> the variation of C(x) is of 0(e 2 ) smaller than the 

r o |k| 

corresponding variation of e Cx) , implyinc that this leading 

factor predominately determines the variation of A 1 Cx) with x. It 

is also clear that the variation in x of the lesdino factors in A Cx) 

Y 

2 2 

determine the variation of A Cx) with x to Ole or I ) for the more 


■*"9 1 ' 


.'SeM* *5- - 
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cm ike in the previous case, the proportionality constant differs 
for each value of y. Since we have already presumed that e and I are 
email, the y • 2 tana in each Hamiltonian is of Ofe^e* fiUlj'U'lj 
smaller than the y ■ 1 term. If w r rict the Hamiltonians to the 
y m 1 term, then the action variables x’ and x are again proportional 
by the factor Still, this is not enough, because the nixed 

secular tsrm v(x,x*) in H has a different proportionality constant 
than the y - 1 term <cf. (2.6.5), (2.6.10)). If the aya^jm is to be 
described by a one dimensional Hamiltonian, then the mixed secular 
term must satisfy one of the following criteria* 1) The mixed ter r 
v(x,x*) is negligible compared to the unmixed secular terms and can 
be ignored. 2) The proportionality factor ' either very large 
or very small, such that the effect of one of the satellites on the 
other is negligible. Of course, this implies that the fluctuations 
in x or x* are negligible (if K^p is lurge or small, respectively) 
and the one dimensional form is ironed lately obtained. It should be 
pointed out that the Enceladas Dione case satisfies the first 
mentioned criterion* If this Is the case, then this Hamiltonian is 

H (*.-'!>) =«(.)< •■(>:) 4 AjUJcos-f, (2.6.13a) 

and the equations of motion a.*e 

dx 3H 
dt " 3* 


•A^x) sin*; 


b) 
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Ms shall use the above Hamiltonian as an approximation to the motion 
of all two-.'xxiy resonances, whether or not tne indirect pirt 
corr ibutes. This is not the end of approximations. Several more 
steps must' be taken before the above Hamiltonian is reduced to the 
desired form. 
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2.7 FURTHER APPROXIMATIONS AND STABILITY ANALYSIS 

Th( next strp is to shrv that the secular term in (2.6.13) can 

2 

be expanded in a rapidly decreasing polynomial series in x to 0 (x ), 

without increasing the order of appro.. station already established 

(roughly Mr—). Again, the secular part of the disturbing function 
L 

O 

can be formally separated into two distinct pieces: 

,.2 _ ,2 

{Scalar part) - £ (J) - ♦ {£ 5 %, 2 ♦ (*> • » 

The first two terms are the zero-order part of H. Recall t u 

o 

and are shorthand for \i Q + i» and m q ♦ M*, respective ./. The 
para^neter X equals one un.':S 8 the indirect part contaira a con- 
tribution to the resonance. In that case K is given by (2.6.10). 

The function s p (x) formally represents the sum of secular 
perturbations acting on both partners of the rescnar.ce. There is no 
problem expanding the first two terms, as long che fractional 

fluctuation in L9x) and L* (x) is small. Whether s (x) can be 

P 

ixparded is less certain. Recall that the secular part of the 

satellite-satellite disturbing function was a polynomial series in 
2 2 2 2 

a, 0 , e* , J and 1* , or, equivalently, a polynomial series .n L, 
r and 7, (cf. 2.2.13). Part of the motion of <f is derived from the. 
partial derivative with respect tc x of B p (x), and the effect of 
ti.ls derivative is to reduce any given term In & p tx) by a factor of 
L T or Z r ^.llc) . The reason that this is important is that if 
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ICO 


to contain teems with lectors of, say, as*, than thoaa 
tams would ptovidc relatively lor 9 # contributions to tha action of 
*3 or fl' if althar a or e* wars snail, respectively. Hons ol tha 
possible gravitational perturbations hava secular tanas of tha above 
type. A Jfc.-tor c? 00 * Is always associated with a Jong period tern 
whose coslro argument contains Ut - ft'). The conclusion we can draw 
is that > rj fr' is n-jt qualitatively different from the sero-order 
terms 4 ml can be oscle-ted without seriously aff toting ths notion of 
sither x or 4. The only sou ««tlon is ths case for angles 4 which do 
not conti in the mean longitude of either rartner. Tor these angles, 
the sero -order part Is independsrt of x (that is, J ■ j' ■ )) , and 
the contribution from e^(x) is uniquely important. Ignoring this 
except 10 : a 1 cast, we find the expansion of ths secular part to be r 
to a goo 5 approximation. 


Jar b'art* - A + 

o 


xA 


OK 


. x' A 
* T a oxx' 


(2.7. ?a> 


wh«e 2.1.7, 2.4.2) 

/. . >-F . , 

■7 . X* .i.' 


'i"o - l' n c' 


A . ♦ >-, V‘. 

OJ -< Al- j V * 


b) 

c) 

d! 


Comensuraulity is associated with the vanishing of A . Ths 

ox 

reeonanc ► ingle * car be constructed so that i and {-A ) ^rs 


positive as the system is pulled towards coassensuratillty by tie 
tidal accelerations. 

The secular part can ba expressed as a parfect square plu» a 
constant 


Utcular f irt\ * ^.\ oX)r (ic 



? 





The HemlltoriJerr then equals 


R(x t «) 


'■ A 7XX ( * ,1 


)‘4K 


.? 


tiCJL..*; (x)o>a4* 
o r.y 


Bines ths constant term <A —■) does not effect the er^tions cf 
oxx 

sot ion of x or 4 , It can oe Absorbed lrto »l wi*.ncut cha^i-ig ***«• 
equations, rven if this constant term wire tlr* difpervl'mt, it ttiU 
would not affect these ecuaUons (altht»uch it x'uH effect K' ar*i 
again this term cc be eliminated b> si htrsrti*.g it from H. Tha 
coefficient A^^ can be effectively foctortd out ty defining « new 
time variable* t related to t by 


W* 




whM. t >»• th« unit. mwl.nflth " 2 .»»« ti I l ). Th. «(uat. or. <t 
motion then t,ike the forts 

ik „ £i £] . Al.^.oLl) . 

nt dt e» ~ 


t vf.’O 


ii - Ji!f li - . 

dt it ** ox 


*’» * ’ - 4 


- ■ ft* ei " 
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0(6x 3 ) in the secular part is equivalent to neglecting terms cf 
0<u4k), which is of the same ordsr as the approx imat ion Urpoatf cn 
the function bU) . Note that the above Haaaiitorian is idantlcd to 
that of pendulum plus torque axcspt that tha parameter c hara is a 
constant and not a llnaar function of tha time, and tha coefficient 
b(x) depends on tha action, not tha time. Tha next major exercise 
is to Include ‘he affect of tha tides, and make tha parameter c tino 
dapandsnt . First, let's examine tha above Hamiltonian fnr Um 
existence and stability of libration canters. 

Tha exiatence and location of libration centers can be determined 
eithur analytically or graphicall-; fron tha phana-sraca curves 
generated !y the tide-free Hamiltonian K(x,4). A typical phasi-ipace 
curve in shown in fig. 2.7.1. Tha asyrvtetry stout tha 4-exir is 
caused by the momentum dtpendence of tha potential-! ike tana ii 
Also observe that in figure 2.7,1 there exist stable libration carters 
at ♦ equal to both even ard odd multiples of *, and la tha retult of 
the momentum dependence of the potential term. By the way, this 
kind of behavior is not general and depends on the form of th« 
momentum dependence and the vali e of the ra r <a**t*f* of the systecs. 

In the limit in which the moa tnt\sn dependence vannhes, ll f »s 
reduced to the Hamiltonian of a simple pt-ndulum. In this case, the 
libration centers are at either even or odd multiples of *, dop* nding 
on the sign of b. 

The analytical definition of a 1 Iteration center is as fo.'.Iow** 

Giver there ex Lets a point (x ,♦ ) such that fx ) is aero for all 

P P P P 



1 03 


i 



PIOD1U5 3.7.1 
PRASE SPACE OXA^XRAK 


1C4 


time, and if there exist closed curves about generated by 

H<x p + ♦ $♦> In sene region surrounding the test point, then 

the*- point is e libratlon center. Inspection of <r.C.), ?.$. 13b, c> 

shove that x^ can ba sat aqval to aero, and $ is «■ .Vtr * or :>*. 

P P 

Next expand the Hamiltonian in a amall ration about thu tost psint; 




i - r « ix :.i . *0(P) . 


This ipjadratic form defines a conic section (OaXley, 1749, p. 109), 

2 

and an ellipse ia defined by the condition * 0* For the 

s-llb ation canter, the above co ndition reduces to 

’'’" b XX> >°* U 7.9) 


For r. i. 2. -car tar tha aqulvalant condition la 

. i (l*b xt ) < 0. (2.7.10) 

If Ib^l l* laaa than on., than ba mu.t ba po.it l v. for a a- 
libratlon canter and negative for a 2s center. The sign of 1 
determines tha "nurtnal" libratlon center, as it did for tha t Iwple 
pendulum examples discussed in (1 2). If Jb^J 1s greater titan one, 
and Ms the opposite siqn of b, then the system can litrate m both 
its normal ft.Td inverted positions. Put if |b i hJB •«"* #1 9 n 
as b, then niethrr llbraticn cent* r U :*.ablo. TMi unusual 
phenomenon is caused by the momentum d«;w*r> donee of th« potential. 

We c.»n determine which resonance variables can exhibit either two 
stable libration centers or perhaps nom- by expirritly expardirg thn 


lu ’ | 4 ( ** < 41 .. 
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coefficient, using (2,2.12). Recall that xn e and I type resonances 
the function b(x) has loading factor* of e^ ard respectively. 

Comparing similar • and I type resonances for which *< * 1 (2.2.17, 

2.6.1) , we find that tha Kavlltoniani for aach art fnmally tha ra»a* 
Thorafora, b(x) no ad only bo expanded for a single e*t/pe resort nee. 
rM relative stability of tha resonances ;an ba dnt<cTntned as a 

if motion o:' k. Tha coefficients of tha expansion am (2.2.12, 2.4.4, 

2 . 6 . 1 ) ♦ 


K’C,* ,kl 


(3.7. 1U) 


b) 


0 } 


Crivativa* of C(x) with raapaet to x hav. boon MflntMl «inc. th«l* 

2 

contribution to b and b are of 0(« ) smaller than tha contribution 

X XX '•> 

from th** earemaion of a(x). Tha factor a ia cf 0(u >, ard 

> OXX O O 

therefore :he co*»fficiant is of 0(jj— e^l*! **), If tha naan 

value of tho ac can tricity, a , is snail enough, |b ’ can ba 
o XX 

greater tha.» n*. /or values of |k| < 4. Tha ooeffii lent b^ 
vanish#. i for | * 2 case, vhila b^ hat tha a am sign as b for 

! k| > 2. Onl*' r „*» |k| - 1 ease dost b^ hava tha opposite sign 

of b. Tn* jk| - 1 caaa la tharafera aapacially lntaraatlng, and 
tha pherwr.enu connected with two 1 ibration canters shall ba thoroughly 


investigated In (3.3). ~ 

At thia point we stake a brief digression to tha subject of tides 
and their secular affect on tha notions of satellites. A procedure 
is then outlined which introduces tha tidal affect into tha tide-free 
Hamiltonian juat developed. 







2.0 wracr or tide* 

•to understand the mechanics of tides, let's restrict ourselvert 
to a specific example, the earth-noon system. Because the moon Is a 
finite distance from the earth, there exists a gradient in the lunar 
force at the extended position of the earth. The earth is airoet 
spherical, but not perfectly rigid or elastic. Thus the earth's 
shape is distorted into a football-like object which attempts to 
follow the moon's apparsnt motion. The maximum response to the 
distorting force lags behind the applied force, as with any oscillat- 
ino system in which there xs frictioral lone and in which the forc- 
ing frequency is much less than the natural frequency. The magnitude 
of the energy dissipation may depend on the relative difference 
between the frequencies associated with the applied force end the 
natural frequencies of the effected body in a complicated wsy (Saule, 
1964). The quantitative si feet of c.dal friction on the earth's 
rotation is uncertain, but the qualitative effect Is well understood 
(Hunk and Macdonald, 19(0, p. 136) 

At present, t)» sa tj‘s MUt.cn Is much faster than the soon's 
orbital motion, an: the ‘f',t j.." - 1 caried ahead of the earth- 
moon axis. Sine.* .h« -.rat.Dnt xn K utr gravitational fome la 
eynuetric along thi nacth-soo'! axis, it attempts to realign the* 
distorted earth along that Mil. The lawltlng torque dsspins the 
earth and accelerates tho m>on in its or ait. The earth also raises 


e tide on the noon. Bit the noon's rotation is synchronous v.tn it f . 
orbital mean notion, and the "redial" tide raised on the noon < nly 
weekly affects Its orbit. This effect shall be discussed latts . 

The torque on the moon due to the tide raised on the earth no* 
only affects the mean motion, but the eccentricity and inclination as 
well. If the lunar orbit ia eccentric, then the torque on tie moor 
(and aarth) is stronger at perigee than at apogee, causing tie orbit 
to become more eccentric as It expands. This positive clang# in th»i 
eccentricity depends or the ear tn ' i rotation being faster than the 
moon 1 a orbital notion. The sign of this effect me/ be reverted onc*> 
synchronous rotation with the moon's motion is 'chiwed. Ttw 
"averaged" torque exerted on the moon tends to be ncmsl to <h<t 
ecliptic plena and tends to increase the normal consonant to the 
ecliptic plane of the lunar orbital momentiss, without changing the 
component in the plane. Therefore the lunar orbit is driven towards 
coincidence with the ecliptic plane, and the lu.v-r orbital inclina- 
tion le decreased. The results of a detailed calculation (K.r.~a, 

1964) show that the fractional rates ot change in a. e, and ! are of 
the same magnitude. Explicitly, 

S TV W e 3T * (_4) r X? * (?.fc.l) 

The numerical coefficients in (1.0.1) were derived under th# 
assumption that the Individual phase legs associated with each 
frequency in the tidal disturbing potential have th*» same vidua. 
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Pro* <2.4. 2) and <2.9.1 1 the secular o)*n<f*e in the notified Del aunty 
variable* L, r, and 2 art similarly proportional, obeying th* 
relation ' .t a and I ar* smaH) i 


i r» v. 




(Z.b.2) 


Wv» varulilt Z vt» promotional to tht* compont nt of orbital angular 
r#o«*ntum .n th* oclipb..c plan*, and fhi» cocap o net it tarda to ba 

-1 dZ T 

conserved. This la why A ■ it ernparatively small in (2.8.2). 

If the liucvr orbital station were fnan*r than th* earth's rotation, 
th* "footbnll" **?uLd te behind the enrth-moor lino, leading to a 
spin-up of th* earth*.’ rotation ar-d l sverslag th* secular effect* on 
tne lunar ul aments Jus.t discussed, r either c#se, the tendency is 
to synchroni** th* ro ation of th* earth with tie lunar motion. Ti* 
earth hast not achieve! that state, bit obviously tii# moon has, as 
have many of the sate.litet* of the other pla wtn (Goldreich and 
So ter, 196b). 

Cnee synchronous rotation is acl-eved, :hers It no longer an 
exchange or* orbital a td rotational nonenuia If, at ths aaas tine, 
the noon s orbit i« circular. Any tidal def>rmatior. of the mocm 
would re * *.n stationary with respect to a fixed set of axes vhich 
co-rotatx) with it. 5 -nee there if no tine-varying distortion with 
respect to a fixed mass alngcr.t of the moon, there ran be no 
friction v’lch would result in a dissipation of energy. 

But in neon is .to ar eccentric orbit. An observer on the 
esrth 3**3 a fjx«d po .nt cn the wcor. move teward and away from him. 


• .1 , 


no 

This in-and-out motion result# in a time varying dissipative ’radial* 
tide. Strictly speaking, this radial tide does not cause a torque 
along the direction of motion of the noon. Since there is no torqus 
acting on the :noon. orbital angular rxxnentjm ir conserved whi.e the 
orbital energy decrease!. This energy is given by <2.5.1) 

E or , Jlt t " 

Sins* E Qrbit decrease*, th* semimajor axis must decrees#. Ths 
relation between the tidal change in the xunar eccenti icity and that 
In the sunim/jor axle can be derived from the constancy of the lunar 
orbital nomentum 

•<r i/mU-o*). CM- 3) 



where T. refers to the radial tide component. The imp rtant pclnt 

to recognize la that the Induced change In • is opposite to that 

da 

caused by th* earth tide. “Jt haH *** ***** *«hctional form as ^ 

earth tide, except t).at it is multiplied by a factor of th* square 

of th* lunar eccentricity. All other thing# being equal , dt is 
da 

much smaller than ~ . Kaula <1964) indicates that the s^cuUr 
charge in e due to the radial tide is about two-thirdri tha* caused 
by the earth tide, at the present time. For other satellite systems, 



de 

apfears to lo larger than (Qoldreiah, 1963), But the 

support; ng .arguments Cor either case ore not absolutely convincing, 
especially since both depend on the relatLve dissipation in the 
primary and their satellites, of which llttla is krown. Fortunately, 
the quertioi of whether the tots l rate of chinge of the eccentricity 
is positive or negative is irrelevant for e-type resonances during 
transition ind after capture into llbration. The tidal tore je acting 
on the naan notion induces a change in the eccentricity threugh the 
ir.terventioi of the resonance. For an e-type case, this resonance 
induced change in e is much larger than the tidal ly induced change. 
This assertion will be supported later. The basi^ fact xs that the 
inelastic tidal response produces a secular change in the action 
variables, onlilce the purely gravitational f ertuebatior. which tends 
to cause onLy periodic variations in the action elements, ,f his 
point has not been proven, and the major theoretical eupporn for 
this conclusion is "Poisson's theorem on th*« invar i« bill ty o f the 
semina jor axes" and the Laplace-Lagrange theory of stcultr 
perturbations (Haglhara, 1972, pp. 164-05). 

have not discussed the secular chant/e in the angle 
variables tec* use the tidal terms are not qualitatively different 
from those produced by conservative gravitational forces and are 
much smaller. By far the largest change occurs indirectly through 
the implied secular change in L, 
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112 

which secularly changes H q (cf. 2.4,1), This suggests that a method 
of introducing the tidal effect on the tide-free Hamiltonian would 
be to define new action variables which are constants in the 
absence of any othsr forces except the tides. 



1X3 


2.9 INTRODUCTION OF THE TIDES INTO 
THE TT5E-FFE3 HAMILTONIAN 


In order to introduce tides into the Hamiltonian, let's return 
to the original Hamiltonian formulation before it wan reduced to one 
degree of freedom (2.4.3). First, the secular tidal chango in L, r 
and z shall be formally added to the part due to conservative 
^ravitatlonnl interactions already discussed. The result is 


IT 


(2.9.1a) 


dr 


: Kf ar T 
Ti + “IT 


b) 



C) 


where R is the tide-free Hamiltonian for the disturbed body (2.4.1), 
dt dr dZ 

and and are the rates of change of L, T and 2 due to 

tides alone. The equations of motion for the angle variable are 
well approximated by the tide-free equations of motion, because the 
secular tezms which ailse from the inelastic tidal interaction are 
much smaller than the secular terms due to conservative 
gravitational interaction!!. 

After adding the secular tidal change, the next step is to 
:n -reduce i new set <f action variables, L, f and 2, defined by the 
re La- ions j 






1 = X + J dt 

r = f + f at 

2 - 2 + r dt . 

The equations of motion of £, F and Z are 



df = JHf(J.w) 

41 a 

IT <5U 
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(2.9.2a) 

b) 

c) 

(2.9.3a) 

b) 

c) 


The elements L, r and 2 are constants in the absence of any 
gravitational perturbations other than the tides. Their explicit 
dependence on the tides has been eliminated, but the tide-free 
Hamiltonian depends on the unbarred action elements, rot on the 

barred variables. We can substitute the barred variables for the 

dL T 

unbarred variables in H^, but the secular tidal motions (-j^, etc.) 
of the elements are al so functions of the old elements, although they 
are approximately constant if the tidally induced changes in L, r 
and Z are relatively small. If this approx imation not possible, 

then the unbarred elements still can be successive l;, el. nr ted in 
the Hamiltonian by a process of iteration. 

There is another problem with the barred action variat'oB in 
that the equations of motion of the angle variables depend upon „ne 
partial darivatives with respect to the corresponding unbarred action 
elements. We can replace the unbarred partial with the barred 
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partial, with the restriction that tha partial darivativa does not 

dr„ T ar aa T 

operate on > , -^or -JJ 7 *® it occur* in K f * The equations of 
motion for each variable ar* now derivable from the same Hamiltonian. 
The procedure outlined in section 2. 1-2. 7 can be used to reduce the 
Haniltorien to one degree of freedom. But this one-dimensional 
Hamiltonian is no longer a constant of the motion since an explicit 
time dependence is retained through the tidal interaction introduced 
into the Hamiltonian. We cn see how this time dependence occurs by 
taking the time derivative at K f (J,w)i 


It L 


^ /dj 



Now replace J with J using (2.9.2), (2.9.3 ) 1 
d H t d J T r) H f ()Hf 

~Tt = ITT + "“5T* 

r "t 

Since J - ™^ dt always occurs in in association with J, wo can absorb 
the first sun of terms directly into the partial derivative of 
with respect to the time. 

The secular part of tin Hamiltonian has already been expanded 

in a power series in x. Tferefore, a simpler procedure for 

introducing tht tidal interaction into the tide-free Hamiltonian is 

to replace all the constants of integration (L t T ,Z , etc.) with 

dL 000 

their tidal counterparts {l. Q + /— } etc., in the one-dimensional 
tide-free Hamiltonian (2,7.11a), 

The next step is to determine just how these time-dependent 


llii 

terms alter the Hamiltonian. The approximate tidal acceleration 

at 

of the mean motion of a satellite, caused by the inelastic tide 
raised in the primary, is (Allan, 1969) t 






It is possible that t)» effective change in n during transition 
c*n be substantially increased if the ratio of the tidal notions of 
n and n* obeys approximately the same common sur ability latio found in 
the resonance variable $. Rather surprisingly, this occurs in the 
Mimas- Tethys and Enceladas-Dione corasensurabilitles, increasing the 
effective change in the seaimajor axis during transition by an or d«r 
of magnitude (section 3.1). This will not substantially affect the 
process of transition from rotation into lihratlon, but It can lead 
to an appreciable change in the long-term behavior of the resonance, 
especially in the amplitude of lihratlon (Allan, 1969) as a function 
of the age of the resonance. 

The effect of the time dependence of the "constants" is 
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greatest in the zero -order part of the Hamiltonian. Any tine 
dependence occurring in the coefficient of the potential term is of 
0(p) smaller, bet's return to the Hamiltonian K after the secular 
part had been expanded in a polynomial series in x. Replacing the 
constants of integration by their corresponding time-dependent te^nas, 
the coefficients of the expansion are (2.7.2) * 


A o = 

-*, 1 V 

i'<L. » j'JT. at)* *■ >:* i:V. *J 


A cx ‘ 

■ -3n 0 

b) 

A_ 

- • -'-v - . ** -ft - 

c; 

oxx 

-u)* * 


Earlier, these constants were manipulated to obtain a Hamiltonian 

which had the desired form. Me shall now show that these operations 

are still legitimate. 
dL T 

If , etc., axe replaced by their mean values, then each of 
the above coefficients depends only on the tlmr. Purely time- 
dependent terms can be added to the Kanlltoniar. without changing 
the equations of notion for either x or $. Again, absorbing the 
term (A - A 2 A* 1 } into the Hamiltonian is a legitimate operation. 

o OX 030C 

Introducing a raw time variable t, (2.7.4), related to t by 

t = Ja^ at, U.9.6) 

is quaivalent to introducing a nonlinear time and is also a 
legitimate operation (see 2. 7. 6-7). 
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Near a ccaraaensurabil it y , tha constant tarn la A is of OOi^^) 

ox 

which is abnormally snail. Therefore thr fractional change in A^ 
over an equal time is of 0(u ly ^) larger than tha change that occurs 
in A q or A^. Therefore, the coefficient c(t), (2.7.8b), is well 
apjr ox mated by 


* -*o»« n o ♦ i'K> 


dn dnj, 

,J 3r* # , 

o o 


(2,9.9) 


Observe that the tine derivative of c{t) la proportional to tha rat 
of the secular accelerations of the mean notions acting in the 
rtao ance variable. Explicitly, the tine derivative is - 


.2 dcCt) - , dn T 

* - * “ D TT”’ • 


dt 


it 


(2.9.10) 


The next question we shall resolve is the effect of the . 

dL. 

momentum dtw^enc*, contained in ■— , on these coefficients. The 

greater t dependence would seem to occur in the lowest order 

coefficient A^. But this coefficient has no effect on the equations 

of motion for either xor ( since the partial derivative with 

dL. 

respect to x cannot legitimately act on -gjv etc. Thus the function 

(A - A 2 a’ 1) can still be abrocbed in H without chinging tha 
o ox oxx 

equations of notion. The coefficient in which the momentum 

dependence is largest is therefore A , A depends explicitly on 

°* °* dn T 

the t.dal torques as they >ccur inside terms such as Jdt-r^(x,*t) • 
dn 

acpressirg (x,t) as a : unction of I*(x,t) w© find (2.9.4; i 
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where 

dL 

Ux,t) » jx + J-~ dt «► L- 0 
dn^ 

Instead of expanding in tems of x, we shall include the lowest 
order time dependence and expand In 

dL 

(jx + XT’ (t - t }). 

d. Q o 

To the fir at order terms in (2.9.9) we must add the term 

o 4 L 5 

The dependence on x in the coefficient means that the operation 
of redefining the time (2.9.6) is no longer explicitly valid. In 
this case, the dependence on x will be ignored so that the simplest 
form for the Hamiltonian can be retained. Besides, it can be 
demonstrated that in those cases in vhich this depenaerce plays an 
important role in the evolutionary history of a resonance, the 
contribution from this term is insignificant. 

The time dependence in the coefficient of the potential term 
plays a relatively minor role for e and I type resonances. This 
can be demonstrated by making a transformation of variables from 
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1* of 0(n a^) , or equivalently of 0 (L ). Unless the secular tidal 
o o o 

accelerations are commensurable, the first tens in (2.9.14) is of 
2 

0(e ) smaller than that associated with a(t). The same result applies 

to the element Z. We have already shown that the x dependence in 

1/2 

b(x) can be approximated to 0{y ) by the leading factors of T and 

Z, and the tidal chango in these elements has just been shown to be 
snail compared with that related to c(x,t). We shall now show that 
the above Hamiltonian equation agrees reasonably well with the 
second order equation of motion derived by Sinclair and alsc ty Allar. 

If the fractional fluctuation in T and Z are small, then the 
coefficient of the potential term b(x) can be expanded to first 
order in x. The equations of motion for x and $ are 

— - -{b + xb )sin$, (2.9.1Saj 

dt ° * 

^ • *x - c - b cos$. b) 

dt x 

Taking the time derivative of — and substituting for — , we have: 

dt dt 

^ - (b . xb )»in* - .. ain*2t - — (t,x> - 0. (2.9.16) 

at 2 ° * * dt dt 

To 0<v : ' 2 ), St equals (-x - c(t,x)), and the two terms In (2.9.16) 

dt • dc 

dependent on x and tend to cancel. The function — Includes the 

dt 

nearly constant torque tern (2.9.10) plus a term <#hich is 
suspiciously like a J term. If we eliminate x in favor of $, change 
the time variable from t to t (2.9.8) , and replace b Q by etc.. 
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the result is: 

“ ° 

at at 

12.9.17) 

where the factor P in the dissipative term is defined by 

e - 16 jV 1 ^ + SH. 16j'^I»'~* 37^ " 0. (2.9.18) 

J o dt JJ ' o dt 

o 

If is small compared to b , then the tine dependence of 

dt ° 

the torque will have little effect in the libra tional phase. In 

addition, the secular change in c(t,x) is very nearly c(t>. The 

validity of the expansion is restricted to a time interval for 

which c(t)b is s’lall compared to b during transition. The first 
x o 

three terms (with c(t,x) replaced by c (t) ) constitute the equation 
Sinclair (1972) derived in his njmerical calculation of transition 
probabilities for the Mimes-Tethys resonance. This restricted 
equation is identic*, to the second simple pendulum example treated 
in section (1.2). Allan neglected to consider the effect of the c(t) 
term on the problem of capture into libration, although he did 
implicitly include it and the i term in his theory on the evolution 
of the Mi-Te system after transition. 
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2.10 SUMMARY 

After a long sequence of approximations, we find that the 
satellite-satellite interaction can be approximated by the one- 
dimensional Hamiltonian 

H(x,4,t) - 1/2 Cx ♦ c<c t)) 2 ♦ b(x)coa$, (2.10.1) 


The equations of motion for x and t are 


dx 

3T 


iH 



and it is understood that — does not act on a(x,t). The explicit 
relations of each of the variables in H to the parameters defining 


each orbit follow. The angle variable is*. 


= 3 A + 3**' + k«+ + iCWi'A* . (2.10.2) 

The most powerful restriction on the integers j, j', etc., derives 
from the fact that the interaction is independent of *'.e 
coordinate syst Since each of the above angles is measured from 
a cormon reference, the sum of the integers must be zero. See 
(2.2.10) for further restrictions which appjy in the two-body case. 

Next, the momentum x is related to the elements e, I (e and ( 
small) and a by (2.6.1) and (2.9.2): 

L(x,t) = t t 

r(>:,t) * -ie°L 0 = kx+: 0+J ^lT dt , 


(-.10.3) 
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Z(x,t) 2 -iI 2 Lo 5 ix*Z 0 +j|jftdt, 

L*{x,t) etc. 

Thus, the variable x is proportional to the fluctuation induced in 
the elements a r e, and I, etc.# by the term in the disturbing 

function which containr the angle $. Observe that the fractional 
2 2 

fluctuation in L is of 0(e or I ) smaller than that of T or Z. This 
implies that the variation of b(x) with x is principally determined 
by its leading factors of I*(x) and 2(x) - if the inclination and 
eccentricity are small. 

The function b (x) is related to the ^-dependent term in the 
disturb. ng function by (2. 7. Sc)# i.e.i 

Mx) = A~f ;> A 1 (x), (2.10.4) 


where (A, tx)cos*> is the corresponding term in R (2.2.10,11). A 

1 oxx 

is the third coefficient in the Taylor expansion of the secular part 
of K and is approximately (2.7. 2d), 



Km 




( 2 . 10 . 5 ) 



^he parameter K equals one if the indirect part of the disturbing 
junction does not contribute to the resonance term. Otherwise it 


i# -iven by (2.6.10) . 


The function c(x#t) is proportional to the acceleration of the 




notion of 4 and is (2 9.*, 11)* 
c(x,t) * cvt) ♦ pfxdt 


( 2 . 10 . 6 ) 
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° m *L a - ^sr 4 *'3?» (t - V' 

o o 


♦ Si 16V 2 L* _1 


0 4 «3V-i 4 4 7 


The x -dependent te/OT in (2.10.6) results from the fact that the 
tidal acceleration decreases as a function of the planet-satellite 
separation. 

Finally, t is related to ordinary time t by 


The most important restrictions and approximations imposed in 
deriving (2.10.1) are the following, 1) The disturbing function 
can be expanded in *.ne ratio of the semimajor axes, and is valid in 
both the librational and rotational pnases. Therefore, the above 
Hamiltonian (2..V0.1) does not apply to a two-body resonance of the 
synodic type whir’ is restricted to the lil conoensurability . 2) 

The fractional fluctuations in the senimajor axes (or in L and L*) 
caused by the perturbations are small. This allows us to exp».v3 
the secular part of H in a Taylor series in x But there are no 
res trie tiors on the magnitude of the fraction f x actuations in 
either * or I. 3) There exist no terms in the disturbing function 
(such as the thort period terras) which cannot be removed or ignored 
because of- their high frequencies or comparatively small coefficients. 
4) Another condition is necessary in some cases (2.9), namely, that 
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the orbital inclinations and eot wr.tr lc it lea ba relatively srall, Tha 

reason ia that we impose the raatr lotion that fractional fluctuation 

In LCx) ba amall compared to * . ; nf f(x) or Z(x). Cna conaeouence 

is fhat tha fluctuation in A. .x, io 0 , governed by lta 

1 L o 

leading *<ictori of ; and a. This greatly simplifies tha functlcnal 

fora of btx). Tha raatriotion to snail Inclinations and 

sccantr icitlas car ba iiftm*, but it means that t .a tarns in bfc) 

o' C<~) must ba npUruy determined to find tha ooafflelant U . 

L o 31 

This Is a non- trivial ..arc i as in tha firsb-ordar expansion of 

t may ba possible to Improve tha ordar :f approximation to 

0(uS if tha indirect part of tha disturbing function do as not 

contain a 4-dependent tars. Also, i Jia mass and/or angular 

moms n turn of cna partnsr la vary touch largar than that of tha serond, 

than a bat tar approximation may bo obtained by naglaotinq tha B)t..on 

of tha f .rat partner dua to tha second, Jnfortuna tely, tha ord»r of 

approximation may not ba sat by tha axpanslon perimeter, but by 

those vary long par iod srma which had to i* ignotsd bacauaa thiy 

could not b# removed by any tochnlquo such a a Brovn'a mat hod . 

Although wa hava usad the apacific example o t tha satelllta- 

»«tel 11 te gravitational interaction, tha asms mtMi can ba 

amp Icy ad on gravitational raaonancaa involving aosa than two 

satellites. Of course, tha manybody raaoranca muat satisfy tha 

same general raatr let ions outlined for tha t^-body case. In 

addition, gravitational raaonancaa with the geopotnntial can b« 

similarly reduced to the one dimensional Hamiltonian form if t)*y 


t a ****' ‘* e 1 


also satisfy the relevant or 1 tar la. 


120 



119 


3.1 JllTRODUCTIOr'l 

*cw that the necessary preparation has bean completed, va are 
raady to turn to the problem o t transition tn orbit-orbit cora»*naur- 
abilities. The cueorelical basis of our dissuasion of transition of 
penduUn-lik* ayatemt subject to torques va* l^id o^rliar with a full 
investigation of two simple pendulum examples. Recall that th« 
corresponding Hamiltonian oj th# first (X> and second til) example* 
were (1.2.1, 1.2. 2^1 r 

lx f. r (?*$*6) « ip* + b( t)cooO, 

ir. H(x,*,t) - 1/2 (a * c<t)) 2 ♦ b«t)cos*. 

From a study of I, a picture of transition was developed which 
involved the notion of th# turning points or roota of (R(p) - 0) 
obtain *3 from the Integral i-olution for Explicitly, 

" R(p' *> 

The ei mi f leant r#*»ult* worm 1) The roota were labeled in t’ i 
rotation phase according to tha value* of t at that root (either 
mod (* ) or oodW*)) and to the sign of p(- -♦) . 2) The relative 

pos^ttn- of tha s# roota in zba complex plane oualitatlveiy 
determined whether the system executed rotation* or iLbrations 
(fig.I 2.Xa-c) . 3J Transition involved the motion of tho interior 
roots towards th* origin aid then ?ut el.ng the Imaginary axi^. 


J (k(p> 
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From this picture, a "transition phase” (fig. 1.2.21 %up defined in 

which our attention was focused cn th# motion of the interior roots. 

We found that their motion was well-defined accept for t> very email 

sot of possible motion* near the sticking motion. 

Th* same methods were employed on II, with th# fol moving 

resultst 1) Wc found that its transition ?hae« was appro: Ana t#: / t.he 

same a* that defined for I to fir#-; order in tho smili j iramoters 

b~* ™ and fc“ 3/ ^ 2 37 . 2) A critical angle >. was der:v*d which 

dt dt 

defined the sticking motion in which th* pendulum went <rv«r the top 
for the last time, reversed .<ign, and again #lcwly Appro* rhed th* 
top, tending to "stick" there. This particular motion itpar#t«d 
transitions which led to libration of the pendulum from transitions 
^ ^.ich led to reverse rotation. 3) Also, a probability for cap' ire va* 
derived, g'ven arbitrary kinetic energy of the p#ndulu«r a* it went 
over the top foi the last time. 

The other major development involved the orbit-orbit interaction 
and how it related to these dimple pendulum examples. A one- 
dinension-a Hawiltonian was derived that is a vr id aj* rexime tior , to 

ly. 

0(u-^j, of t,.s very long perLod motion, including the i*<.ular 
** o 

motion due to the tidal influence. Tils third Hamiluot lap (rii) has 
the form (see section 2,10 for definition cf terms)* 

III* H(x, + ,t*. • l/2> ♦ c(x,tH" ♦ b(x/cc*$, (3.1.1a) 

where it is understood that th# partial derivative wit', respect tc * 
doe* not operate on th* x dependence In c(*:,t) (section 2.9 . TJut: 
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equations of notion trot 


<U 

dt 


3h 

34 


-bl*)slr.4 , 


0.1,1b) 


<*'*> - h txJcoe*. c) 

dt * x * 


Moreover, tho valiHty of lit extends to both tha Ubratlonal and 
rotational phases »>f t> * resonance variable 4, so that it can fca used 
to discuis transition. Khan restricted to small Inclinations and 
eccentricities, the fractional fluctuation of the coefficient of 
the potential tarm(~) As determined to CM^) by its leading factors of 
e and I. In addition, the simple e and X type resonances f*>r which 
k equals 1 had identical functional forma for b(x>. This suggests 
that tost of the physics of transition can be discovared from a 
thorough investigation of a few veli-chosen examples. 

The |k| • 1,2 case* of the simple e type 
(i.«. o ■ j) ♦ J')' ♦ h*»> wi'.l be these examples. The reasons arei 
1) C - |k| - 1 can i comprise:; most of the observed resonances; 2) 
the | h ; • 2 case 1 * equivalent to a first order Taylor extension 
of b<x and la the; afore a practical generalisation of all possible 
tvo-bcdy resonance!- 'except synodic) In the snail fluctuation limit 


(«x *<• L ,r ,Z >. 
o o o 


The principal topic of this chapter vlll.be to understand tha 


transition mechanitss operating in the orbit-orbit Interaction in 
terms of the pictuie develop'd in section 1.3. Although such 
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affects as the secular change In tho roots />• tho syster approaches 
the transition phase are relevant, ve shall defer isir «xpUc..t 
determination until section 3.4, along vlt> other topics related to 
the recular cMp;t in the elements in tach phase, Including the 
amplitude of Jibraticn. Tha principal analytical tools * .11 be 
1) action integral and 2) averaged equation* of notion f <> * tits 
turning point:; of x. 

There are tevaral steps we can perform preparatory to an 
explicit examination of transition, among these being the derivation 
of the equations of motion, and a specific labeling set ere to 
identify roots. To further simplify the problem, we shall adopt ~ne 
following conventions, which in no way restrict the geuor alr^y of 
.he problem. 

a) The "noreal" liberation canter (obtained from blx) - b(3)> is 
at mod(2e), or equivalently, h ■« 0. if, for a <;..van 
resonance variable 4- the nomal center is a nx»d>>, then a 
new resonance variable can be defined (1 - 4 + s , for which 
the normal center is at itodf?*?). 

b) The angle 4 i* constructed such that U;e ,jdal torque 
secularly decreases i . This seers that 4 ie putative 
before the conwnsur ability is established. It 'oilovs t)^t, 

in the positive rotation phase far fren transit . ?n, Ue 

— <1#~ 

function cUft) is less than taro ard — Is greiter i-han 


zero. 
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c) In th« positive rotation phase far fevn transition, choose 

r o such that the tins average of x vanishes. Recall that the 

action J p OJ<rot< else vaniahes given this particular choice. 

Another consequerce is that H . * l/2c* , 

pos • rot • 

Vrcw the ir-egral solution for x, we find that the ro t* of the 
quartic polynoaJ R(x) bound the notion of x. This statement is 
trve even if these roots are time (or momentum!) dependent. Hie 
roots are obtained from f.(x) • 0, or 

RU) • b 2 (x) - (H - 1/2U * o) 2 ) 2 - 0 (3.1.2) 

For the simple e-type resonance, the function b(x) has the fora 
(2.7.8a, 2.4.4)! 

b<X> - J- A o ^ C0(x (3.1.3) 

• '"I*) ~ (2i,’ l {-to« ■ r ; 

o o 

where 

-t * l/3e :, L >* 0, (3.1.3) 

o r o 

?: » J rn»l tr * >* nanr.pulatlon, the veriablet x, c, H, 

b{x) au i t wilt be re's. m<<3 sui h that the roots of (Rtx) ■ 0) are 
dimentionless, '»hilt thr sedations of motion remain canonical, 
rirst, divide R(x) by and retefine x, c(x,t), p (which occurs in 
c (x, t} ) r and K by 
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*»(»:»- « 0 t-fcx -♦ 0 . 1 . 9 ) 

7hr four roots, ol tn lned fro® R(x) - C, shall be labeled 
*« cording- to the folloviirj scheme. 

e) First, if a gtnen root satisfies the aquation 

f l/2(* r ♦ O 2 *• j(x r i# U. 1.10a) 


h root x r is an implicit function of the variables x ml 4 thiough 
the functions H and c. The situation of motion can ba acp.^ratisod asi 




ci.i.m 


Thi< partial derivatives of x with re* pact tv c end H cai as 
expressed as functions of tha roots using 0*1. 10) i 


then that roct ii a 2*-roet Those roots which 

satisfy 

n * i/2(* r ♦ «:> 5 * »:x r > b) 


shall ba labels i-roots (x^) , 

: } The following rulas shall ba used to distinguish between 
w-roatsi 

1) If there exists only one i-root, tha label is unique. 

2} F >r the cere where two n- roots are reel, label them 
x and x with x > x m . jf these two roots are 
complex -'I jugate, again label then x^ and x fl- , with 

*,♦ * » x ,.- 

3) There can txist three »roots for the )k| •• 1 case. 

Label zh teas x^, x f+ , xj^ if x^ < Re x^, or * w4 , x w- , 
x*_ if x^ > 9m X vm . For three real roots, label then 

*.-• x ,.' «:♦ if *.-• **,+ **u • 

Similar rules can be used to uniquely label the 2*-roots. 

The next step is to derive the equations of motion of ths roots. 


»\ 

w 





l 


alee, the fell ti*s derivative of H ngusls its partial time 
derivative. Carryinq mt these atepn, ths equations ot motion of 
t w root* *ft 


dx 




Wh.r. 


♦<V • - (x , ♦ 0 * W>- 


♦‘V ■ "* x 2« * c ♦ b * u 2. ,> - 


(3.1 .12a) 

b) 

(3.1.13a) 

b) 


For the esse Uhere two * -roots ere couple* conjugate, their 
equations of motion can be separated intc real and Una< inary parte. 
Tha rssult is 
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x^ • Re ♦ 1 XA 


U. 1,14ft) 


IT ’ - dt^ ( <* * *• Val?" *, 


The partial derivative* of Bn x and Re x w sen be o talned from 
(3.1.9b) after -separating these equations into their real and 
imaginary parts. This procedure works veil for the |k| » 2 ease, 
but not for the jk| • i case. A better wet Hod for this latter case 
is to separate the quart ic polynomial (R<F.e x ♦ 1 Bn x ) * 0} into 

V If 

its real and imaginary parts, and from then# relationships determine 
the partial derivatives as functions of the real and imaginary 
parts of X*. 

Pe'ore proceeding, ve shall consider thj effect of 
approximating clx.t) by c(0,t) during the transition phase. To 
estimate the effect of tha action variable in c(x,t) on the 
probability of capture into libration, let's approximate b(x) by 
b(0) l t , so ttatt the pendulum motion ie approximately that of XX. 
The des<rlption of the transition phase will be qualitatively 
similar to that fer X, to lowest order in the small parameter 
e" 1 (fig* 1.2.2). The equation of motion of the imaginary part 
of x^, during transition, is 
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expressing «nd (x e c) in te.*3os of ♦ , using 12 . 1 36), vc 

find 

2 

dim x dc « * 2 

1/2 • *>♦ • (3.1.16) 

l 

At transition, the term pc(t) is of C (g| V 7 ) and w:..l be neglected. 
The transition phase integral determines the e.ntic-1 ancle * ic for 
which the sticking motion (i ♦ 0-) o* curs. T'* conditior is that 
Jm 2 x(f) * I» 2 x(i) * 0, or, after int-tgr .cing (2,1.16) < 

0- . H«* f - * li} ) - »L[ <3.1.17) 

To lowest order, ^ « and H i -8 during transition. Therefore 4 

is 

4 2 * -26U ♦ eosi), 
s*al (3.1.17) reduces to 

tf 

0 " |t <1r * *10* * 2p6 ^ti 4<al * 0Q9 * U (3.1.10) 

Since 4 lc lies in the range « ’•* * 3s, and 6 • C, ve can 

deduce that capture occurs only if j> < 0. The above aquation 
(3.1.10) is identical to the relation which defines tie sticking 
motion for XX, if the factor (j|) if replaced by (2pB). therefore 
the probability for capture can be immediately derive* sy replacing 
the same factor in (1.2.3a). I’he result is 
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It is unfer stood that if * > 0, then * e - 0, And If > 1, capture 
occurs (ct all po#»lbl© ijJtial conditions, or has unit probability. 


- 1 / 2-1 dc 
dt 


is, fcr the sinpl»-e case {3,1.3, 2.10.6), 


-1/2-1 do 

P It 


p « i A ‘ 
o > o oxx 


'1ST ♦ i’s? 1 ' 


16(j2 3^ * 3 ' 2 r1L ^3; e 

dt e u' L o at o 


'J.1.20) 

||* ||(.| 

Th * n^nltjd of tb. «S)uv. t«m la of o{< , which ia, by 

*0 

assure ;lon, vary larqe. But this Implies that P , dus to ths x 

c 

dependence in c(.M,t) is very small. 

Inci^sntly, cf ve chcose c (x,t) to equal 


: (x#t> * - t o ) ♦ p/tx + c <x,t) >c3t, 


( 3 . 1 . 21 ) 


than ths second order equation for 6 is* 


♦ - bsi«4 - - o4> 


The left hand side of this equation la formally Identical to the 
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equation whioh describes the spln«orblt interaction, while the rlqht 
hand aide of this equation ie one possible font for the onaentasr- 
dependent torque actinq on the affected planet. I’ainq tala equation, 
Goldrelch and Peale (1966) derived a probability for cap taro identical 


to (3.1.19). 
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3.2 TRANSITION THEORY FC« SIMPLE ECCE N TRICITY- 
DEPENDENT RTfOWAHCE* U * jk ♦ j'V + k*) WHERE |k| • 2 

Because of its mathematical rteplicity, tha |fc| - 2 cam will be 
exaninod first. Tha explicit solution of tha roots as functions of 
c and H Is partlcul arly simple, sinca fc(x) is llnaar in x. The 
sQKplicit solutions for each root, along with their idsntif ication, 
ore i 

*»± * * (c * kS) 4 kB)* ♦ 2{H ♦ 6 - 1/2 q 2 ) 

• -<c ♦ k8) ivXi ' (3,2.1a) 



Vl , dc " *2*t } 

dt dt */A~- 4BCkcT l) 


Since H is of 0(l/2c 2 ) in tha positive rotation phase, far from 
transition, both the roots x and x are of 0 (-) while x_ ^ and 
M are of 0(-2c). In the positive rotation phase, x is bounded 


IflU . MS ***** « 


143 

between the mailer pair of roots jc^ and x r- . insp-iction of the 
aquation of motion for x reveals that the moat negati ’*• value of x 
occurs at $ • mod (2v). Thus we can deduce x. < x The relative 

0 V— Tt — 

location of x 2 ^ can be discovered by ir spec t ion of their 

difference i 

x 2w+ - x fl+ - (3.2.3) 

for k * -2, X 2v+ > in th * Positive rotation phase* If k • *2, 
then > x 2w ^ far from transition J v >> Ib^j). Put near 
transition, these roots may interchange their relat.va positions. 

We should point out that in the positive rotation plase, the value 
cf x for which b(x) vanishes (or a •• 0) equals k' 1 aid must be lass 
than x- when k < 0 and greatar thAn x. wh**n k * O independent of 

<6"~ r- 

the magnit Ae of k From an inspection of (1.1.9), *f« can deduce 
that for |k| « 2, a r and a 2* root are equa- if and only if they 
equal k~ l . Conversely, if a * (or 2*r) root equals V * , then there 

must exist at least ona 2s (or s) root which equals k~*. Therefore, 

• J, 

ws can conclude that 1) for k - -2, no root equals V in either tne 
positive rotation phase or the transition pUise; 2) if K - 2, then 
the roots x , and x. , interchange where b(x) vanishes. 

Ht 

This information is sufflclant to construct the pouslbls 
rotation phase diagrams for each case (fig. 3.2.1a~c). Incidently, 
t)>e principal qualitative difference in the positivn rotation phasa 
diagrams for this example and XI is that in 71 (f ig - ■ 1 .2. la) the 



Im x 
* 


- 1/2 

(e=0) 


| J 

-(c~£) 
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*2^ 0 


Mfce *W 


*2v+ 


a) Positive rotation phase for k • -2. M and M, are the aid* ,ints 
of the * and 2* pairs of roots, respectively. * 


Im x 




c) Positive rotation phase for k - +2. In this diagram the roots 
* 2ir+ and are to the right of x • +1/2, where b(x) vanishes. 


FIGURE 3 . 2.1 

POSITIVE ROTATION PHASE DIAGRAMS FOR |k| » 2. 
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midpoints of the pair of * and 2s roots co'ncide, while in this 
example their respective midpoints, M fl and are offset by an 
amount (-2k$) . There exists only one possible diagram for k « 2 
(fig. 3.2.1a), while there are two possibilities for k • +2 
(fig. 3.2.1b,c). Figures 3.2.1a,b are qualitatively similar to the 
diagram for the simple pendulum (fig. 1.2.1a). in addition, the 
transition phase for these diagrams will involve the f*roots which 
coincide, and then develop an imaginary component. 

Figure 3.2.1c suggests thac there exists a radically different 
form of transition directly into the libration phase, in which the 
roots x and x* exchange their relative positions on the real 

*T“ 2 *+ 

axis. The equations of motion of these roots (3.2.2) show that 
this automatic transition is an allowed motion of the roots. 

Whether it occur e depends explicitly on the value of the parameters 
et H, and 6 at transition. Subtracting x end «... at 

n m 2*1 + 

coincidence we findt 

x_ - x w - +2kP+A“'-~4B(kc + 1) +/ST - 0 . (3.2.4) 

Zw+ 

If k - -2 the above cannot vanish since $ is negative. Transposing 
(2k0) In the above elation and squaring, we fiivi that H must 
satisfy the relation 

H + 1/2/Z /A - 48 (kc ♦ T) • 20 2 C3.2.S) 

at coincidence. 

The condition that separates these two typrs of transition is 


£ •* ,A 
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for the two v-roota to coincide at the sane moment that x^ ar.d x^ 

2 

coincide. This requirement is satisfied if A - 0, H - 26 and 
c • k3 • -1/2. It turns out that this information is not sufficient 
to unique!/ sepcify the parameters. Another relation can be obtained 
by evaluating the action in the positive rotation phase for the 
special case where x^_ • x^ - x^ ■ 1/2 (see B. 19,20). The result 
1! 

8 - -1/8) c » -1/4; H - 1/32) - -1/2. {3.2.6) 

If !e| * 1/8 and k » +2, the system automatically makes a transition 
from the positive rotation phass into the libration phase. For 
| Sj < 1/8, we expect that capture into libration depends on the 
initial conditions, much as it did for H. Transition of this type 
begins when the w-roots coincide and move into the complex plane. 

The w-roots are complex if A < 0. From inspection of {3.2.1a), 
the real and imaginary components of x^ are easily identified and 
are 

Be ■ -(c ♦ kS) (3.2.7a) 

Im x w+ - ±»^A . *>) 
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tr-roots first coincids and the time t £ that x - Re x ff . The result 

is 

l/2«» 2 x n <f) - a/x/il) • |f *„-*>, > o 

( 1.2.9) 

2 

By construction, Im x^U) * 0. The condition which separates 

transition into libration from transition into negative rotation is 
2 

that Itn x^(f) equals zero to first order in the mall parameter 

—1 dc * 

6 ~ . Recall that this dsfines ths sticking notion where ♦ »> 0-. 

If Ira^x^ff) > 0, then the imaginary coemonents are nonzero as 

x »> and 4 reverses sign again. The conclusion is that the 

system has successfully entered the libration phase. But if the 

integral is negative, the implication is that the system has rade 

the transition into the negative rotation phase. 

The integration variable in (3.2.11) can be changed from t to 

x. Then the condition that the integral vanirhes will determine the 

initial value of x, x ic , which leads to the stick 4 no motic. t »> 0-. 

Unfortunately, it turns out that explicitly depends on the value 

of the roots at transition, which in turn depends, in a complicated 

way, on the parameter $• The dependence on 6 must be determined 

using the action J* Instead of pursuing this cou.se, we shall 

determine an approximate condition, accurate in the small fluctuation 

limit (-j << 1) and delay finding a eore accurate relationship until 


ths next section. 
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An instructive transformation is to observe that x - Re X is 

s 

related to $ by 

(x - Re x w ) - -* ♦ kfld + cos*). (3.2.10) 

The transition criterion is then 


0 - ffl* lc - * k Bf t * dta ♦ coat)]. 


0 . 2 . a> 


To lowest order , t(f) as and $ ic lies in the range s * 3s. 

Thus for k - -2, the system completely evades capture into the 

libra tion phase. Observe that this relation is similar to (1.2.26) 

if the parameter ~ in (1.2.26) is replaced by ke||. If the 

fluctuation in b(x) is small, then the motion is nearly that of a 

simple pendulum. Given this approximation, the probability P Q for 

capture can be immediately obtained by replacing 57 * by ^ 857 * in 

Qt o . 

(1.2.40). Furthermore, we can find a general result applicable in 
the to all fluctuation limit if we replace -k$ by & . The 
probability for capture for a system 1 th an action-dependent 
potential isi 




c 1 ♦ ./2,|a|" l/2 b j' 1 


4b 

1 . 


(3.2.22) 




S 


V 

*1 


The restrictions are that p c * 0 if < 0 and ? c » 1 if 

|B!* 1 / 2 b x * 2 /v. Por the e 2 type resonance, the condition that 

P * 1 (i.e. automatic transition) exactly corresponds to the 
c 




requirement* that x fl _ - x ff+ - ~2v+ ftt transition ' ot 8 " -1/8. The 
above formula predicts that P c - 1 when 8 ■ -* 2 /16. Therefore, the 

applicability of the above formula must be restricted to values of 

1 1 ?b 

) 6 1 much smaller than 1/8 in order to keep fix and hence ■ y smal- . 

For the |k| - 2 case, the "small fluctuation limit" is defined by 

the condition that | 8 * « 1/8. Conversely, the "large fluctuation 

limit" is defined by the condition that \&\ > 1/8. At the end of 

the following section a more accurate probability estimate is 

derived (see Pig. 3.3.9,10). 

There is a more serious question concerning the validity of 
the Hamiltonian in the limit where b(x) vanishes. This problem ir 
connected with the existence of terms proportional to e in the 
interaction which can cause very liirge fluctuations in u as n 
approaches zero, and is discussed more fully at the end of section 
3.3. The indication from that analysis is that the automatic 
transition mechanism for the k - 2 case is seriously affectec by 
such terms. 



150 

3.3 TRANSITION THEORY TOR SIMPLE ECCENTRICITY- 
DEPENDENT RESONANCE j U • jl ♦ j’A* ♦ kfl) WHERE |k| - 1 

The case in vhich b(x) is proportional to the first power of e 
is particularly interesting because of both its exotic behavior and 
the fact that the majority of the naturally occurring resonances are 
of this type. It includes the Bnceladas-Dione (1^ - 2^ + 0 D1 ) and 
the Titan- Hyper ion (-41^ + 3X Ti ♦ exas^les. The bag of tricks 
used on this case shall be more qualitative than the rather 
straight forward method applied on the simple pendulum and |k| • 2 
cases. The principal difficulty is that although the quartic 
equation fR(c> - 0) can be solved for its roots, the solutions are 
too complicated to derive tram them the position of the roots in 
the complex plane as functions of the parameters c, $ and H (see 
App- C). Fortunately, there are other ways to answer the important 
questions. Our past experience with the k • 2 case should suggest 
that this example should behave similarly in the small fluctuation 
1'jait. That is, capture into libration may occur , depending on the 
initial conditions, if k is positive, while it is evaded if k is 
negative. The first really interesting question is, how does this 
system behave in the large fluctuation limit? From analogous 
behavior for k » 2, we expect to find that for k - «*1 the system 
automatically enters the libration phase if the parameter 6 is 
large enough. A second question of concern is, what condition 
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separates the two types of behavior? For (k{ - 2, the condition was 
that three roots equal 1/2 (3.2.6). The final question is, how does 
the system enter the inverted libration phase - if it can? As 
with previous examples, the starting point of this discussion shall 
be an investigation of the motion of the system in the positive 
rotation phase. 

Far from transition, the motion of the real variable x in the 

complex plane must be bounded on the left by a 2-a-root (x^^) on 

the rig/jt by a e-root while the magnitude of these roots is 

of O(-). The relative magnitude and position of the remaining pair 
c 

of roots can be discovered from an inspection of the equation 
(R(x r ) - 0). For |k| »i, we find 

R(x r ) - D« * 1/2 (x f ♦ c)V- 2 <**x r + 1» - 0. (3.3.1) 

Since H ie of 0(1/2 c 2 ) , the first tern In (3.3.1) is snail erly if 
is approximately equal to either 0 or -2c. The first case 
corresponds to the relatively small bounding roots, * 2ir _ arid 
If c is small compared to 1, then the pair of roots of 0(-2c) is 
also real, since the second term in (3.3.1) is then negative. But 
if <-2kc) > 1, then both terms in (J-3.1) are positive. The 
conclusion must be that if (-2kc) » 1, then the large roots are 
complex. Thus, in the positive rotation phase, where -c » 1, the 
large roots are real if k • -1 and complex if k » +-, while in the 
negative rotation phase, the converse is true. 

If the roots are ceai, then the’r labeling in the positive 


'A 
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rotation phase nust be similar to that of the simple pendulum. That 
is, the large roots are x^ + and x 2>+ , and the two w-roota are 
interior to the 2ir-roots, the same as in 31. This description always 
applies to the fc » -1 case and would appear to lead to uninteresting 
behavior during transition in that it always enters the negative 
rotation phase. Although this is true, as shall presently be 
demonstrated, this peculiar pendulum can, for a time, librate in an 
inverted fashion. 

For k - +1, these large roots are complex, if -c >> 1. They 
can identified by solving for the roots in the positive 
rotation phase for typical values of c, H and 6, and then substitut- 
ing them into (3.1.9), which defined the 2s and e-roots, respectively. 
The result of this exercise is that these complex rents are both 
s-roots and, i*. the positive rotation phase, can be labeled 
and Another piece of information that ocld be useful later 

would be to have some idea of the magnitude of the imaginary parts 
of these "exterior” t -roots. We e;<pect the imaginary part to be 
small compared to the real part (this shall be rigorously 
demonstrated later). If this is so, then we can substitute 
x - Rex**- llmx directly in the relation defining the n-roots 

i tr it 

(C.1.9), expanding the factor (-kx^ + l) 1 ^ 2 and separating out the 
lowest -order imaginary component. We find Im is, approximately, 

2c tm x t * 1/2 8 [Re \\~ 3/2 , (3,3.2) 
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which is small if |Re x ff | >> • 

For the case where all roots are real prior to entering the 

transition phase, we expect that transition is initiited by the 

coincidence c i the interior w-roots, If the roots x and x* 

t+ 

remain complex through transition, then something analogous to the 
automatic transition from rotation into libration found for the 
k - +2 case must occur. Pecall that this phenomenon happened when 
in the large fluctuation limit both e(t) and 2 were relatively 
large parameters. 

To discover the equivalent phenomenon for the |k| - 1 case, we 
shall assume transition occurB vhen |c[ >> 1. The maximum fluctuation 
of x must be of 0(1), since b(x;c* <-kx ♦ l) 1 ^ 2 * 0. This implies 
that the Hamiltonian can be approximated by neglecting the x 2 term 
for values of |c| >> 1, since the real part of the complex roots is 
of 0 (-2c) . Physically, this means that the variation in the mean 
longitude of either resonance partner is small compared to the 
varition in the per ic enter of the unprimed partner. The result is 

H - V2c 2 * h - xc + 8(-kx + l) 1/2 cos*. (3.3.3) 


If one of the turning points in the motion of x occurs at x - k, 
then h equals kc. This condition defines the value of h for which 
b(x) periodically vanishes whenever x equals k. The values of h for 
which the system either librates or rotates can be determined by 
substituting ack for h in (r.3.3) and solving for cos*. The result 
is 


w ,.,. 



154 


cos* 


kc(-kx + g) 
et-tor ♦ l) 1 


13.3.4) 


vh«r« o ij an arbitrary dimensionless parameter. Tbs physical 
solution must be such that Jot < 1 for all allowed values of x. if 
the pa ram er a is < 1, then It appears that the right hand side of 
(3.3.3) can t or. both positive and negative values, with the 
Implication that the system has either rotations or librations of 
amplitude > 90V (later we shall discover that only the former is 
allowed.; But i>' o > 1, then the left hand side cannnt change sign 
as x varies. Si «ce |* > 0, the turning points in the motion of x 
must occur at cose - h. The irplicaUon is that this system librat 
about the 2* -center for Jt - +1 and about the * -center for k - -1 
with amplitude * 93*. 

To better understand the behavior of this unusual pendulum, 
let's reconsider the integral solution of xi 


/ sign(-Bsin$) » t - t . (3.3.5a) 

/k[x) 0 

The function R(x) now equals 

R(x) • B 2 (-hx ♦ 1) - (h - cx) 2 , b) 

and la a quadratic polynomial in - instead of a guar tic polynomial. 
That is, there are only two bounding roots, and x_, giver, by 



- ~\^ch - kB‘) ±|e|/? ♦ 4c(c - kh)3. 

2c 


( 3 . 3 . 6 ) 
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In the positive rotation phase, the roots x + , x~ can be equivalently 
labeled and x 2ir _, r,l| P* ctively. Since the labeling is unique 
until h ■ ka, we must conclude that the system must rotate up to the 
critical value of h - ka. 

The 4 ta-iral solution for x, assuming that the parameter c is 
conet*' ?n particularly simole in that the integral is 

pi >po c arcsine of a linear function of x. The solvticn 

tn .* .xon of t 1st 


< x - x )sln|c| (t - t ~ 


< 3 . 3 .*>) 


Observe hat the frequency of the nation equals |c| and becomes 
small as |c| «♦ 0. But transition is already presumed to occur 

when Jc | » 1 . Therefore, the secular behavior of the system as a 
unction of c c«n be obtained from ti ction integral J, since J is 
sdiabatically conserved when the oscillation frequency is large 
compared to the changes in the system's parameters. 

In the positive notation phase, the appropriate integral 1st 


2 * 

J . - L xd$. (3.: 

pos.rot. J 0 

The first step in evaluating this Integral is to express x as a 
cunction of $ using (3.3.3 ) t 


x 
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* *?.-:• th# turning polrts In she motion of x occur at )oo*+) • l, the 
t«o jolutlo't t» - implicit in 0.3.6) mutt be nebeh*A at either 
* * - ->t v..w. .he argument of th« radical in 0.3.9) vanishes, 

i ha first choice appllea If th* pendulum routes and the second if it 
.'ifarate'. Tn the positive rotation phas«t, the mtching of the two 
solutions lr*ds to the folltvlrg result i 

*pos.rot. * **^l 2ch " ►® 3 co» 2 6 ♦ 6coe6*4^o'e% ♦ 4e(c - kh)]. 

(3.3.10) 

Tn the 11 brat ion phase, t*n argument of the radioal vanishes wheie 
o reverses sign, Also, * J1J| is a minimum at 6 - 2», 4 > ,0 for the 
y • 1 caie, while it is a mexlsnss at 4 «• a, ♦ > 0 for k ■ -1 . Thi^ 
mplis# that the contribution from the aecond term in (3.3.9) must 
be multiplied by sign Br/licitly» 

x Hb. ’ - kB 2 cos 2 6 ♦ signt-ki) | Boost /& 3 cos 2 t ♦ 49(0 - kh)]. 

-C 

(3.3.11) 

Ofe.arvc that th# flrit t.rm In (3.3.10) f or * . la an avan 

pos . ro t * 

function whi* the second ie an odd .unction in 6 over the interval 
{0, 2w). "tils means that the integral of the seosnd tern, over the 
range {0 < t * 2*) must vanish. Thus, ^ 9 lven by 



-—r(4ch - \0 2 ) . 
2c 


U.3,12) 






isa 


amplitude of Iteration 4 m at a function of the parameter c. The 

extreme In tl.a notion of * *** obtained from the 

condition that the argument* of the radical in <3. *.9) vanish. This 

c^~u3ition reduces to th«. fallowing relation for 4 » 

in 


-sin* e 
m 


2c 

T ' 


<3.3.17) 


or eln* m decreases linearly with e(t). 

The folioring der icipUon of tranaition can )#e deduced from 

these facte. Tor k - ♦!, the root approaches - 1 ae c approaches 

the value, l/2d. This n.;ens that the function b(x> which it 

proportional to <-x ♦ l) 1 ^ 2 x» smaller in magnitude at * - modfir) 

than at * a mod (2*). At transition, when c e 1/26. the function b (x) 

vanishes at * • mod(2v), whl ) for |e | < 1/2 |s |, the pendulum 

libra :ea with maximum amplitude * • 90* about the 2* center. 

m 

Figaro 3.3.1 shows this ssquonc# of events for the k • *1 ease, 
while figure 3.3.2 shows the equivalent sequence for the k • -l case, 
whi.h for |c| < Iterates about the « center. 

By the way, it is no accident in figure. 3.3.1,? net the path 
traced out by the pendulum in each phase is e circle o' constant 
radius. Greenberg (1972b) observed this phenomenon in his analysis 
of a similar problem. Furthermore, he found the angular velocity, 
as measured from the center of that c.rcle, is constant and equal to 
|c|. To prove these assertions let's consider figure 3.3.3 which is 
a diagram of the path traced out by the pendulum in ti e positive 
rotation rh<is«. The center of the flours is at D *here d « Cf> is given by 




a) Positive rotation phese. Here |c| > 3 |b| end (-c) is e 
monotonicalJy decreasing function of tlm*. 



b) Here the root equals ti, while the function |b<x)| f vanishes 
et c - mod (it). w In addition, the root * 9m changes from n to i 
2* root* 



,c) Li brat ion phase f the pendulum Iterates about * w mod<2*> with 
amplitude 4* * 90*. 


FIGURE 3.3.1 

TRAN* 1TX0M IN LIMIT ;fl| » 1 FOR k • el. 
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a) Positive rotation phase. 



b) The root x equals -1 and changes its label from a 2s to a * 
root w ~ 



c) The pendulum llhrate* about ♦ • • vl th sncl Lttd i 4 90*. 

ID 

FlGUKt 3 , 2,2 

TRANSITION IN T.IXIT |fi| » 1 TOR fc • -1, 


4 i >»- *i 
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d • ^(Wc« + ♦ ]) V1 - l-kx^ ♦ t) va l. u i.lt) 

Thus, 0 is the sldpolnt of tha m ax l as j» and minimum values of L>U)cos#. 
Substituting h ti given by 0.3.12) in the equations defining * 4 and 
x_, we find 



t 2 • d 2 ♦ 0 2 >k* ♦ i> ♦ 2d(M-kx ♦ U l/a cos6. CJ.2..1) 

Eliminating thi oos6 dependence using 0.3.2}, we find that the 
explicit x dependence occurring in nich tens cancels such that 
r 2 ■ 0 2 ■ constant. That Js, th# fijurs is a circle of constant 
radius. It is no ,>r'-blsm to show thit this result »lso holds in 
the llbrsticn phase. This means that the absolute variation in btx) 
remains constant although its mean value grows as ]c(t)| decreases. 
Inspection cf ':ha figures denting the motion of the pendulus 
through transition reveals t.*afc the relative derping of the 
fluctuation of bis) Is proportions* to c(t)» vanish Irg as c(t.) * 0-. 



The fact that T La a , linear function of the tin* follows directly 



FIGURE 3.3.3 

Path traitd out by perv3ul\dn governed by (1,3.3) in 
positive rotation phase f.r k • H cats. the pa rasa* ter 
d is the midpoint between the minimum and maximum values 
of b(x). v is the angle made by vector r with z aspect to 
tie origin at D while $ is the angle node by vector 
btx) with respect to the origin at 0. 


free the law of ilnee. Explicitly i 


slnT ->Jn4 
*3 (-kx ♦ 1) " ~ 

Using the equatisn of motion for * (3.1.1b), find (3.3.6), 


0,3.22) 
we fine! 


IdlsinY •• • f l30 slc|(t - t ft >, 


(3.3.23a) 


or 


f • | c ) ( t - t Q ) f mod (s ) » b) 

Eventually |c(t)| decreases enough to that the Kairlltonian 
h(>:,*,t) (3.3.3) le ro longer a valid approx JMt Ion. s expect that 
the Mpproximai.it n breaks down w)ten the x 2 tarn In H la of ths sane 
metfr .tuds .*t thn cx term, or until x - |c|. From (3.3.19), the near* 
vai.ua of x is 


(*> • l/2<«. ♦ » ) • . 0.1.24) 

4o 

Therefore, the lelnlmia value ef |c(t)| for which the Haalltonlar h 
is a useful app:oxlm>tion of the motion is 

III’ 1 '* • (2.3.25) 

71.1s means that 1 0 1 lust be a very large master if there Is to be 
significant dasping of the icnplitudi >f ' ibreUon via this 
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owe ha r. is*. 

So far wt have a qualitative picture of how the pendulum 
behaves in both the small and large fluctuation limits fo* the 
M • 1 case. The next stage of the development will be U> determine 
the critical values of the parameters which separate the tvo types 
of behavior and to develop a wore accurate probability argument. 

First we shall discuss the consists qualitative behavior for the 
k • -1 case, since it can be easily deduced from the previous 
arguments. 

Far from transition* the relative position of the four roots 
will be as shown in figure 3.3.4a. Ai» the system evolves* two 
possibilities for the transition phase exist. If the left bounding 
root does not reach the value -1 prior to the coincidence of 

the -oote, then transition uill involve the temporary motion of 
the *-roots off the real axis. As with the k ■ +2 case* the 
variable x will move pest the point where x - Re x^. Thereafter 
the components iln x ff will move back toward the real axis* reaching 
it before x returns to the value Am x^. This means that the system 
haL entered the negative rotation phase without any possibility of 
librating (see fig. 3.34b*c). The other possibility is that the 
root * 2 *- < *°** r**ch the value -i prior to the coincidence of the 
interior w-roots. The pendulum then begins librating about the 
la* rted position, eventually the r -roots will coincide and become 
complex. The motion of i will vanish and reverse sign somewhere 
near the top ^-position (fig. 3.3.5* ,b) . Since x is increasing an 
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. a) Diagram of poei~lv*i rotation 
phase in complex x-flane 



D.\agr«m of path traced out in 
polar coordinates (jM*)!,*) In 
positive rotation phase. 



( b) Diagram of the mo .Jon of x -Motion during period when i 

just after x-ri»ot » coincide. reverses sign. 



. c) Negative rotation phase. 


FICUFJt 3.3.4 

-mVEITTOW DIRECTLY INTO NEGATIVE ROTATION 
PHASE FOR k - -1 AND |fi| « 13. 





t ■- 


P1CURT 3,3.5 
•mr* nf 


v 
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it novea past the position where It aqua la Re x^, the coefficient 

t/2 

of the pendulum potential b(x) which la proportional to <x «■ i) ' 
*l*o increases, after ♦ rtvtrHi sign* After x reaches th« root 
Xj. (+ / the variable x than move* bach toward the positio. Re 
Meanwhile the imaginary component cf x fl haa vanished and aga'.n the 
system antera the nogative rotation phaaa. 

The epecial situation which separate* thane two types of 
behavior is the foliowingi Let the root x 2f? _ equal -1 at the same 
instant that tie two interior ir-mota coincide. Proa the action 
integral (App. B.23i the valuea of the parameters for this epecUl 
case aret 


B - -13.89/ - -6.670/ H - 217,3/ 


x 2l| _.-J> » B - *.780, x^ - 18.12. 


If \e\ < 13.61, the pendulum behaves normally during transition. 

Por |d | > 13.89, the s" steal temporarily sntara the inverted 
libration phase. 

How than the h - -1 caae hat bean disposed of, we shall 
concentrate on the more subtle behavior of the k ■ *1 case during 
transition. With just one aore piste of information, we con outline 
the complete qualitative behavior of the pendulum during transition 
for the k • ♦! caa<i. We assert (and will prove later) that the 
function itO (3. 1.8a) vanishes wh^n the exterior n-ioots, and 
Xj + , reach the rea* axis and coincide. Recall that i(*) or^urred in 
- » f • fi-roota we l ve 




$£. 

dt 


♦ <*> 


C3.3.27a) 


♦ <»> - -* t - c - ^d(-kx v ♦ 1)~ 1/2 . b) 

Proa the above equation, the only allowed notion after coincidence 
ia for the two root® x^ and xj + to separate. One root (x^) noves 
towards the left bounding root * while the othir (x* ) noves 

w* ffT 

towards the x » *1 position on the real axis. Prom our earlier 
discussion we expect that when the root • +1, it changes from 
a r to a 2* root (x^) and thereafter decreases. We also know that 
when the interior v-roots coincide, again vanishes. This 

means that this function must be double valued in x • Prom the 

II 

equation {' {») - 0), the following relation can bo derived i 


§» — l-»» M + ♦ c). (3.3.28) 

Pigure 3.3.6 is a graph of 1 S | versus x ff+ for the special value of 
c - -1/2 and k • +1. This special value results from evaluating 
the action integral at the instant and coincide (B.23). In 
addition, the posit on of the left bounding root x fl _ is graphed as a 
function of 8, using (C.2a) and (B.16a). Observe that the position 
of x w _ i* to the left of the second -.cro of ♦(c). This ia to be 
expected since t v e motion of x^_ lo toward the position where 
(x^.c) vanishes. {Dote that this will occur ftr a more positive 
''slue cf c than c » -1/2). According to figure 3.3.6, there exists 


a m a xi mum allowed /alue of | 0 [ for whi cb th« equation (J.3.2B) is 
satisfied. This critical value for 6 and can la detarmined by 
maximizing {3.3.26) with respect to x ff . The results are 

6 , - — - -0.27221 X . - | • 0.8333. (3.3.29a> 

cl ^ ncl 6 

Incidantly, the remaining values of the parar*ten« for this special 
case are (B. 21) t 

H m V “ " l/2 ' b) 

For |fl| < th * value of x^ where the exterior e-roots reach 

the real axis is in the range |- * x ff * 1, while t svc^.kI zero In 
♦ (e) for c - -1/2 lies in the range {1/2 * x fl * . Clearly, at. the 

critical value of $ * 0^, the three e-roots coincide. Por value* 

Of | 6 | > |$ cl | :he implication is that the exturicr e-roots never 
do reach the real axis during subsequsnt evolutlor of the system. 
Therefore transition is qualitatively similar to the earlier dis- 
cussed approxinatlon. The situation is more complex for |fi| < |8 cl ). 
As the system evolves, the roots x^_ and x ff+ move towards coincidence. 
If these roots coincide before x* - U, then the system can either 
temporarily evolve into an inverted libration fhase or remain in the 
positive rotation phase as s* own in fig. 3. 3. 7-8. Which occurs 
dspends on tho motion of the imaginary components of the interior 
it roots during the transition phase. If the imaginary components 
first move of* the real axis and than return while x is to the right 
of Re x fl , then the motion of x is eventual ly trapped between 
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X 


figure : i.6 

Graph of |0 1 versus for the special value of *1,2 
and k - el, obtained from the {♦(»> - o) e, ration. From 

the solid line ( ) , one finds the value of the roots 

and at coincidence as a function of |s|. The 
dash line (— ) gives the location of the second aero of 
the (its) - 0} equation at c(t) * **1/2, while the dash- 
dct line (— -) represents the location of the x root 

IT** 

at c(t) « -1/2. During subsequent evolution, the dash 
and dash-dot curves move toward coincidence. 



a) Transition is initiated by .The pendulum moves toward the 

coincidence of interior n- top « position where $ vanishes at 

roots and continues as and reverses direction. The 

they move off real axis. angle then increase! to the 

The motion of x is first maximum j* where i again 

towards x*_, next towards vanishes and reverses direction 

Re x tfl , near where i van- and then moves toward the 

iahes {near top n position) bottom ir position. After 

then bgck through Re x passing through the r josition, 

where 4 again reverses sign. $ decreases to the value (-f*) 

This motion involvos tran- where it again reverses n 

sition from positive rota- direction. Final 1/ it roves 

tion phase to yet another toward the * position and ? 

"temporary” positive rota- vanishes for the fourth tire, 

tion phase. completing one revolution in , . 



b) Second transition, in which Here, the pendulum moves into 

the system automatically libra tion phase as blx) tends 

evolves into the libra- to vanish. The maximum 

tion phase as the * root amplitude of llbratlon is in 

equals +1, then reverses the range 90* * $ * IPO*, 

direction, becoming a 2w- w 

root. 

FIGURE 3.3.7 


TRANSITION INTO LIBRA! ' PHASE FOR 8 


IN THE RANGE {**0.2104 * 8 *-0.2722) AND k - +1. 
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and x*^. Therefore, the pendulum has made the transition from 
positive rotation into inverted libration. But If these components 
never return to the real axis, the motion of x is trapped between 
x 2w _ and x*^. Although the curve traced out by the pendulum is more 
complex, and 4 reverses sign four times during one rotation, it still 

executes positive rotations. In either case, subsequent evolution 

\n 

of the pendulum tends to increase x* (and decrease B(-x* +1) 0) 

until it reaches the value 1. If the system has made the transition 
into the inverted libration phase, then as |cj decreases the pendulum 
eventually enters the negative rotation phase (see fig. 3.3,8b). But 
if the 3y*tem has remained in the positive rotation phase, it will 
automatically enter the libration phase with amplitude of libration 

s 180*. 

Whether or not the pendulum executes these exotic motions during 
transition depends, intuitively, on the parameter 6. If its 
magni”ide is too small, then x^ equals +1, then decreases and 
changes labels prior to the coincidence of the interior w-roote. 

The condition which separates the ordinary transition from ti s t»- 
stage type Just discussed is for x* + - 1 when and x^ coincide. 

From the action Integral (B.23), we fi »l that the values of the 
parameters are* 

B * -0.2104» c - -0.4705; 

C2 

* *0.6463/ x 2w _ - - 
H - 0.1407 . 

* ** v >» * •*' *•• ' 


•0.4108/ 


Cl. 3. 30) 



a) Transition in which interior 
ir-roots first move off, then 
return to real axis, trapping 
x between x and x* . This 
is a state of inverted 
libration. 


The penoulum moves towards top 
^-position (x ) whers $ reverses 
sign, next moving through the 
bottom r-position (xj^> and 
again towards and through the 
top Tt-poaitlon. Tie amplitude 
of libration i * about the — 
center is £ 90^. 



b) Here x* » +1 and becomes -The libration amplitude 

a 2ir-r5ot. At this point increases to 90* and thereafter 

the system automatically executes negative rotations, 

enters the negative rotation 

phase. 


FIGURE 3.3.8 

TRANSITION INTO NEGATIVE ROTATION PHASE THROUGH A 
"TEMPORARY** INVERTED LIBRATION PHASE FOR S IN RANGE 
{-0.2104 £ |B| £-0.2722} AND k ■ +1 . 
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Zn sunoary, ve have found that there are three dJ.atj.net nodes of 
transit ion depending on the magnitude of the parameter 0. Zf 
|fil < |e c2 |, then transition is qualitatively similar to II in that 
it involves the coincidence of interior s-roots. For 6 in the range 
{ ) S C 2 1 S |8| * |e / . 1 D. transition is a two-staged mechanism; the 
first involves the coincidence of the interior e-roots, after which 
the system evolves either into the inverted libration phase or a 
comp] ex positive rotation phase. In the second stage the pendulum 
evolveu directly into either the negative rotation phase or the 
libration phase, respectively. During this stage b(x) tends to 
vanish {i.e. e (*) -> 0). The final mode involves direct transition 
from positive rotation into libration if |ft| * | 1 * Nov that the 
qualitative behavior for the k • +1 case has been thoroughly 
discussed, the next step is to back up some of our assertions and to 
develop a more accurate pxobabillty estimate than that developed in 
section 3.2. 

First let's obtain a more explicit form for the equations of 
motion of the real and imaginary parte of the x^-roots (3.1.9). The 
partials of the real and imigir^ry parts with respect to H can be 
determined as explicit functions of the 'parameters c, B, H and the 
v-root parts by aeparatirg the equation (J?(Re p ff +ilm p^) - 0) into 
its real **.d imaginary components. After some simple manipulations 
of these two components, the following relations are obtained i 

j 372 ' 

ft! p * ImCx ♦ c) » In X - ±/-2H + + Re p i (3.3.31a) 

n H If KG P T! 


Is*****- 
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W* should point out that the above equations are valid only if 
** *** > G. Otherwise (3*3.27) governs the notion of the real roots. 

Inrldently, the condition that la x^ vanishes iraplie* that ♦(») 
also vanishes. This can be simply demonstrated by substituting H 
evaluated at x ■ x^ into (3.3.31a) and observing that the resulting 
e^refcsion is consistent with J(n) • 0 (3.3.27b). 

It appears as the ays tea evolves towards transition (for k « +1), 
that the imaginary part of the exterior jr-roots firjt gradually 

increases (sinew they are *1/20 |Re x fl J“ 1/2 for (-c) » 1), reach a 

2 

maxiauffl and then decrease towards zero. For c >> n and n < 1# the 
notion of tm is definitely towards zero. The value of the 
parameter n when the x w+ roots coincide, after using (3.3.28) to 
elimi:.ate the 0 dependence, is 

2 <-*„ ♦ 1 ) 

t, . (3.3.34) 

Since * 0 lies in the range * x^ * 1), n Is bounded between'.! 
and aero, laving its largest value when the three v-roots coincide. 

One can also show that the parameter n is greater than or equal to 1 
when the interior v-roots coincide. This is to be expected, since 
otherwise the imaginary part of x could not move off the real axis 
while the variable x is to the left oi ?e x y+ (see 3.3.30a). During 
the normal transition phase (i.e. coincidence of interior ir-rootsl) 
the parart* er e is of 0(0 1 and is small compared to n. Therefore 
the equation governing the motion of Im during transition can be 
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approximated by 

i/* *• *-■»>* t3 - 3 - 35) 

which is similar to the equation of motion for the k « 2 case. The 

parameter o la a slowly varying quantity. Therefore the factor 

( — n ' ) is also slowly varying unless (n - 1) very s «• -ly vanishes, 
n ■ l 

Therefore, the transition integral can be defined as before, and 
evaluated to lowest order bv replacing the sJowly varying parameters 
n and Re x^ by their values at coincidence. After integrating the 
above equation we have 


l/2IIi» 2 X ii <f)- l/inAt^U)), ||(— ^y)J t *0W *, - x)dt. (3.3.36) 


fly construction, Im x^(i) ■ 0. This time, let's change the 
integration variable from t to x, keeping in mind that t is a 
monotonicully increasing function. The right hand side of (3.3.36) 
then becomes 




From (3.3.36) we 


greater or less than l$ c2 |). Therefore, the above integral reduces 
to the following t 


to eliminate the specific dependence in P on . 

c ic 

finds 


2^<— 
at n 



|dx | sign (He - x) 

Ax - x^ Y I ix 4 - X) 


(3.3.39) 


From the transition phase diagram (fig. 1.2.2), we find that 
x f - Re x^ t and x ^ lies between and R * Although this 

integral can be explicitly found in terms of arcsine functions, 
this step wii.1 be deferred until later. 

Recall that if th« above integral is positive, the system has 
made the transition into the libration phase. But if it is negative 
the converse is true. The value of x^ for which the tbove integral 
vanishes (x^) separates these two events. Recall that the 
probability measure was defined in terms of tho value of i 2 as the 
system went over the top for the last time, and was directly 
related to the function Izn^x^ evaluated between fixed limits 
(1.2.36*38). Specifically, 


♦«3» 2 ♦-Sir 

p is_ 

c - *-3 t 

* ♦-w 


(3.3.40) 


Note that of the terms is evaluated in that part el the 

transition phase where ♦ e 0. The corresponds twits in the 

var table x are x and x. . Again the appro; path Integral la 

*T 1C 

(3.3.40) over x is understood. The relation which defined x^ c can be 


2 2 • 2 • 

hut (f) - im x <♦ * 0)| +1® x (♦ <0)1 - 0. 13.3.41) 

* IT . S , 


Therefore the probability P c takes the form 


Im 2 x TT (C^0) +Im x^(C '*0)1 

P S 3n. 

, • |3tt 

0.3.4 2) 

(3.3.39) can be used to evaluate each of the terms in 
substituting the appropriate limits and the value of ♦ 
former equation. We find 

(3.3.42) by 
in the 

2 * 3* 

Imx^ (i 4 0)| f - ♦ arcsind; 

(3.3.43a) 

5 » _ 

xm (♦ * 0) | 3y “ 2 " 

b) 

where 


. *», * *2T- * *2.+ 

x 2r+ * JC 2ir- 

c) 

Th. retore P is given by 
c 


p - 2 

C 1 + ^(arcsini) 

C3.3.44) 
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This for^ila is valid only fe- t ha k • 1,2 cases with d in the range 

{0 s 6 i i). p vanishes if k - -1.-2. 

c T 

If the 2 it- roots were symmetrically placed about tne coinciding 

»-roots, then S and the arcsin of 6 would vanish, implying that P 

c 

would also vanish. But this case corresponds to b(x) - constant for 
which capture in *o libration never oc-urs. The position of the roots 
at transition is implicitly a function of the parameter B. There- 
fore, the interesting relationship is how depends on B. Again, 
the action integral can be used to uniquely determine the parameters 
of the system as a function of B (B.16-18). Figure 3.3.9 is a plot 
of P c versus <B/0 cl > for both the k - 1 and k - 2 cases. Observe 
that the small fluctation limit Tor P £ where it is proportional to 

| B 1 1/2 (3.2.12) is oi ly valid for p < 0.1 for the k - 2 case, and 

c - 

P c i 0.5 for the k • 1 case. Another interesting fact directly 

derivable from (B.18) is that P - 1/2 when 8 - 6 „ and decreases to 

c c2 

rero as 6 approaches 8^. Observe that the probability that the 

pendulum may temporarily enter the inverted libration phase (k « +1) 

is reasonably la.ge. A more revealing graph (fig. 3.3.10) is a plot 

of P versus the mean eccentricity e far from transition or, 
c e ° 

equivalently, I-®-). From 13.1.3), this ratio ie 
* •» 



(3.3.45) 
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FIGURE 3.3.9 

Graph of probability for capture intc libration. P c 

versus 6/B _ for the k » +1, +2.*ceses. 
ex 
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FIGURE 3,3.10 
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Note how dramatically P decreases as a function of e for the k • +1 
C Q 

case. P is less than 0.1 for e no more than 3e *. Furthermore, 
c o cl 

the tw-stage transition occurs only if lies in the narrow range 


(3.3.46) 


Before turning to questions concerning the secular behavior of 
the s>r*ew before and after transition, there is a serious problem 
concerning the applicability Df the theory of transition developed 
for the e-type resonances, which must be examined. Recall that in 
developing the one-dimensiona. Ramil to an which describes the 
resonance interaction, the sectnd order effect of non-resonant and 
non-secular terms on the i*' traction was of 0 1 *i * ) smaller than the 
first order terms. These ;ii»ded an infinity of terms for which 
the coefficients of their respective cosine arguments were 
proportional to the first trover of the eccentricity %. If e were 
very small, then such terms would tend to produce very large 
fluctuations in the motion of the perihelion, £. But, as we have 
seen, a major part of the motion of the resona ce variable $ is due 
to the motion of & if e Q is very small. In fact, e must very 
nearly vanish if transition involves the automatic entry into 
libration pkise. The question is. could these other terms 
proportions to e effectively inhibit this automatic transition 
mechan*.™? To get a qualitative idea of their effect, let's add 
the following term to the Hamiltonian equation* 



M(-k* ♦ l) 1/J oo# (•*><»»♦. 0.3.47) 

Here » 1* constant and la assumed large compered to the me an value 
of 4, ami X la a dimensionless parameter of 0(1), Tha Hamiltonian la 
^ow given loy > 

M * 

H<x,4,t> • 1/3 (X ♦ ,(t)) 2 ♦ bt-kx ♦ l) 1/a iC-l« • 1) * * 

Icoa (wt)3cos+. 

Pa see immediately that tha effect of thla extra tana on tha notion 
of tha roots la to add a high frequency oscillation. For |k| • 1 
rase, this rapid oscillation vould not radically change tha avaraqa 
position avan for x vary naar tha value k* 1 . Tha indication is that 
tbt coupling of tha abort period tarns to tha resonance term does 
not radically affect tha |k| • 1 case. But for th* |k| " 2 case, 
this extra tarn clearly dominates tha motion of 4 for x vary naar 
k*\ leading us to be Have that tha automat io transition machanim 
described in section 3.2 for this case is seriously inhibited by 
similar terms occurring in tha expansion of tha disturbing function, 

A more rigorous discussion of this point Is beyond tha scope of 
this thesis. We should also mention that (3,3.47) does not faithfully 
represent the complete effect of a high frequency tarn occurring in 
tha disturbing function. In addition, tha following tana 

VB(-kx * l) l/2 sin(wt)cos4? V* A 


VX'I, -v*'- /^r % • ‘ • .» Sff i, **>• 


its 

also occurs In the expansion of the disturbing function, associated 
with the frequencies (4 t w> , 




la? 

Thor* art two Units in which intonation concerning the 
adiabatic damping et ia readily obtained, and they depend on 
whether the magnitude of 0 la >* or «« than 1* The toner case ie 
oddly enough/ the simplest, and many of the important results have 
already been obtained* 

For the caae where |s| >> 1/ transition ie governed by the 
approximate Hamiltonian M*, + ,t) 0,3*3) • Recall that "transition" 
involves the secular evolution of the system front positive rotation 
directly into Ubration (♦ • 90* > without the possibility of the 

system entering the negative rotation phase. This unusual mechanism 
is distinguished by two facte. Unlike the elmple pendulum the 
instantaneous frequency associated with the pendulum motion remains 
finite and varies smoothly during transition. In addition, the 
eccentricity (and the function b(*M tend to vanish during transition. 
This evan occurs whan tha parameter c(t) • 1/30. The llbrstion 
amplitude thereafter decreases rapidly as c(t) ■* 0-. explicitly 
(3.1.17) 

,ln *m " **8^ tJ ' 4,n 

Conversely, the mean value of the action rapidly blows up as c(t) *• o*. 
rrem (3,3,7/13) we find 

4c 2 


Incidantly, the above relation is valid in both the positive 
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rotation and llbratlon phases , indicating that <x> tends to vanish 
as c + -* . Of course, this is the expected behavior. 

The eccentricity is related to the action variable by the 
following relation (3.1.3) t 

£ - <-* ♦ 1) 1/2 0.4.3) 

o 

whnre e Q is the seen eccentricity in the positive rotation phase , far 
frivi transition. Zn the llbratlon phase, the average eccentricity 
1/2 '*«x ♦ *■ inversely proportional to c(t), the exact 

relation being (3.3.18) i 



On the other hand, the fluctuation Be 5 - e . ) remains 

mmK min 

constant in the llbratlon phase, and equals 2 s q . By ths way, the 

benavior of e ftvo and 6e is reversed in the positive rotation phase. 

The average eccentricity is then equal to e . while Be ie given by 

o 


(3.4.5) 


Hus behavior directly follows from (3.3.18-19), and a careful 
inspection of the appropriate diagram describing transition (fig. 
3.3,1). 

The validity of the above results extend into the llbratlon 
phase until 
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2^1. •» |er l/3 (3.4.8) 

The parameter 0 is a function of ths mean orbital elements of the 
resonance partners evaluated in the positive rotation phase far 
from transition (3.1,5). But the orbital elments available to us 
are those that the system possesses at the pressnt time in Its 
evolution. What we would like to know is hew the above limitation 
on thle (3.4.6) approximation' s validity translates in terms of the 
parameter 6, evaluated with the presently observed mean orbital 
elements. 

We have already assn that in an e-type resonance, any changa 
in the eccentricity caused by the applied torque is of 0(*~*) larger 
than similar changes in the semimajor axis of either partner. 
Therefore, the important queatlon is how does the parameter £ scale 
as a function of left)). Recall that in deriving the dimensionless 
form for the Hamiltonian (sec. 3. 3. 4-6) which described the motion 
for an e-type resonance, H was divided by the factor 

<-r > 2 - u/?.?»o 2 

o o o 

lk! 

Since the coefficient of the pendulum term ie proportional to e 1 

lb 1 .4 

the parameter 6 tends to scale like e* * For the k • ♦! case, 
this implies that the relation between B, evaluated in the positive 
rotation phase ($ o ) , and tlte same parameter evaluated with the 
present mean value of •<B now ) 
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' B o‘- 


-£2S) 3 

I - 

S 


(3.4.8) 


Thu», by 0.4.4, 6, 8), h(x, + ,t) is a valid approximation for 


I 8 now I " 0(11 • 0.4.9) 

This la significant sines for tha Enceladus-Dlone resonance, 

* 10, which nsans that the complete tidal evolution of this 
example can be detex.,insd from the simplified Hamiltonian. 
Furthermore, it appears that its presently small amplitude of 
libra tion of 0(1*) can only be explained as the result of tidal 
damping via the mechanism just outlined. 

The last question connected with this topic iai How dees the 
dissipative term in ^ Cx,t) affect the damping in the limit 
[fj >> 1? In this case it happens that the contribution of this 
term is easily obtained by taking the time average of the equation 
of motion for h over one libration period. The equation of motion 
for h is x 


dh 

dt 


f£ <t.x)* 


The tidal, torque is (2.10.6) i 


(3.4.10) 


dc 

dt 


(c,t) 


•—•(0, t) ♦ px - (1 ♦ fcCx) + • ^*))r 

at at 


1 



< 0. 


We have added an subtracted a term proportional to W in |~(x,t) 

dimes we want to emphasize that it is the fluctuation In x and not 

its actual value which contributes to the damping of ♦ . The 

parameter 1 ie a dimensionless regative*definlte constant of O(e^), 

which means It is small compared to one. Therefore the term 

(1 ♦ l<x) ) is equal to one to 0(s 2 ). 

o 

Taking the time average of (3.4.10), we find 

; §| ,<*> . . <*> 2 )>. (3.4.12) 

Th. .la. avaragaa of x and x 2 ara obtalnad from (3.3.7). Tha taaulta 
are 

2 

<x) . - 2* , U. 4.13a) 

4a^ 


<x 2 > 


• 1/2 + 



b) 


Thus, the average equation of motion of h(t) after changing the 
independent variable from t to the dimensionless parnreter c(t) isc 


fjf * - « - 21k) 

de 4. 2 

Integrating the above equation and choosing the integration constant 
so that h « kc when c m 1/26, we find 


ii HeriTi i ir > - ' 


*44.. .i** I H 
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h(c) » Cl 



<► Ikfi 


(3.4.15) 


Going /jack to the relation vhich defines sin* in terms of h 

m 

(3.3.9.17), we find that 


sir* 

m 



♦ Xk<2 - % 


0.4.16) 


Am long as 1 is small, the effect of the x-dependent in the tidal 
torque Is minimal. 

At the opposite extreme, the magnitude of 6 for the Mila- 
no* 

Tethys case Is % iO* 4 . This fact, coupled with its presently 
large libration amplitude (97*) indicates that 3 q #as also small at 
transition. For this case, a different approximation can be invoked 
to determine its secular behavior which depends on the fact that the 
fluctuation in x is always small, compared to one, if 8 0 is small. 
This follows from the observation that the maximum fluctuation in x 
is of 0(|s| X/2 ) in the small fluctuation limit. The next step is 
to consider the time average of ♦ in the libration phases 

<*>■-<*)- <c) - <b ( (x)cos^ (3.4.17) 

This time average identically vanishes in the abeence of a tidal 
dc 

torque and must approximately vanish in the adiabatic sense if 
the changes induced in b^Cx) are slow compared with the libration 
frequency. Therefore to 0(|e| X/f *) we find, 

- c(t) 



(3.4.10) 
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Since the fractional fluctuation in hCx) is small, x can be 
rfplaced by its mean value Cx) in the function b(x). The resulting 
Hamiltonian is then Identical to example II of simple penduJum plus 
constant applied torque discussed in section 1.2. Therefore, in the 
libration phase, H(x,+,t) for the k - +1 e-type resonance 1st 

H(x,*,t) - 1/2 (x ♦ c(t)) 2 ♦ 6 (c (t) ♦ l) X/2 cos* , 

While the equation of motion of H imt 

♦ c(t)> + 1/2637(0(6) ♦ 1 ,-VW ( 3 . 4 . 19 ) 

dt dt at 

Observe that the dependent variable can be changed from t to c (t) . 

We can either take the time average of the above equation or 
use the action integral to determine the secular behavior of H with 
c(t). The action integral represents the simplest approach for 
this example. If the x-dependent term in c(x,t) is included, then 
it appears that the averaged equation of motion must be used and 
that numerical integration is required to fir.d H as a function of 
c(x,t> (Allan, 1969). 

The evaluation of the action Integral in the libration phase is 
a straightforward exercise and, rather than repeat it here, we 
shall ol srve that it agrees with the result of Best (I960)* 

J... '• 46 (c + 1) I/2 [B(K) * a - Jt 2 )K(K») (3.4.2C 1 

lib. 


Where the parameter K is given by 



193 




196 


K(K) and E{K) art tha complete elliptic integrals of tha first and 
aacond tvpes, and ara dafined as follows* 


Tl 

K'K) -J (X - K 2 «ln 2 0 )"*d 9 
o 

TT 

7 

e (k) - J (i - K 2 ein 2 e> + *ae 

O 


(3,4.21a) 


b) 


Tha most interesting properties of tha above functions ara tha 
following (Byrd and Friedman, 1971) t 


K(0) - E(0) - y 
E(l) - If 

lira K -» 1 K < K * " 

Tha amplitude of llbratlon 4, as a function of c(t), is 
a 

obtainad by evaluating H whara i vanlahas. Vfa find that 6 can be 

a 

directly related to tha parameter Ki 


(3.4.22a) 

b) 

c) 


0.4.23) 




Using (3.4.20.23), $> n is graphed in Fig. 3,4.1 as a function of 

a (c(t) ♦ l) 1/2 . 

In the limit of small Iterations the simple pendulum can be 
approximated by a harmonic oscillator. The adiabatic constant for 
a harmonic oscillator is proportional to |b(^x) ) | 1 ^ 2 , or 

ID 

* « (Ctt) ♦ l)' 1/8 i 6 « 1 (3.4.24) 

m m 

This ooana that adiabatic damping in the small lifaration limit (when 
! 3 1 « 1) is an extremely slow process compared to the adiabatic 
damping found for the approximation obtained when 1 3 1 » 1, For the 
Kjrwis-Tethys case we find that 4^ eq. » ) 30* wnen the function b(x) 
was approximately five times laroer th. n at present. This implies 
i>.*& the dimensionless parameter B Q was approximately 5"’ times 
larger than 8 ^ ^. Including the moment ksd dependence in the tidal 
torque increases the rate of damping as a function of the 
inclination of Mimas such that the function blx) was about three 
times smaller than at present (Allan, 1969) . In any case, the 
tidal evolution of Kimas-Tethys can be adequately determined using 
this approximtion since the implied magnitude of $ o at transition is 
of 0(10~ 2 ). 

Again, we should observe that this contribution to the damping 
from the dissipative term is magnified in the Mimas-Tethys 
resonance by the unusual situation that the ratio of the tidal 


>><<»>> 

b(0) 


torques of Mimas and Tethys approximately obey the same ccanensur- 

ability relation as the resonance variable. The effect of this is 

to greatly magnify the ratio of the x -dependent term and the 

effective torque, acting on the resonance variable. 

If $ q happens to be of 0(1), then no approximation appears to 

be available. It happens that for the Titan-Hyperion resonance, 

which we shall examine, has 8 • 0.058 and 4 (now) - 36°. To 

now m 

treat this case we shall use the action integral to determine the 
secular behavior Of the system as a function of c(t). The action 
integral in the libration phase takes the following form. 


J lib. " * -k** 1 /(-kx ♦ l)d* + k" 1 ^* (3.4.25) 

We know that (-kx ♦ 1) is positive definite since e(x)>»(-)oc ♦ l) i/2 
and e(x) is real. For a libration, the initial and final values of 
4 are identical. Therefore the term fd4 vanishes. A more explicit 
form for the above integral 1st 


-lib. 




♦■in'* * 0 


(-kx ♦ 1)4, ♦ kj 


W* ‘ 0 


(-kx ♦ 1)4, 


(3.4.26) 

Since the Hamiltonian is symmetric about 4 • mr (n being an 

integer) , the roimiaua and m ax i m u m amplitude are equal and opposite 

(i.e. ♦ - * ■ -4 v ) • Corresponding to 4 equal to * , the 

a max min » 

action variable x equals * where x_ is defined by the condition 
m td 
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i * 0 or 

2 d 

m ~*m * a ~ ^2^* Wc * ♦ l)* ,k/2 co«^ (3.4,27) 

If the libra t Ion center is at * • mod(2v), then x lies 

n 

between the two roots x^l* > 0) and x^C* < 0) . We can change 
the integration variable f a * to at in (3.4.26). The result is 

3{ X 

J IU>. ' •*/ " ( " 1 “ ♦ ** + 2k J 2,,+ (-lot + l)-r^r<l* 

* 2 .- ' 

0.4.28) 

Th. function can ba obtained from i||l and expressed 
as a function of just x,H,3 and c (see B.6). 

The roots x 2f _ and can be expressed as functions of 

H,B and c. (App. C) . The Unit ran be expressed in terms of 
c and H by using the Hamiltonian evaluated at 

a - 1/2 tx, e c) ♦ ♦ l) V2 cos* a 0.4.29) 

to elininate the coe*^ dependence occurring in 3,4.27. From these 
two equations* we find that p fl U ^ + c) satisfies the following 
equation! 

K(1 - - Ike + l)p n + h » 0 (3,4.30) 

The solution for he function (-kx^ + 1) for the |k| - 1 case 1st 


Mw n ♦ 1) » 1/3 (kc ♦ 1) ±j/~0tc ♦ l) 2 ~ (|)H (3.4.31) 

In fig. 3.3.7, there are two distinct angles, and for which 
4 vanishes. Observe that as the system evolves towards transition 
into the libration phase, the angle ** increases to 90*. Thereafter 

IQ 

4* no longer corresponds to a value for which i vanishes. On the 
m 

other hand, 4 does correspond to a real libration amplitude in the 

iff 

libration phase. The important question 1st Which solution (i) 

corresponds to this "normal * solution? The right hand side of 

(3.4.31) is positive definite and the above equation must be valid 

for all possible values of 3. In the limit 1 S ] >> 1, the function 

c (t) , at transition, equals 1/23 which is presumed to be a large 

negative definite number. Therefore, for the k • +1 case, we 

should choose the ♦ sign in (3.4.31). In the limit |b( << 1, 

x m w -c(t) in the libration phase. Again the + solution of (3.4.31) 

agrees with the expected behavior. The "nonaiX" solutions for 

(-x fc + 1) and cos* arei 
m m 

(-x^ + 1 ) . 1/3 (c ♦ 1) + //9<c ♦ l) 2 - <2/3)Hi <3.4. 32«) 

won, - 2p a (-* a ♦ b » 

The action integral can be solved in terms of standard 
elliptic integrals, but its form is exceedingly complex. A simpler 
procedure is to numerically evaluate the integral for a given value 


rffrr 


. . -A »*• 
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of c and a test value of H. Then H can be varied tip til the 
numerically calculated value o t agrees with its appropriate 

initial value. J lifa is easily calculated at transition and is 
found to be (B.30,31); 

j . { 4»e > iBl < lB cl l 

'’lib. X -2vk , |0 1 > |0 cl | 

Figure 3.4.2 is a plot of the initial value c(t^) versus the 

magnitude of 8. The break in the curve occurs at 8 cl * -0.2722. It 

results from that fact that transition for [b| < U ci l involves the 

coincidence of two e-roots while for | B I > U cl U it involves the 

vanishing of b(x). Figure 3.4.3 is a graph of si,i$ versus the 

m 

parameter (1 - 2c{t)/0' for several values of the parameter 
|fl| 5 i® cl l- Thi 8 graph clearly supports the analytic 
approximation developed in the limit | &f » 1. Note that the initial 
slope of the curves approaches the straight line generated by 
ploting 2c(t)/0 versus (1 - 2c(t)/0). Furthermore, the value of 
c(t) for which the slope begins to flatten out is approximately equal 
to Figure 3.4,4 is a graph of 6 versus (c (t) + l) 1 ^ 2 for 

IQ 

values of |s| * 2. The important observation is that the curve 
rapddly approaches the limiting curve generated by the approximation 
obtained when |b| « 1. 

Zn the course of this development , various kinds of «>. ilxlity 
have been mentioned or implied. In studying the neighborhood of 
stationary solutions of the Hamiltonian (Section 2.7), we .inferred 
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that these solutions may or may not be stable against small 
perturbation*. If dynamically stable, these stationary solutions 
were designated as libration centers. The system of simple 
pendulum plus constant applied torque suggested that librations were 
possible only if the magnitude ot the applied tidal tor rue is smaller 
than maximum value of the pendulum torque or is H tidally stable." 

Finally, in investigating example III in which the coefficient 
b(x) of the pendulum term was momentum dependant, we found that if 
the system made a permanent transition into libration, the 
magnitude of b(x) had to tend to increase thereafter. Since it is 
the s)ow change in the parameter c<t) which indirectly causes the 
magnitude of b(x) to adiabatically increase or decrease, we shall 
use the terms "adiabatically stable" and "adiabatically unstable" 
in alluding to this kind of behavior in the next chapter. 



FIGURE 3.4.2 


Graph of the initial value of c at transition, c , versus 

o 

the magnitude of 6 for ] B | ,^ 2. 


fV i >r*i >»*A. .aSafoam v**\j*’ '*** 




FIGURE 3.4.3 


plot of the sine of the amplitude of libration versus 
the parameter (1 - 2c ffi) for several values of & above 
the critical value of B cl - 
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0 5 10 15 20 


[c(t) + l ] ,/2 


FIGURE 3.4.4 

Plot of versus the parameter (c(t> + l) 1 ^ 2 for values 
of |0| *■ 2. 
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4.1 THE SATELLITE- SATELLITE RESONANCES OF SATURN 

There exist three examples of two-body resonance interactions in 

the ten-satellite system of Saturn. Observations of Mimas 
*Mi 

(- — . 3.11) and Tethys (4.94) reveal that the ratios of their mean 
motions is 2tl and that conjunction of the two satellites tends to 
librate about the midpoint of their nodes with amplitude 48.5*. An 
investigation of their mutual gravitational interaction (Tisserand, 
1896) shows that this phenomenon can be explained as a gravitational 
resonance in which the <*ngle $ = 4\ - 2A_ ♦ fl • lihrates 

Ml Te Hi Aft 

about zero with maximum amplitude 97* and with a pe*i.od of 70.78 
years. The coefficient of the term in the expansion of the disturb- 
ing function specifically responsible for the observed behavior is 
proportional to the product of the inclinations: I Mi I Te * Thus the 
resonance can be classed as a mixed * type (2.3). 

The con-functions of Enceladue (a. 99) and Diene (6.33) are 
observed to librate about the per ic enter of the inne** satellite with 
period of approximately twelve years (4 £ X gn - 2X D 4 The 

amplitude of libration is very small, quoted values rarging *rom 
20* (Goldreich, 1965) to 1.5* (Sinclair, 1972). 'jfc fr.sof 
the resonant perturbations in the mean ? ongltu-.es of u \ir are 
even smaller j 1.4* in Encelodus and 0.9* •*». Ji one -Brouwt id 
Cl erne nee, 1961, p. .193) . This means that the librati- the 
resonance variable is principally governed by tne -i of the 
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per leer ter of Eaceladus. 

The pair of sate!' iw far thee t rsoovod from Saturn is also 
involved in an e-type resonance. In this case, conjunction of Titan 
(20.48} and Hyperion (24,83) libra tee about the apoc enter of the 
outer satellite with amplitude of 36* and period equal to * 2 years. 
The cocnaenaur ability ratio is 3x4 and the resonance variable let 

* 3, « - 4 *«y + V- 

One novel aspect of this resonance phenomenon is that it tends 
to Keep the satellites as far apart as possible at conjunction. This 
behavior is seen as enhancing stability among the participating 
satellites. We should also mention that the masses of the satellites 
can be determined from a knowledge of the p riod of libratlon and the 
ratio of the libration amplitude of the mean longitude of each 
satellite (Jefferies, 1953). 

This is the s«*t of information available concerning these 
satellites. Any speculation concerning the existence of appreciable 
tidal ly- induced torques acting on any of Saturn's satellites has 
not as yet been supported with visual evidence of either a secular 
change of orbital periods or, in the case of resonances, a dis- 
placement of the center of lihration **ay from mod (v) or mod (2*1. 
Goldreich suggested that tidal torques are acting on the sa> 
and offered two arguments to support his thesis: 1) The exi 
so many resonances cannot be a chance affair and must bo due tc sox*, 
mechanism. He suggested that significant tidal evolution of the 
inner satellites of Saturn and Jupiter must have taken place over the 
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g 

age of the solar system (% 4 * 10 years) . A lower bound on the 
dissipation function q" 1 was calculated by integrating (2.9.4) and 
demanding that the closest satellite (ftimas) was at the planet's 

n 

surface 4 x io years ago. The tidal dissipation function Q is 
defined by the relation 



where E is the maximum energy stored in the tidal distortion a. ' 
o 

the integral is over one complete cycle. 

2) A theoretical alculation of the dissipation function was 
attempted for Jupiter and Saturn and was found to be in rough 
agreement with its lower bound (Goldreich and Poter, 19661. Q 1 
for Saturn is estimated to be 1.5* x 10^. Accepting this estimate 
of Q and assuming that it applies to all the other satellites, we 
find, for example, that Titan, wnich is the most massive of Saturn's 
satellites, has increased the radius of its orbit by only 1/4% over 
the age of the solar system! 

We should mention that ving -his value of Q and (2.9.4', the 

-22 2 

tidal deceleration of tie mc*n motion of Mimaa is 1.4 x io sec , 
or equivalently, 0.04° century The magnitude of the romance 
torque can be estimated from the period of libration Ci e.» 70.78 its. 
for Mi-Te ) • We find that the ratio of the tidal terque acting on 
Mimas to the parameter B is ^ 10~ 5 . Becaus the torques acting on 
each body tend to cancel in the resonance variable, the ratio of 





b(x) to the sun of tidal torques in $ (2.9.10) la an order of 
magnitude larger. Sine* the torques acting on the other satellite 
resonances are weaker , the ratio Is even larger for the ? systems. 
Therefore, the hypothetical tidal evolution of these satellites 
should be well described by toe theory developed in chapters one and 
three. Recall that its quantitative accuracy is set by the parameter 



As mentioned e. --tv the effect of dissipative tides raised by 
a given satellite on its orima.-y is to cause a torque parallel to its 
angular velocity. This is true if the spin of the planet and the 
orbital motion of its satellite are in the same direction anJ the 
planet's rotation period is shorter than tae satellite's orbital 
period. This torque tends to increase the size of the orbit and 
decrease its period. One expects that in a many-satellite system, 
after a time, some pairs will approach a ccsmnensur ability and then 
evolve throrgh a succession of related resonances. Table 4.1.1 is a 
list of the strongest resonances associated with a 2:1 
cammensur ability. lUnprimeJ variables refer to the inner 
sat*'lite, primed variables * - o the outer body.) 

They have been ordered in the same sequence in which the pair 
would ncounter them under the following assumptions i 1) Tho tidn* 
acceleration of the inner satellite is at least twice that of the 
cuter 3ne. That the tidal torque acting on the inner satellite 
determines the sign of the torque acting on the resonance ~*ariable. 
Therefore, in the absence of a resonance $ tends to deertase 
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2) The motion of the perihelion ui is prograde while that of the node 
ft is retrograde. This kind of behavior is caused by the secular terms 
in the disturbing function. 3) The motions of 8’ and ft* of the 
outer satellite are smaller than the corresponding motions of 0 and 
ft. 

We should observe that all the resonances listed in Table 4.1.1 

are "adiabatic ally stable" U.e., tht .. lal torque tends to sec .i&rly 

increase the coefficient b(x) in the libratio.* ohase and decrease 

the llDration amplitude). This pathological result is related to 

the two body resonance. We have not rejected any tidally unstable 

resonance variables cutright. The strongest tidally unstable 

resonance must have at least three leading factors of e and I. For 

example, the angle $»X-2X-w + 2w' has a leading factor cf ee 1 

It is adiabatically unstable if, in the libration phase the tidally 

induced decrease in e is greater than the corresponding increase in 
o 

e* . In the three kx>dy interaction, both kinds of resonance 
variables appear with comparable coefficient-. b(x) (see Table 4.2.4 
for specific examples) . 

Table 4.1.2 lists the: pertinent data now available on the masses 
> ,& orbital elements of these satellites. This information shall be 
called upon during the course of this discussion and is collected 
here for convenience. 

Table 4.1.3 lists the observed angular frequencies of these 
satellites, the libration period of the resonance and the observed 
periods associated with the secular motion of w arvd ft. The 



observed secular action of A and ft includes both the contribution 
from tho •ocular and roaenanoo terns of tho llaturblng function, 
explicitly 


Loading Dependence 
of b(x) 


Reeorv.nc# Anrle 


i 



I — 



* 



2X-4>'*2fl' 

2h-4X*oD>ft' 

2X-4>'e2tt 

X-2X’*tf 

x.ax^ci 

2X-4V+2*' 


2A-4XU2& 



TUU 4.1.1 

table of tho wo jot coso nance frequencies aseociatod with a 
2»1 conaonour ability of the inner (inprimed) satellite with 
the outer (primed) one. The fre*w<mclee are ordered in the 
•atse sequence in which thie pair of satellites encounters 
then. 


A& 


ob 



♦ <**> 
* w 


rss 


(4.1.1) 


The separate contributions from the sseular end resonant parts of 
the motion of the appropriate per icon ter or node have been calculated 
by Jefferies (1951), Zn rone cases data la missing either because It 
Is inappropriate or because the corresponding eccentricities and 
inclinations are vary mail and variable and tha corresponding 
observed avs rape motion of ths peri center or node has not been 
determined. 

Accept tr>x Hyperion, the sign of the motion of perihelion end 
node Indicate that contribution from ths secular term in the disturb* 
lnq function is much larger than tha* of the resonance term. 

These frequ ncies appear to be well spaced for both tha Mimes- 
Tsthys and JCncslaf us -Dions examples so that ths evolutionary picture 
involving a single resonance variable tentatively applies. Zn the 
eeee of Hyperion, the motion at the perihelion is retrograde and 
lares compared to tha expected prograde motion due to the secular 
tenia of the disturbing function, Compare ^ ** ♦ O.l'/yr, with 
||t(y * * lfVyr. both should be approximately the same if the 
prlneipel contribution to the motion of fl came from the seedier terms. 
The preeont ordering of the resonance frequencies for ‘-a Titan 
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Hyperion example for the 4*3 comaeneur ability indicates that the 
first resonance variable encountered is the presently observed 
example. This is due entirely to the large negative (retrograde) 
motion of resulting from the resonance tens. Before we proceed 
we should ask if this sanw* situation prevailed at transition * or 
the other two examples. 

In the 1 ibration phase , an ediabaticaliy stable r* 00 nance tends 

to cause a retrograde notion of either perihelion or node, depending 

on whether the angle variable is an e or I type, respectively. 

This can be seen by inspection of the appropriate equation of 

motion for & or 0. In the case of an e-type (~) 1st 12.4.3* 

at res 

3.1.7b): 



If 4 librates about mod(2ir), then B is negative. For an adiabatically 

stable resonance, k is positive, hence (tt) is negative. The 

Qt reo 

angular frequency associated with this impressed notion, in 

dirtens leal ess tine unit t (5 A x t), is of 0(b). In terms of reel 

oxx o 

tine t, this angular frequency is approximately given by* 
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a n 2LL{-£) c#' 
dt res n o p V 7 c 
o > 


{|jc| - 2) 


M.1.3) 


We see explicitly that the inpreseed notion of 0 tends to blow up as 

e o * 0 only for ths |k| • 1 esse. If e Q is very null st transition, 

then the notion of nay be both large and retrograde such that 

the first resonance variable encountered is the k • *1 e-tyvS. 

The mixed- X type can also lead to a large retrograde wtion 

of (ft ♦ O') if one but not the other inclination la very small. Bvei. 

so, we expect that the impressed motion of the nodes is 0(1) smaller 

than the elmiler motion of the perihelion of the lighter satellite. 

This means that the e-type <k - *1) resonance variable can still be 

the first sneountered unless is vary much anallsr than s . 

O o 

from Table 4.1.2 we can see that the motions of ft and ft, in 

which the secular term predominates, tend to be equal and opposite, 

and that the retrograde motion of the node of the inner satellite is 

greater than that of the outer one. (Spoclfi.ally, Inspect the 

Mlmas-Tethys case to assure yturself that the above statement is 

true.) if this is the cate, we should determine the critical value 

of e such that the impressed retrograde notion of 0 (or for that 
o 

matter, S') ie greater than twice the retrograde motion of the node 
ft of the inner satellite. 

The present valu* of the eccentricity for Mimas ie 
exceptionally large compared to that for other inner satellites 
(Compare e^ » 0.020 with the next largest eccentricity found among 
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satellites is presently captured. We shall discover la tar that a 

••4 

reasonable eetimato for a at transition la * 10 , 

°En 

We now see that tha earlier description of tidal avolution 

through an ordarad aaquanca of well-spaced raaonanca variables only 

appllas If tha eccantricitias of either partner of the resonance are 

not too small. Otherwise, it may happen that a variable such as 

<1 - 2X * ♦ u)} is anoountarad first. Whether or not capture occur e 

shall be our next topic . But before we proceed, we should observe 

that the one dimensional model so carafully constructed may fall 

under certain circumstances . 

Consider the implication that the relative order of tha 

raaonances may be interchanged depending on tha paramatara a , I , 

o o 

etc. It may happen that two resonance frequencies may nearly over- 
lap. It is no longer necessarily true that tha resonance system can 
be described by a one dimensional Hamiltonian. Normally one expects 
that an e-type resonance (i.e. b(x) - e) is stronger than, say, a 
mixed I type <i.e. b<x) « II* )• Then, hopefully, a reasonable 
approximation would be to ignore the mixed I type. But if e^ is very 
small, then the dominant resonance may be the mixed 1 type. 

Another possibility is that e Q is so small that the corresponding 
resonance variable is well -spaced from any I type resonance a t 
transition . If this value is less than the critical value, then the 
system will automatically enter the libr&tlon phase with an initial 
amplitude of 90*. Since the resonance ie tidally stable, the mean 
value of tha eccentricity will thereafter increase. Aa the mean 
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value of e increases, the magnitude of the impress**! retrograde 
me>t*on of the perihelion must decrease, eventually it will overlap 
with an I- type resonance. It may happen that this T. type disrupts 
the established resonance, depending on their relative strengths. 
Qtactly what may transpire in either case would require a rigorous 
examination of a two-reaonance variable system subject to a constant 
applied torque. 

Does either case have relevance to the previously discussed 

examples? In the Mimas-Tethyn case, compare the value of e for 

°Te 

which these two var:*ble« overlap <4.1*4) and the present value of 
I hi I Te* We * ind that ***• lattor is of 0110) larger. For this 
example, it appears that if the two resonances overlap, the mixed I 
type still predominantly determines the fluctuations in the mean 
longitudes, X and 1*. On the other hand, if we make a Similar 
comparison of the variables ( X * + and 

(21^ - 4? Dt ♦ ♦ fl^}, w * find that the coefficient of the 

©'•type is 0 (10) greater than the mixed I type. 

.e 

It seems unlikely that e is or was ever as email as 10 . If 

°Te 

this be true, then Mimas-Tethya evolved through the sequence found 
in Tablo 4.1.1. The first resonance encountered was the 1 resonance. 
Presumably because of unfavorable initial conditions at transition, 
the system evolved past this resonance and later approached the IT 1 
resonance. The corresponding resonance term in the disturbing 
function R has the form* 


j| ' ' v ' j •; • * ’■ l ***e*d> a*.*** W ‘W* -r . l * *-*-*►» -< ■ s** 
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{R) r... ■ ^7 I M1 U)1 T. U1C “**' M * l ‘« 

where the dependence of the inclinations on the action variable x 
is (2.4.4; 2.6*1) t 


Jttii. . <-* . i,va , 


^ ■ T * ( - - (^> * * i) l/ ' 2 ( 4 . 1 . 7 ) 

I 1 

The mass ratio * Thus the variation in I_ tx) with * is 

nu. i ' *e 

1 

(j^>) times smaller than a corresponding change in X^. This also 
applies to any secular change in the inclinations after libration 
is established. The factor C can be expressed in terms of Laplace 
coefficients (2.2.15) and numerically evaluated (Tisserand, vol, 14, 
p. 100, 1896)t 


C = - iabj^(a) ■ - 0.+ OHfe. 


(4.1.8) 


The probability of capture is determined by the dimensionless 

parameter f which occurs as a factor in the pendulum-like term of 

the Hamiltonian (3.1.5), The appropriate parameter in this case is 

(obtained by 4evi4lngoy Z»- i/2 1^ L )i 

°Mi 


d. 



+ il 



Compart this with the critical value 0 ffl * - 0.27 for which the system 
automatically enters the libration phase. Xncidsntly# the 
Inclination of Mimas would have had to toe ?1 for this to occur. * The 
present libration amplitude ;97*> and the small value of B 

now 

indicate that 0, evaluated at transition# will also be small. 

Therefore the approximation developed in section 3.4 in the limit 

| P o | << 1 can be applied to this resonance. 

Neglecting the effect of the x- dependent term in the tidal 

torque (2.10.6)# we find that transition occurred v**an • 0,20^, 

or when the mean inclination of Mimas wae one-fifth its present 

value. Allan (1969) included the effects of the x- dependent tern and 

determined the value of the sereixajor axis and the inclination for 

each satellite when 4 m 180* numerically. His results are* 

3D 


,- 3 ' - .**.•< 


..j 'Z*r »v>v~ 
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a » 0.9 922 a (now) f 
°Mi Mi 

I - 0.277 I^Cnow) - 0,«415»6, 

(4.1.11) 

a • 0.9922 a- - 0.9922 (now) 

°re T * Tt 

I - 0.945 1-Artov) - 1?046. 
o Te 

Ta 

Calculating $ using the above values# we find 

0 » -4.4 * 10~ 3 . (4.1.12) 

o 

Since the value of the parameter 0 Q is quite urall compared to one# 
(3.2.12) can be used to approximate the probability ? c for capture 
into libration. 

# c * i|s|* A/2 b x (o> » lu ♦ y^-) I s 1 1/2 * «•«*• M.i.ij; 

Net. that b(x) i. proportional to (-* + l) 1 ' ,2 (-y^* ♦ 1> 1 ^ 2 and that 
the mass ratio of the inner to the outer satellite is * jy. Thus 
b^tO) eauals 1/2(1 ♦ I e l whete & "P* 1 * V Sinclair found 
through numerical calculation that '« 4%, which is in agreement 
with the above analytic result. The probability for the first 
resonance encounter can be estimated by comparing the appropriate 
function b(x) for an I 2 resonance with that for an II* type. For the 


k ■ +2 case btx) 1st 
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C* COS* . 

*Di 80 


C4.1.16) 


From Jefferies (1953 ) , th« factor C approximately equals - 0.753 
while the coefficient C* belonging to R' equals ^ -1 H . The 
relevant parameter B for this type 1st 


v -e- A oxx a Di a En 


3 **■* 12 

« 2 En ^ 4l 


The parameter \ is obtained from 2.6.10 and equals 1.59. The mass 

» . 

ratio of the satellites is * Calculating 6 we find 

bl. now 


8™. —10. (4.1.18) 

now 

Since the critical value $ cl * - 0.27, and |B q | must have been 
larger than | ® rvo%r l » w ® can conclude that the system automatically 
entered the libration phase with maximum amplitude of 90*. Since 
1 6^1 >> 1 # the relevant approximation of tne Hamiltonian applies 
(3.3.3). Using Sinclair's quoted value for the libration amplitude 

(4 “ l.’S*), the change in e since transition is obtained from 

m- *• 

(3.4.1), We find 


®o 


e 


now 


sin* 

m 


~ e a 1,1 x 10 
42 now 


(4.1.19) 


The relation between the initial and the present values of the 
parameter c(t) occurring in the Hamiltonian equation (3.4.4) is 
given by 
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- . -± c 

~now 42 o * 


(4.1*20) 


The change in c(t) is directly related to the change in the 
comcensur ability relation by (2.9.9). Explicitly* 


OSO€ O 


diU diL 

"‘Sn ' *\i ' T( ^L' J %7 ) ' 


(4.1.21) 


T is tb« time since transition. Of course c(T) * c 

now 

The change in the comnensurability relation can be expressed in 

tanas of tije present average value for e_ <e - 1/2 (e + e , )) 

en flV8 • mx min 

using (3.4.0, 4.1.20). Namely, 


n Bniu^*a? >C *Bn <noW ^ ' 


(4.1.22) 


where we have assumed that the mean motions (and samimajor axes) 

have changed little since transition. Observe that the right hand 

side tends to vanish as e_ increases. The implication is that the 
En 

approximate coranensurability tends toward an exact one as th system 
evolves. Goldreich was apparently the first to make this 
observation (1965) . 

The next step is to determine how .much the mean motion and 
the semimjor axis of Kneeled us have changed since transition. 
Assuming the same dissipation function for both satellites we find 
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(4,1.23) 


Since the resonance only weakly affects the motion of the mean 
longitudes in the limit |fi| » 1, we find, for example, that the 
change in the mean motion of Enceladus since transition, An^, is 
the following. 


— S.edctTJA^t-y * - 0.032 n^m (4.1.24) 


where Ac(T) equals c(T) - c ■ 41 c (T) . This corresponds to a change 

in the semimajor axis Aa^ of Enceladus given by (using the 

approximate relation* — • - ~> i 

a 3 n 


^En * °* 021 a En lT) * (4.1.25) 

Using Goldreich ’e estimate for Q we find that the transition into 
90® llbration occurs whan 

T ■ 1.4 * 10 9 years. (4.1.26) 

Again, this appears to be an estimate within tha age of the solar 
system. 

Titan and Hyperion are also presently engaged in an e-type 
resonance. The relevant part of tha disturbing function is* 



(4.1.2T) 
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Xncidently, the value of e^ corresponding to 6 cl • - 27 jss 

e t - 0.057. .1.29) 

Cl Hy 

The above value for 0^^ ia certainly close enough to the critical 
value for automatic transition that we should expect that it did 
occur- This node of transition can be inferred from the following 
argument. 

Assuming for the moment that the limit j 3 1 >> 1 can be applied 
to determine the mean value of the eccentricity at transition, we 
obtain the result 

e * sin (4 ■ 36. 2* Jaj, (now) - 0.061 (4.1.30) 

°Hy * niy 

- e_ (now) sln(4tj “ 36* ). 

Hy m 

This limit represents the fastest possible damping of for the 

smallest possible change in the eccentricity of Hyperion. The 

actual change in e^ since transition must be greater than the 

above value. Thus this value (4.1.30) represents an absolute 

punHimiffl for the value of e at transition. After comparing 
°Hy 

e , with (4.1.30). we can infer that capture into 90* libration 

C V 

occurred. Unfortunately, we also see that neither approximation to 
the Hamiltonian (i.e. |B| » 1 and |s| « 1) can be rigorously 
applied to determine the evolution of the system. Brt we can 
attempt to match the solutions found for the two limits to obtain 
an estimate on the evolution of the system since transition. The 



most reasonable value of 8 to match the two solutions is 6^. That 
la* we shall demand that the system evolve backward In time according 
to the solution obtained fox |8| « X until b(x) - 8^. Before that 
time have this resonance evolve according to the solution found in 
the limit |s| ?> 1. A more rigorous derivation using the action 
integral is given* following this argument. 

The present values of 6 and # are 0,0575 and 36** respectively. 
Going back in time* we find that when btx) * B ci at time T, b(X> : , 
has decreased by a factor of 0.21. This means that* since b(x) is 
proportional to e^ (x) , the value of e^ at time T is 

• Q.211e (now) * 0.022, (4.1.31) 

ny By 

Proa figure 3.4.1, the 1 Duration amplitude at time ¥ is 

? - W*. {4. 1.32) 

a 

The relation between e * e {now) and the parameter a (T> is 
ny By 

r a<n0W> - lam +l) 1/2 « (a«)> 1/2 , (4.1.33) 

e Hy 

where we have chosen a <T) *♦ 0 at time T. Since (now) is much 

Hy 

larger than e Hy , (a(T) + 1) *» a(T). 

The maae of Hyperion is much smaller than that of Titan, 
implying that the resonance has not affected the tidal evolution of 
Titan. Therefore, the change in the mean motion of Titan since 
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FIGURE 4.1.1 

Plot of e and sin6 versus ) for the Wtan- 

ave m d 

Hyperion resonance. The 8 dashed line indicates the 
present values for these parameters. 


- >rr vrT^'Vr^ 4, -i **«&♦ W* ’«* 


Rather surprisingly, the above value for e agrees reasonable well 

Hy 

with Greenberg's (1972) estimate of 0.015. The change in Titan's 
semimajor axis up to the time T, a . ’ the value for T, can be 
determined with arguments similar to those applied in the Enceladus- 
Dione case. The result iat 

Aa Tl a 0.0125a Ti f T a 2 x lo 10 years. 

Thus the age of the resonance and the total change in Titan's 
se amajor axis are t 

da Ti a 0.07a Tl (now>; T a g x io 10 years. (4.1.40) 

Greenberg estimated that transition occurred 4 * 10 10 years ago. 

A more accurate estimate of both the age and the change in the 
orbital parameters can be obtained via the action integral. To 
make use of the action integral, we must first reduce its 
dependence to a single unknown parameter, if we evaluate 4, H and 
e at $ equal to • 36° where $ vanishes at time T (corresponding 
to the present), we find 

♦La - 0 - ** ' c (T) + 1/2BC-* + 1) 1/2 cos* . 

v m y m m m 



(4.1.41) 
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The change in the parameter c (t) since transition »s related to the 
change in by 

in 7i - i(c(T) - c o )A oxx r o 

' % * "n * 

The corresponding age T is 6 * lO 1 ^ years which is in reasonable 
agreement with our earlier heuristic calculation. The evolution of 
the parameters as a function of c(t) can be found using the same 
procedure as outlined in (3.4). Figure 4.1.1 is a plot of the average 
eccentricity (h 1/2 “ e^^)) and the amplitude of libration as 

a function of the parameter (1 - 2c(t)/8). 

The reasonableness of these calculations is naturally 
conditioned b* the initial ac sumptions. Clearly the most crucial 
of these iu the dependet*.e of the tides on a constant Q which is 
the same for all satellites. If Mimas could have risen from the 
seas of Saturn , Fnceladus and Titan would have been nearly motion- 
less spectators to the event. Only a 6% and a 1/4% change can 
occur , respectively, in the orbits of Enceladus .nd Titan. Thus 
"significant" tidal evolution is limited to the closest satellites. 

Perhaps one way out is to say that Q has an amplitude 
dependence (let's guess that it's proportional to the height of the 
tides raised on Saturn's surface). Because of its greater mass, we 
might expect that the Q for Titan is significantly larger than it 
is for other satellites, rhe amplitude of the tide is roughly 



proportional to he gradient of the force at Saturn's surface. 
Explicit: 71 

Tide Height «>;{£) 3 . {4.1.42) 

Comparing the tide height raised by Titan and by the next 
strongest case, Tethys, we find 

Tide Height by Titan . 

Tide Height b - Tethys 10 (4.1.43) 

frtxn which we might conclude that it is at least possible that 
Sitan 18 8l * nificantl y end that the age estimate for this 

resonance is within reasonable bouoc . But if Titan's tidal torque 

is greater, then why not Tethys* - The tide height of Tethys is 
about four times greater than that for Mimas. The problem is that 
the two tidal torques nearly cancel in the commensurability relation 


,_3f2s_ . 

20dt Hl 


with Mimas just barely winning the battle. If were just a 

'’Te 

few percent larger, the resonance variable becomes tidally 
unstable) 


We should mention that the probability fur capture Is 

increased slightly C* 1%) by _he x dependent term in ^(x,t>. But 

at 

it also happens that the ages and the capture probabilities are 
significantly affected by the assumption that both tidal torques 
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obey (2.9.4), independent of the value of £. 

In the Kimas-Tethys case the values of the orHul parameters 

was significantly affected by the fact that the x dependent term is 

an important factor during evolution. The relative damping due to 

this term is large only because the two torques tend to cancel m 

I dn • 1 dn^ ; 

the cowmens irability relation (41.142). If 1 >> ] j* then 

the x-depe -lent term is much less of a factor. The evolution of 
the resc;iance is then accurately described by the solution obtained 
for the limit |o| << 1. Of course this means that the initial 
inclinations were less than predicted by Mian. Going back and 
repeating the calculation we would find that the probability for 
capture for the II* resonance is increased to 7-8%, while that for 
the I 2 resonance is increased to about 10%. Contrary to 
expectations, the at sonance is dec eased significantly 

by about one fifth. tuation holds in the Enceladus- 

Dlone case il 


I^T 

»|!^l 

i dt £n 

1*01' 


We've already noted that the dissipative term plays a minor role 

during capture and *volution in the limit | 8 | » 1* But the 

cancellation o* the tidal torques does affect the ccmnnensurabillty 

dn T 

relation in the original problem. However, if >> c j t * 

En Di 

the age of the resonance is decreased by a factor of 0.28. Thus, 
the age of the Enc ladus-Dione resonance ray be as small as 
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4 * I0 8 years. 

Although we’ve generated many numbers concerning the ages and 
evolution of the various resonances discussed, we find that good 
numbers are hard to find. This must be accepted as an excerise 
which demonstrates that Goldreich’s hypothesis is probably correct 
but that it generates more pussies than it solves. 
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4.2 THE LUNAR- PLANETARY RESONANCE HYPOTHESIS 

The muon, like the satellites of Saturn and Jupiter, is spiraling 
away from the earth due to a tide lly- induced torque. If a simple 
tidal model is invoked, and if the present value cf the t'dal 
acceleration is used to determine a constant dissipation function or 
Q-number, then several investigations have sttovn that the moon was 
within the Roche Limit less then two billion years ago (see &• reich, 
1966). As the earth-moon system appears to be much older, pome ..mg 
else must be invoked to resolve the time-scale paradox. Tht most 
plausible solution is that complex factors influence the tidal torque 
and that, contrary to expectations, the energy dissipation factor may 
have been considerably less in the past. Pannella, KacClintock and 
Thompson (1968) examined the tidally induced p iodic i tie j in the 
daily growm structures of various type? of shells of widely different 
ages. The implication they drew (fig. 4.1.1) is that the tidal 
torque was both variable and also probably less in the past, beyond 
approximately 70 million years ago. 

R. J. Hlpkin made the novel suggestion that pert<aps the time- 
scale paradox could be resolved if the moon were trapped for an 
appreciable time in the past in a resonance with Venus. Goldreich 
had already shown that partners of a resonance, subject to tidal 
torques, tended to maintain their near commensur abilities through 
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MftllOMft Of VfAII 


FJCUXB 4,3.1 

Plot of the number of days por synodic lunar month versus 
entiAeted ege of each sample. Tt is graph U taken from an 
article by Pennells, Haclintock and Thoprioun ( Science . 

Vol. 162, papas 792-96, 1968). Nasal and Waller determined 
that sons of tha oga estimates were in nrror, and this was 
acknolvedged by the above authors fSciatca . Vol. 164, papa 
201, 1969). The dashed Una la tha result of these 
corrections. 


i. J (rt'S . nv^l. 1 ' lO**,,*. ' oi >t ‘ - 
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a transfer of angular momentum from one partner to another. Therefore, 
e hypothetical lunar -Venusian assonance could transfer angular 
momentum from the lunar orbit into the much larger orbit of Venus. 
Hipkin reasoned thev *uch a process would negligibly affect the 
Venusian orbit, and would effectively trap the lunar orhjt at a fixed 
radius over the lifetime of tne resonance. He is presently working 
on this problem, basically using tha asms approach as outlined here, 
except that his solution is apparently much more precise, it should 
be mentioned that in hte solution (and here also) the planetary orbits 
are assumed to be both circular and copianar since otherwise the 
expansion and evaluation wf the relevant coefficients become 
Inordinately difficult. But in calculating the effect of second end 
higher order coupling, Hipkin does take into account the coupling of 
terms in the expansion proportional to the eccentric ity of me of the 
resonance partners before setting the planetary eccentricities to 
aero (Peels, private communication) . we shall find later that the 
approximation of circular planetary orbits effects the maximum value 
of the resonant torque. On the other hand, th* .^Ja of the tidal 
torque in the past is not well established either. Therefore, the 
first step should be to calculate the relevart gravitational torque 
with approximations which simplify ths calculation as much as 
possible, and c ampere ths magnitude of ths resonant torque with the 
present value for the tidal torque. If the tidal torque is much 
greater than a given resonant torque, then there is no need to 



241 

•urther refine tita calculation. The 1^11 cation la that tha tidal 
•volution of tha lunar orbit could net ba arrested by tha given 
resonant interaction. 

Venus la chosen aa the moat llkitly partner in any rasonanca 
because it induces larga perturbation* in tha maan motion of tha 

noon dua to lta relatlva nearresa to tha aarth. Although Jupiter la 

•« o 

much mere massive than Venus, a coun :er balancing factor of (—-) 

a > 

enters in the development of the disturbing function, where the 

integer c la the ratio of tho synod.: mean motion o* tha moon to tha 

synodic maan motion of tha disturbing planet. The ratio c is • 14 

for Jupiter and • 20 for Venus. A trivial calculation shows that 

^-t.2) << *-*(.72). A plausible partner, not considered by 

»© 

Kipkin, ie Mercury, which has the smallest ratio c of any planet 
f- 4) . Hipkins* original argument, applied only to resonances of tha 
synodic type, bet as ve have seen, the simple e-typo should also be 
considered, because of captuze cc ns i derations if for no other reason. 

The method followed xn ttte:'minlng the first and second order 
contributions of the disturbing fuix tions acting on the moon is the 
procedure outlined in sections 2.2 ind 2.5. Naturally, the first 
step ie the expansion of rive relevant disturbing functions. They arei 
Disturbing Function Acting on the Moon by a Planet. 




U.2.1a> 
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Disturbing Function Acting on the earth by a Planet. 

RfFitl-Jlptl - b) 

r p 

Disturbing Function Acting on tha Moon by the dun. 

1 

Figure 4.2.2 shows all tha rslsvant radius victors and anglss. Tb 

simplify tha expansion, we shall make th» approximations that all 

orbits are cop lan at and, exespt for the lunar orbit, a*re circular. 

In addition, the motion of tha earth about th* barycenter shall be 

neglected in each of the disturbing function. This ewtion esn be 

Included in the expansion by expanding the relevant vectors about the 

barycenter of th.i earth-moon system (PI unmet, 1240), but the expected 

1 

error introduced by neglecting it is of 0(jj^ * jp •*** ie smell 
compared to other approximations to be invoked U’«*r. Finally, the 
interaction of the moon on tha planetary partner will be ignored 
because of the planet's relatively larger mass and angular momentum. 
Since the interaction between the planet a rtf moon will tend to 
conserve angular momentum, any change In th« lunar angular momentum 
produced by tte planet will be balanced by an ecpial and opposite 
change in the planet's momsntua by th# moon. Therefor# the ratio of 
tha fractional, change let 




FIGURE 4.2.2 

vector diagran of the planet, earth and lunar petition* 
with respect to the sun. 


T, «p "» !BB Jn. 


The above ratio lndlcatei that the pertuifcatior of the planet's orbital 
elements by the mo<vt is <{ulte negligible. Therefore, we only need 
determine the first and iscond order contributions of the previously 
mentioned disturbing functions. Starting with the lunar •solar dis- 
turbing function, R (0 ■* 7 ), let's first axpanc the direct part A* 1 
(2.3.1) . 


A 



*■0 



The indiriot part cf Rf® *♦ will exactly cancel the l » 1 
term in A** 1 , will# the 1. • 0 tern does net contain any of the lunar 
elements and cert be dropped. The result Is (2.2.1) 
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The final expansion is of the function imi in terms of the 

mean anomaly M (2.2.0}. The complete expert# it n of R(© ♦ J ) i#* 




<4.2.2? 


The expansion of R(p ■* •) is even sixqpler, estcept that o’ 1 shall 
be expanded in terms oi Laplace coefficient* (2.2.15) i 




shore 


rote that all terms in tl.e above expansion, except j - 0 are of short 
period. 

Tne expansion of R(p -♦ 'P ) involves the same ^-ocedure, with one 
extra step. The result of the sequence of expansions just outlined lei 

R(p~> t)*Ap 

■ *5’>V w ..<o(V<-<v«> 

i v ID "*2** 

<a* T >|Q ~K~ 
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The angle 4^ is just the angle made by the vector D vlth respect to 
the reference frame vector in the orbital plane, ft (fig. 4.2.2). This 
angle la related to and through the following relations i 

ooo4 d « D*x « 

eir.^ * fi'fc - a^TT 1 (ein>^-^JinX p ) f 

from which can be d*riv«»d, after a binomial expansion i 

e-^D «(a (J D-^ ml y:(-) 8 (| 0 ) 0 (;?)e- i * X P 
* o*0 & 5 

The final rform of R tp * P ), after expanding D " # la 

1 ( LzbI * 


(4.2.4a) 


r(p->d)» I? W^lV^Up) 

J ---0 


where 


(4.2.4b) 

Prom the nyimnetry pro;?ertie* established in section ** 7 the above will 
collapse to a cosine series The exponential con be replaced 




is 
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~2\ » i ^-f — ) * ^ 


which vanishes for (ft - m) odd. Given l • 2, |a| is restricted to 

the values (0,2). For a synodic type resonance (q ■ 0), the t • 2 

term 1/ the above esq urn si on cociespor ", % to twice the cowraenBur rtbil ity 

ratio, or a double angle resonance. The disturbing function contains 
a 

a factor of for a given 4 ^ _ and would be raised to th<* 

a^ a^q^p 

eighth (|'-8) and fortieth v *40) powers for Mercury and Venue, 
respectively. The lowest order contribution to the single angle 
resonance is contained in the l • 3 tern of the expansion. A touch 
comparison of the coefficients of the single and dcubla angle terms 
is given by 



Calculation of th. sbove ratio indicates tint it is within an order 
of magnitude of unity for both planets. Therefore, both terms shall 
be explicitly calculated for thn synodic case. For the strongest 
possible e-type reionance, q - 1. As |ra| - G.2 for the £ - 2 term, 

the lowest order contribution for the e* ype corresponds to the 
single angle case. 

Tie lowest order contribution to * iven angle can be written 
down directly, usiig the explicit, lowett orde* -expansion of the 
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| 2 0 1 -e 


»— 

I 2 2-1 _*e 

l 

\ 2 2 0 1 


Definitions 




t T(-i k | 2l -r]i -f i 

I * 


-i , ,Jl+1 
2 <-y) l (-j 


» 2l 'V 


y»cos| f <r»sin^ # ^*ra-i +£p .-2k t 

and k is sunned ov«* all nonnegative factorials. 


_ _ j 

afD > *i_. *>«(*) 5 5 < 






TABLE 4,2.1 

Table of reliant Hansen** coefficients, x^ w * 

Included are the explicit definitions of the inclination 
function F, _{X). 
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Here (w,J> represents the complete set of conjugate action and angle 
variables. The functions 6lj , , etc., are the short periodic 

first order perturbations of the canonical variables and are obtained 
from the appropriate generating function through the relations (2.5.9) t 


f w) 

57 ? 


3 o( J , w) 
7vv 


(4.2.11) 


The generating function is ralat >d to the appropriate short 

periodic part (R^) of each individual disturbing function (2.5.2) 
acting on the noon. If each R # is formally expressed as a transcen* 
dantal ser.es in 4 (2.5.4a). 


xCV 14 *, Btc «. 


then the generating functions are (2.5.17) i 


3 R 




This allows us to write down 6 h and 4R in teens of the above 

o s 

formal expansions for 8 and R # , and determine which contributions 
are significant! 
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'VVe 

♦*- 3^ A *«ST.<^ + i^V 



recall that V is a function of L (cf, 2.2,5). Therefore the partial 
deriv«tive with respect to L explicitly acts on the frequency v in 
the above equations. 

There is the obvious requirement that each of the above angles 
V' *s" r or ha ® nonzero frequency, which is derived from the 
short periodic nature of R^. Because of the assumption of ooplanar, 
circular (except *) orbits, the only action variables which anter in 
the above are Lp , and L^» In the case of synodic resonance, 

the set can be restricted still further with the assumption that the 
lunar orbit is circular. The explicit angles which occur in the 
above are (4. 2. 2,3) 


+q B )Mj + ( j'-e 1 ) ( * 

V*V= n " l VV +<, -" K » +r( VV* C4.2.K.) 


Given an angle ^ ^ p, the second order coupling corresponding 
to $ b * ♦ ♦ ,* will be restricted to a single terra. That is, the 
indices ra", q", and J* m are fixed. The second order coupling corres- 
ponding to $ b , + ♦ „ will involve several terras, which can be 

restricted to the terms of lowest order by minimizing the number of 
a > 

factors of (— > and e. which occur in the expansion, consistent with 
a © 

the given angle. The ratio of the relative magnitude of second order 
coupling of perturbations of the lunar orbit with that of the earth's 
orbit in 6 K b is roughly 


(sec Old orders*) _.L.,n,,,2 !■«. 1 

(second orders*) » 


C4.2.15) 


Assuming tliat th. sum of Laplac. co.ff ici.nt* is approximately the 

same. The same ratio tends to hold for contributions from £h . 

o 

Given that the important contributions invol* -* frequence of 0(r j ), 
the approximate ratio is 



Thus it appears that the perturbations in the earth' pc 'it ( on 
give the principal second order contribution to a given angle in 
R. (p > ) » which suggests that the second order contributions be 


t-i jfl * s a** « -*<-** yrikh;^ ** *r ' *»’ -r**i 


***>«. ■ > M* 




( 4 . 2 . 18 ) 
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"type* . i 

2,-1 ,p 


^V' ? ,-1,P* ll 0 ^ ( t» )2 C4.2.21c) 

4x £.K (x> } co8< ^,-i,p 5 

- tn ” ♦o.i.P 8 

«\>o t x.p-^^5 ,a «*<W,p‘ • <»> 


-E« 2 

* 


The factor A la the second order coefficient In the expansion of 

OXX 

the zero order Hamiltonian (2. 7. 2d). Since the resonance has little 
effect on the planets involved, is approximately given by* 


sli*±ai . 


(4.2.23) 


The tidal torcrue term (A 2 _ tt) acting on $ is (2.9.10,17) t 
* oxx at 


The integer p is negative, implying that each of the bracketed terms, 
which are sums of Laplace coefficients, is positive. Comparing the 
above results with the direct part determined earlier, we see that 
both contributions have the same sign for each angle. 

The next step is to compete the lunar tidal torque with the 
maximum torque due to resonance interaction with the given planet. 

If the torque due to the action of the planet is conspicuously Larger 
than the tidal torque, the tentative conclusion is that the resonance 
is "tidally stable.** For this comparison, the second order equation 
of motion for the angle variable $ (2*9.17) is most useful. This 
equation is approximately 


4 A ox/ V in * * A oxx a 0. 


(4.2.22) 


~2 

The units of the coefficient A A. are tine , and the scale of this 

OXX 1 

coefficient is set by the factor 

i 


\ 


2 da 
oxx 3T 


h 0* 


(4.2.24) 


The oldest determination of is that of Spencer Jones. 


“ n T k -2 „23 2 

( — * ) - -21 • century ■ 1.1 * 10 J rad/sec (4.2.25) 

dt s.j. 

Recent investigations by Newton (1969), van Flandern (1970), 

Morrison (1971) and Oesterwinter and Cohen (1972) indicate that 
about twice Spencer Jcnes value is a more reasonable estimate. 

The function A^cos^ is a term in the disturbing function related 
to the resonance, and is the sum of the contributions from the 
direct and indirect parts. Explicitly, 

A 1 = A 1I +A 1D • <V B n>f-2 W * t4 - i - 26) 






Another interesting question is how lone ago the resonances were 
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established. A rough estimate of the tine, T that 4 

*wq#p J*#q>P 

approached zero and then reversed sign is given by (Allan, 1969). 



(4.2.27) 


dtvj , 

T is calculated using Spencer Jones value for ■ ■ ■ *■" ■■ (4.2.25) 

m.qp ^ ^ <Jt Q 

Since it appears that the more recent determination applies* only to 

the present value and rot what it may have been in the distant past. 

Table 2 contains the relevant parameters needed to calculate 

A A. „ Table 3 has the nwerical values of th*> relevant sums of 
oxx 1 

Laplace coaf ficients. Table 4 has the numerical evaluations of *1D 

and A lx , the estimated age of the particular resonance, and an 

evaluation of their tidal atability. 

The resonarce angle 4 has bee* constructed so that it 

m,q»p 

secularly decreases in the absence of the resonance. Therefore the 
analysis of capture and the stability developed in chapter three (3) 
directly applies. At the end of section 3.4, we mention that there 
are various types of stability, the most important here being the 
"tidal" and "adiabatic" stability. Recall that adiabatic stability 
is governed by whether the magnitude of A^ is increased or decreased 
by the long term tidal interaction Such that the amplitude of 
1AK a r Ion decreases with time. For a- type resonances, the leading 
factor is elx), where 


c 

«o 


«(ox+l 


q*.-k 


(4.2.28) 
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Planet 

Semi major 

Kean 

Maas 

Eccen- 

Incline- 

or 

AX*3 

motion n 

Ratio 

tricity 

tion to 

Satellite 

< A - U -> 10-W T 

MpAVicT 6 * 

Ecliptic 

Mercury ? 

0.3871 

8.396 

0.165 

0.206 

7°0* 

Venus ? 

0.7233 

3.2 30 

?.44 

O.OC 68 

3 O-. 3 , 

Earth Q 

1 . 00 " 

1.974 

3.04 

0.0168 


Moon ^ 

2.57 10 ~ 3 

28.37 

0.0376 

0.0519 

5°. ti* 


af * .00845^ m .I13n e (prorratfe). 

% »« . CNMOln^ . 0 5 3 6 r. # ( re t ro *-rad e ) . 

TABLE 4.2,2 

ORBITAL ELEMENTS FOR MOON, EARTH, VENUS AND MERCURY. 
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TA SX*E 4.2*4 

Numerical w tea of the max r. j» torque (jn ^g^ A^ A^ on the 

moon due to the given resonance. These values are to be 

compared with the present estimates of dn * which ranges 
-23 2 dt 

from 1 + 2,x l<f rad/sec , the higher value beinj the more 


recent determination 
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For small librations -~lt • t Q ) , Since i 8 positive 

(4.2.24) , e increases or decreases, depending on whether q is ■ « 
n^oative or positive, rsspectively. Thus . is adiabatic ally 
stable while ♦ _ . is not. The " tidal stability * of a given 

Q, 1, p 

resonance is determined by whether the tidal torque is greater 
(unstable) or less than (stable) the maximum amplitude of the 
resonance torque. 

we should observe that the adiabatic stability of the synodic 
resonances is not governed as much by whether b(x) increases or 
decreases, since this is a relatively small effect, but by the 
asytrmetry of the tidal torque. This a Bynum try arises from the fact 
that the torque is a rnpidly decreasing function of the planet- 
satellite distance. To lowest order, this asymmetry will add a i 
term to the right hand side of the second order equation of motion 
for 4 (2.9.17). This term is 



As the coefficient p is negative, this asymmetry implies a 

dissipative mechanism which tends to damp out oscillations. The 

results in Table 4.2.4 Indicate marginal tidal stability for some of 

the resonances. This qualitat'-ely agrees with Hipkin's results. 

The strongest adlabatically stable resonance involves Venu3 and the 

resonance variable , , Q . Mercury also has surprisingly strong 
2,-1, lr 

resonances. 




Unfortunately, there are seme serious flaws in either the 
approximations or th*i supposed effect of a given type resonance which 
drastically change the results so far obtained. Also, the question of 
the likelihood of capture into libration needs to be answered. But 
first, let’s examine the approximations more closely. 

We could discuss tl»e effect of second and higher order terms 
which have been neglected. Although these terms may be sizable, it 
it unlikely that they will critically change the order of magnitude 
of the coeif Solent A^. It appears k at the grossest approximation is 
the circularity of ilanetary orbit although on the surface it seems 
to be fairly reasonable, at least for Venus and the earth. After all, 
the present eccentricity of the ea/th ts * 0.017 while that of Venus 
is ** 0.007 which are both quite small. (Mercury's eccer riclty is 
*.2). But in the past ! Brower and Clemence, 1961b), these 
eccentricities have varied considerably due t^ the long period 
perturbations of one planet on another (Table 4.2.5), The toon 
itself is indirectly affected by such perturbations, especially by 
the long period fluctuation of the earth's eccentricity. These 
fluctuations have periods ranging from 50,000 years to approximately 
two million year£* The libration periods of all the planetary- lunar 
resonances are all roughly given by 

T period ~ 2irl0+ sec - ~60 f C00 : r rs. w , 2 . 30) 

* A l A oxx 

So the important question is, how ooes the fact that the planetary 
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(4 *. 34a' 

cosUetsinVi) £ J 0 (2te) + 2J 2 (2te)co3(2:/.) f 
8in(2etsinV) s 

b> 

Next, expend C*) S to second order in e. Also use the fact that ■ » 1 to 
simplify the coefficient*. The result is 

(“)° = l+(^) 2 -S0 cos f +(->£) 2 cos(f> . (4.2.35a) 

Next express the above in terms cf M, accurate to second order in et 

(^) S 58 l+(^) ’ — C! o(coj( tV' LM) . b) 

Multiplying the factors together and keepinn only lowest order 
secular terns, we lirve 


3 V 2f:e ' f(2 F> 2 “ l-(te) 2 +i(ie)2. t4 - 2 - 36 » 

In the IvLia - »*netery d i jturbing function, the minimum exponen* of 
th*» ta • a( is approximately |p|, where |p| >> 1. Therefore 

set t - p ii, the above equation for fc (e). Thus, to lowest 
order, we can deduce that taking the eccenU icities inu. account will 
multiply A 1 by the following factor; 


*D.p( e * ) X? Pt -(0~ (l-?(pe )( ) 2 )U_*(oe ) 2 ). 


(4 .2.37) 


The data in Table 5 is taken from an article by Brouwer and van 
Woerkom (195^), a.td gives the maximum variation of the eccentricity 
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TOStE 4.2.3 


ihe secular elements of the planet* an determined by J*rouw*r 
an 4 Van ttoerfcom (1930) using th* l*tple:e»Ingrange approximation. 


due to the Interplanetary perturbation* just discussed (also in 
Brouwer and Clemence , 1961b) using the Laplaoe-Leguange approximation. 
Zt appeirs that the fluctuations in bUc) for Venus must be of at 
least an order of magnitude and may actually reverse the sign of b{x). 
Mercury, on thn othnr hand, fluctuate* only by a factor of two or 
three for the single angle synodic raer nance and the e-type resonance. 

Thle it the potential associated with the taro* 

eccentricity approx Una t ion is now split among several "side-band* 
frequencies which differ by the motions of the perihelion* of the 
eaith and the disturbing planet. The motions of the inner pi me tar y 
perihelion! are of order 3 * l© 3 see. of ere per Century, or 

3t plnn.t. * * ' »° J V«»ntury v (4.2. Ml 

Compere th l* with the maximum valun of 4 that the resonance can 
withstand without being disrupted! 


" Wi * 3 


* 10 


12 


(4.2.J9) 


Obviously th# side band freqianciis differ from the libration 
frequency by soma mall multiple and should smart strong accelerations 
the resonance variable. It is possible that the moon could be 
trapped not at the single frequency, but among the lever* 1 aide bar* 
frequencies, lut such a motion, a. it 'ely, appear i to be highly 

unstable to tidal disruption. The impotent point Is thst an 
eccentric planetary orbit effacUvely destroys thn phase relationship 


4 U< !* S) tlii' H 'i 4.^.1- 
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between the planet end noon at 3 injunction because of thu Ugh 
ccemensurability ratio. This is specially true for Vamia since 
c ?o. 

Tharo it yet enothnr long period olfect associated with the 

varying planetary iccentrlclty. A varying eccentricity implies a 

tire varying notior of the planetary perihelion or a Ion*; period 

•■secular* acceleration. Thia acceleration is of 0{sec. cf arc/ 

2 

tertury ), and because the associated periods are up to several 
multiples of tne 1J hr at;, on period, the> would appear over a fm* 
oscillations as a necular acceleration in the equations governing the 
side band frequencies. A wore significant at feet arises from the 
feet that these long periol variations (primarily that ol the earth) 
induce a correspond ing effect on the lunar eccentricity and 
perihelion which wr uld ten! to disri.pt the eccentricity-depenlent 
resonsneis. 

The next probi en is to estimate the probability p that sny 

c 

one of these resonance variables made the transition from positive 
rotation into librition. Admittedly, the theory developed in 
chapter three ajrpl ed to the esse where the tidsl torque :,* mnall 
cor, pared to the pa, irouo pendulum torque. K^r did we not consider 
how long periodic lucturations ,.n the coefficient bb<) of th* 
resonance term aff -cted capture. Therefore, any numbers calculated 
should be lndiceti/e of t) e gere.tal magnitude of P^. For an e-type 
resonanc • i|k| - 1 , is apt reximatsly given by (3.2.12) 
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For an e-type resonance where b(x) is propcrtional to e, 0 tend* to 
scale like e~ 3 . Therefore, if the lunar e< centric ity were ten tires 
smaller at transition than it is at present., the probability for 
capture would be increased to h* 3%. 

In the case of the synodic resonance, capture ap*tars to be 
mainly due to the ^-dependent tent asso<.ia*;td with the tidal torque. 
An estimate for P c can be obtained from (3 . 1 • 14 , J 5) t 
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P C , i£i^±3iLA^L_p. 


(4.2.44) 


Tor a double angle synodic resonance, (m + q) • 3, whU* for a single 
angl# resonance * ♦ q) • 1, The strongest possible synodic resonance 
involves Venus v;. fch 

^s.c.n- W ? -VV3. 

Tne probability .for capture is ^ 0.5 * lo“ 5 , and unlike tbs previous 
case, thare appears to be no way to substantially increase It. The 
conclusion to be drawn is thnt capture la a hichly unlikely event ir 
only tidal force*. are Involved. 

The l. ng t«nn variation* in the pendulum torque associated with 
a given resonanc. » may cause temporary capture. Since the variation 
in the coefficient b(x) is a significant fraction of its mean value, 
v* should expert that thr> probability for temporary capture is much 
larger than tho numl>ers calculated for pnnanent capture. But the 
maxiim-m tine that the 1 inrat ion can last is of 0(10 5 yrs) # the time 
scile fc» Lrvj set vy the ;?*rio<ii$ associated with the so-called secular 
variations of the planota. There is one final argument to be 
leveled against our earlier result*. It is that only the synodic 
type resonances tend to "lock" the moon at a fixed radius without 
appreciably aff rcting tho other orbital parameters. Pecall that the 
sta: le ^type i .‘jso’tanc'is increase e if tnere is a tidal terque. The 


mean notion, and therefore the semimajor ax; a, is fixed if the moon 
is trapped In one of these resonances, but the mean radius actually 
tends to increase. The short period average of r to lovost order m 

e is 



The present value of e^ puts an uppor hound on tho lifeline of ar 
e-tyj>e resonance. This fraction is qualitatively small compared to 
the change in a^ due to the tidal acceleration over the lifetime of 
the earth-moon systisn. Thus, only the synr.dic type issonanra 
effectively trap the moon at a fixed rad i*r. Put tha synodic 
resonarces are, at beat, marginally stable to tidal disruption, .ire 
sublet to long period disruptive torques due to the- "secular" 
interplanetary perturbations and the associated "secular* change 
the lunar mean motion, and have a very low probibilitv of capers. 

We ire forced to conclude that It. is unlikily that the noon was ever 
trapped in an orbit-orbit resonance with a planet. 

There is one last question to ask. C>uld the orbital elements 
have been signif icantly changed oi‘ her by passage thioi gh t'-ir 
resonance or by temporary capture? Hie change in tho ao«n value of 
the momentum variable due to passage through resonance is 
approximately given by (1.2.47) i 


Ax ^ ®|b| 


1/2 m 





273 

Keive airuiy s*en £i*t 0 is a vary snail parameter for sac h of tha 
ur&r resonance'..* For an e-iype resonance , Ax is roughly proportional 
to tha fraction l change in the eccentricity. Vhe above effect cannot 
** very 1 arge if lean cy<aTo ’/ary small before transitio \ , 

Th> phenor i non of temporary capture for a., e-type resonance 
covltl lead to .» much Wrger change in tho eccentricity, although 
still cnly a trill fraction of its present value. This can be 
deno nst rat <k 3 by comparing toe estimated lifetime of the resonance 

(•v io 5 yrs. ) w:.vh the tune necessary to havo tldally changed the 

2 

fractlc nal rsdi.ia of the wcon ky an amount aqual to l/2ep (4.2.45). 

Thin t me is approximately 

1 7 

^ ?T.!0 veird . 



At nos;. ta*np<rrary capture may have change the lunar eccentricity 
by as iucf as 1%. 

w* should note that the inslination-ti'pe resonances have not 
**en nwiticnel since their strengths are an order of magnitude 
snails r than tiat of tn* e- types. This f o' lows from the observation 
that all : or uixed-I type resonance variables hav« a ccsfflcient 
2 

fc(x) f roportional to X or IX*. 


4.3 OH THE TKW52 -SATELLITE ftESCKKMCE OF JUJ ITER 

The theory developed in the first three chapter? has n*t with 
varying success whsn applied to specific examples. The priiciral 
problem is that whereas the Hamiltonian invoked is o i*~dir*!ifar ra ' , 
thn real world is not. There always seemed to he a ronplJcstirr i it 
wan side-band frequencies in the car- of toe lun.tr a id planetary 
resonance. In the example of the two-body resonance, it was tl* 
possibility that sometime in their evolution two resonance variables 
associated with the same commensurability might overlap and destroy 
tho simple one-dimensional description of the phenomena. Vhather 
«u;h complex phenomena will yluld their a - crets *o readily aa the 
Hanllumian we derived with the analytical tor.ls developed to 
describe transition is best left to future investigation. 

There is a mere specific example of a satellite esomnee ir. 
which the necessity of a rmilti-resonance-variable t! eory to explsln 
Ijo th capture and its present libration amplitide smns unavoidable. 
This is the three hod;, f -*1 race relation satisfied by she tnree 
i.rer Galilean satellited of Jupiter. (1 J3 , - 3 \jII 
(Stagihara, 1972). The following observations about t*is relation 
sre especially interesting. 1) No libration amplit ids has benn 
observed, 2) the satellites also nearly *atinfy a Sil comoT*».s*u:abiH< y 
between the inner pair (JJ and JIX) and the outtr pair (Jil and Jill t 
3) at the present tine, a tidal torque which acts rrincipally or. 
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tha irx.er satellite* crives than tower i thin coon* nsur ability; and 
4) th< objfnH eccenlricitiot ax* snail a.-*2 variable and th* naif as 
urn wiLhli ai order cf mgnitude cf one author. 

:n imar th*'.ry we found that wy^xJU' frequencies could be 

found in the '-xpauslar of the Jlstuzbinn function. It happened 
herause the >ylncipal disturbing fut ction R{p ♦ J. > contained 
-ix?lic: Itly- tm- orbital elements of the *m i t v . as wail as of th* 
perturbing planet and th* moon. Thxs occurs only because thsrs 
exist two distinct "trr Unary" bodies in tho problem! th* sun for the 
planet i end ti* earth for th* moon. In vUpiter's thre*-sa tell its 
resonance tier* exists one coarwon primary. Any inter-satellite 
gravitational interaction contain* the orbital elements of just two 
of them. Sjrt* th* .'supposed tesonance Involves jxpliciUy the Man 
longitudes of three bodies, this angle munt o^cur as a second order 
tern in a pci turtat*or. expansion. Thin inplles that pertur bat ions 
in the mean *onoitu>-e or senixajor axle off any on* of these 
ae tel 1 ites 1‘ of 0('t*) 1//2 • j i). a typical me** ratio i* 10~*. We 
found that li the case of lunar synodic resonances, both capture *n£ 
damping were extreme . y slov , fcuing gcverivx? by th* cagnltude of th* 
fluctuations it tlm wean motion. Sc what is going on in the 
Jupiter sato’lite cane? Did we miss something in th* dsvelopTwn^ 
*hich night adically affect the result*? 

Reral- that in the generating function f^cm which th* start 
pencd.ic p< zturbntions wern obtained, th*r* occrurrsd ir th* 
denar in* tor uf each tern a factor nearly uqual to the frequSTicy of 
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the corresponding cosine argument J2.5.21). it so happens that 
Sou ill art fee* H&glhar*, J972) has developed an analytic theory along 
the sane lines as outlined in chapter three, and he finds that tli* 
dominant contributions involve th* coupling of e-type anal* variables 
which nearly satisfy tie near cwroirrorajllitim arong the inr,«r 
and outer pairs of satellites. W* rhould also poir.t out that the 
synodic frequency occurs as tha difference of two t— type frequencies. 
Specifically, they are 

la * + 

S “ A II" ?A III + * > II S 

*a " *b ■ *1 * 3X II " 2X IX* 

Th* frequency occurring in the second order coupling of these 
dominant term* is 

*>“ n l" 2n ir ( ^ 1I) aec* 

Me should note that th* tru* frequency associated with the related 
e-typa *nqle variable differs from the above by an amount <j~) ecoei) 
and is a natural result of tha pertuebatior expansion, trh* 
fluctuations . which nay occur in v due tp th* refinance appear to be 
restricted to the zoean notions. Another possibility i# that i? the 
* cos# term did contribute a large retrograde notion to 2, its 
absence in v may indicate a breakdown in the perturbation mothod. 

The differences of the moan motions in the tarn (n x - 2n JX > 


* a ^ i i ii;t < 





277 

nj « ?r;?48g/U y* 
n,j ,! 101?174/*?i»y« 
n--:>n n * O?7-'0/1oy « >7C?l/joar. 

one of the possible implication* of previous remarks is that the 
interaction of on# or ware e-type variables with the synodic \ariable 
was crucial in both capture and subsequent damping of the resonance. 
II thesu variables were Lignif .cant, then the impressed retro<rrade 
motion of, say, mu«t have been very largo in the past (greater 

than 2 * 2* T 0 e .1 /year) . Fro® (4.1.3) we find that if wers is large 
as U 2 , the impressed retrograde motion of £ is still cf 0<djgtees/ 
day). My quest is that the critical evolution which « -{plaints tike 
present anal l libraUons appears to be tied to an irt ruction 
involving synodic ar.d e-type angle variables. How those variables 
came to overlap can be exp lair ed by either of the following t vo 
scenarios! 

1) Tier a it ion Into the three-body sync -die resonance occurr-*J 

f ir*t. The system then evolved through the e-type resonance. 
Tie 3© variables* complex interaction led to a rapid damping 
of the amplitude of 1 iteration of the synodic variable. 

2) First, a pair of satellites established an e-type resonance. 
Subsequent y the system evolved toward the synodic ccsnoen- 
surability. Somehow the original e-type resonance wee 
disrupted whil* allowing capture into the synodic to occur. 






. . -31 e titbit* At:- 
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Because of capture considerations the second seems the most plausible. 
Whether the above remarks have any relevance must await a more 
rigorous examination of this type of problem. 
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be one of the constants of notion. Also, the tine dependent part of 
Sta^q^t) can be chosen such that H' is identically sero. if we 
write 

Sfci.qift) * ♦ S’Vj.qi.t), 

then o is equal to H , or 
1 O 

"J m -h • 


if we demand that the new generating function sMa^q^t) is inde- 
pendent of t. 

33 

The next step is to substitute p A for — - in lA.2). The result 

is 


, if, # 2 


■ht&t o 


(A. 7) 


The only form of SMo^q^) for which (A. 7;- is separable is 


s*. S 2 (») * SjW ♦ S t (r). «-8) 

Since 4 does not explicitly appear in (A*?) set equal to 

3 2 • *}*• (A. 9) 

Separate the remaining terms in (A.7) into those which depend on 6 

2 

and those which depend or. r» and equate each to a constant a^. The 
solution for S is 

s * r 1 ir * 


The constants r , 9 are at our disposal. Chooie the valuos 
o o 


r 0 ■ r c®pcricec * n(l - o) l 9 0 * 0 


(A. 11 ) 


are defined by relation {A. lb): 


(A. 12) 


'fait* 


Beginning with (A. 12a), we know that the integral most be equal 
to 6^ + t and that 6^ is a constant. This suggests that the 
integration variable be changed from r to M since M is a li.tai 
function of t. Using the relations 
r » e(l - e C0B K) 

and 

* 

Jr. . e«? Iirv e, 
dr. r 

the integral becomes 

f r , r M _ ea’rirE d*. (A.X3) 

J> fsoL, . **- ((Ma-ti-e 2 ;../;. - 

r <* 1 £ M * r t* 
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FIGURE A.l SPHERICAL TRIANGLE 
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@2 8 u) » S3 - n 


(A. 17) 


Thr» integral in the B^ equation can be reduced by a similar 
tra '.sf ormation to 


*JL 


/U 

I 

'9 



: tr ' 1 f to* 


t in. u) 


The transformation is defined by the relations 


a) (3^ « 







Again, let R ■ 0 to determine the new set. 
we can demand 


(A. 21) 


Since H* • 0 if R * 0, 


But this integral equals i - 6^ by (A. 12c). Therefore, 
tnr.( 4 ‘ - S 3 ) ■ cob i tan u. 

Again ve appeal to the relations derivable from Fig. A.lr we find 

As* -ft. (A. 18) 


Collecting the results, we have 


*2 » v 4a(T - e ? ) 


e cob I^n(l - ©2) 


^2 * S3 - ti 

^3 ~ 


tA.19) 


The equation of motion for the set of orbital elements {a, s, I, c, ft, ft) 
can be derived from the above and (A. 3) as a purely algebraic exercise. 

The above set is not Uie most useful 'or our purposes. Make 
another H - J transformation on to a naw ®*t ^i ,M i^ an<3 

demand that the new Hamiltonian satisfy 


V 


* ^ » / 


£* - -*jt * #-(„ lf m) 


(A. 22) 


where the new generating function is explicitly independent 

of the time. We have some freedom in choosing and shall 

demand that it be equal to the identity transformation in the old 
variables { 0 ^,$^ 




Now we can write down the original SMo^rB^) and use (A. 15) to 
eliminate t. The result is 


SLk 

n 


M « f 3-Oi ♦ 01 j u • 


(A. 23) 


There is cbvloualy quite a bit of freedom in choosing the new 
variables. One choice for is 
6 l> z 8 H 

ut z • d - n (A. 24a) 

^ 

For which the conjugate action variables are 


H* 1 * H* 


K"*l 


(A. 20) 
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J 1 s - ^ - /H2. 

J 2 z 0*2 9 >/***(! - ® 2 ) 

^3 * ^3 3 cos I^afl - e 2 ) 


(A. 24b) 


This set constitutes the well -known Delaunay system of elements. 
Several similar sets of conjugate variables can be obtained by 
rearrangement of the angle variables in {A. 23). A modified set used 
this paper is 


L- y A -M* a 

r *(fi - e 2 *- > a 

i - (1 - cos i) J/4 a(l - e 2 ) ; D . 


(A. 25) 




